Pee) Woligang Paull 


si Oe edited by C. P. nz 
Volume 5 ewer ye 
Wave Mechanics Victor F . -Weisskopt 


- Pauli Lectures on Physics: 
Wave Mechanics 


Pauli Lectures on Physics 


1. Electrodynamics 

2. Optics and the Theory of Electrons 

3. Thermodynamics and the Kinetic Theory of Gases 
4, Statistical Mechanics 

5. Wave Mechanics 


6. Selected Topics in Field Quantization 


PLY 

197 

¥,S Pauli Lectures on Physics: 
Volume 5. Wave Mechanics 


Wolfgang Pauli 


Edited by Charles P. Enz 
Translated by H. R. Lewis and S. Margulies 


Foreword by Victor F. Weisskopf 


Library 
Western Wyoming Community Collegé 


The MIT Press 
Cambridge, Massachusetts, and London, England 


Copyright © 1973 by 
The Massachusetts Institute of Technology 


All rights reserved. No part of this book may be reproduced in any form 
or by any means, electronic or mechanical, including photocopying, 
recording, or by any information storage and retrieval system, without 
permission in writing from the publisher. 


First MIT Press paperback edition, 1977 
Second printing, 1978 
Third printing, 1979 


ISBN 0 262 66037 7 (paperback) 
Library of Congress catalog card number: 72-7805 


oh 


Contents 


Foreword by Victor F. Weisskopf 


Preface by the Editor 


Preface by the Students 


Introduction 


Se wh 


aw 


to 


won 


We 
12. 


Wave Functions of Force-Free Particles 


. Association of waves with particles 

. The wave function and wave equation 

. The uncertainty principle 

. Wave packets and the mechanics of point particles. 


Probability density 


. Measuring arrangements. Discussion of examples 
. Classical statistics and quantum statistics 


. Description of a Particle in a Box and in Free Space 


. One particle in a box. The equation of continuity 

. Normalization in the continuum. The Dirac 6-function 
. The completeness relation. Expansion theorem 

. Initial-value problem and the fundamental solution 


. Particle in a Field of Force 


The Hamiltonian operator 
Hermitian operators 


XI 


xiii 


XV 


23 


23 
27 
31 
33 


38 


38 
40 


vi 


it. 


14. 


26. 
27. 


29. 


30. 


CONTENTS 


Expectation values and the classical equation of motion. 
Commutation relations (commutators) 


. More than One Particle 


More than one particle 


. Eigenvalue Problems. Functions of Mathematical 


Physics 


. The linear harmonic oscillator. Hermite polynomials 

. Matrix calculus illustrated with the linear harmonic oscillator 
. The harmonic oscillator in a plane. Degeneracy 

. The hydrogen atom 


. Collision Processes 


. Asymptotic solution of the scattering problem 

. The scattering cross section. The Rutherford scattering formula 
. Solution of the force-free wave equation 

. Expansion of a plane wave in Legendre polynomials 

. Solution of the Schrodinger equation with an arbitrary 


central potential 


. The Born approximation 
. Seattering of low-energy particles 


. Approximate Methods for Solving the Wave Equation 


Eigenvalue problem of a particle in a uniform field 


The WKB method 


. Matrices and Operators. Perturbation Theory 


. General relationship between matrices and operators. 


Transformation theory 

General formalism of perturbation theory in the matrix 
representation 

Time-dependent perturbation 


Angular Momentum and Spin 


43 


51 


51 


55 


59 
63 
72 
88 


107 


108 
110 
112 
lis 
116 


120 
123 


126 


126 
132 


138 


138 


143 
147 


152 


CONTENTS vii 


31. General commutation relations 152 
32. Matrix elements of the angular momentum 154 
33. Spin 156 
34. Spinors and space rotations 159 
10. Identical Particles with Spin 165 
35. Symmetry classes 165 
36. The exclusion principle 167 
37. The helium atom 169 
38. Collision of two identical particles: Mott’s theory 173 
39. The statistics of nuclear spins 175 
Exercises Liz 
40. Fundamental solution for interval Wetes 
41. Bound states and tunnel effect 179 
42. Kronig-Penney potential 180 
43. Spherical harmonics 180 
44, Fundamental solution for harmonic oscillator 182 
45. Angular momentum 184 
46. Partial waves 185 
47. The symmetrical top 186 
Bibliography 191 
Appendix. Comments by the Editor 193 


Index 199 


Digitized by the Internet Archive 
in 2023 with funding from 
Kahle/Austin Foundation 


https://archive.org/details/paulilecturesonpO0005paul 


Foreword 


It is often said that scientific texts quickly become obsolete. 
Why are the Pauli lectures brought to the public today, when 
some of them were given as long as twenty years ago? The 
reason is simple: Pauli’s way of presenting physics is never 
out of date. His famous article on the foundations of quantum 
mechanics appeared in 1933 in the German encyclopedia 
Handbuch der Physik. Twenty-five years later it reappeared 
practically unchanged in a new edition, whereas most other 
contributions to this encyclopedia had to be completely re- 
written. The reason for this remarkable fact lies in Pauli’s 
style, which is commensurate to the greatness of its subject 
in its clarity and impact. Style in scientific writing is a quality 
that today is on the point of vanishing. The pressure of fast 
publication is so great that people rush into print with hur- 
riedly written papers and books that show little concern for 
careful formulation of ideas. Mathematical and instrumental 
techniques have become complicated and difficult; today most 
of the effort of writing and learning is devoted to the acquisi- 
tion of these techniques instead of insight into important 
concepts. Essential ideas of physics are often lost in the dense 
forest of mathematical reasoning. This situation need not be 
so. Pauli’s lectures show how physical ideas can be presented 
clearly and in good mathematical form, without being hidden 
in formalistic expertise. 

Pauli was not an accomplished lecturer in the technical sense 
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of the word. It was often difficult to follow his courses. But 
when the sequence of his thoughts and the structure of his 
logic become apparent, the attentive follower is left with a new 
and deeper knowledge of essential concepts and with a clearer 
insight into the splendid architecture of reason, which is theo- 
retical physics. The value of the lecture notes is not diminished 
by the fact that they were written not by him but by some of 
his collaborators. They bear the mark of the master in their 
conceptual structure and their mathematical rigidity. Only 
here and there does one miss words and comments of the mas- 
ter. Neither does one notice the passing of time in his lectures, 
with the sole exception of the lectures on field quantization, 
in which some concepts are formulated in a way that may 
appear old-fashioned to some today. But even these lectures 
should be of use to modern students because of their compact- 
ness and their direct approach to the central problems, 

May this volume serve as an example of how the concepts of 
theoretical physics were conceived and taught by one of the 
great men who created them. 


Victor F. Weisskopf 


‘ 
Cambridge, Massachusetts 


Preface 


It is astonishing that Pauli had given a complete version of 
this course “Wellenmechanik” only once, during the winter 
term of 1956-57. This is due to the fact that although quantum 
mechanics was part of the physics diploma curriculum of ETH 
at Zurich, a regular course on the subject did not exist for a 
Jong time, a surprising fact considering that the author of the 
famous article “Die allgemeinen Prinzipien der Wellenme- 
chanik” (Handbuch der Physik, Band 24/1, Springer, Berlin, 
1933) joined the ETH faculty in 1928. 

While much of the spirit of the “Handbuchartikel” is re- 
flected in this course, Pauli took the opportunity to put into it 
much of his extended knowledge of that 19th-century mathe- 
matics symbolized by Whittaker and Watson’s A Course of 
Modern Analysis, for which he shared a liking with his teacher 
Arnold Sommerfeld. This is one of the features of the present 
volume not found in the same detail and rigor in other books 
on quantum mechanics, which makes it a worthwhile text also 
for today’s students. 

But apart from this more technical aspect it is of interest to 
see how the most critical among the founders of quantum 
theory actually taught his subject. As mentioned in the preface 
by the students he did it mainly in mathematical language with 
few hut penetrating comments. And he chose the topics of this 
course according to their conceptual and historical interest. 
Examples are the probabilistic nature of quantum theory, the 
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concept of spin, the problem of identical particles, and the rela- 
tion of the statistics of rotational states of diatomic molecules 
to nuclear spin. For this last reason and because of Pauli’s deep 
involvement in the development of quantum mechanics, I have 
tried in the appendix to comment on some of the fascinating 
highlights of its history. 

Pauli’s second delivery of the course starting in October 1958 
was tragically interrupted by his death in December of that 
year. He had treated the first 16 sections, as well as the hydro- 
gen atom in polar coordinates, when I was asked to take over 
the course. As mentioned in the preface by the students, Pauli 
had intended to correct their notes taken from the first de- 
livery; he therefore made corrections only for the part men- 
tioned. Some of the responsibility for the remainder of the 
course and, in particular, for the wording of the problem sec- 
tions therefore falls upon me. 

The carefully prepared notes by Herlach and Knoepfel have 
allowed a fairly smooth editing job of this English translation, 
in which the work done by the translators was of great help. 


Charles P. Enz 


Geneva, 27 October 197] 


Preface by the Students 


Professor Pauli participated personally in the development 
of wave mechanics. For this reason he was able to deliver this 
lecture course of four hours per week without notes, and he did 
it with particular pleasure, as a complete master of the subject. 

At Professor Pauli’s request, we took notes in his lectures 
during the winter semester of 1956-1957 and worked those 
notes into their present form. He wanted to correct the manu- 
script by delivering his lectures from the notes during the 
winter semester of 1958-.1959. Unfortunately, it was not pos- 
sible for him to complete that objective. However, the very 
small number of corrections in the part which he looked 
through (up to the hydrogen atom) encourages us to publish 
the entire set of notes, and we hope that this will be as our 
revered teacher would have wished. 

In these lectures he wanted particularly to treat fully the 
mathematical foundations of the theory. This treatment of 
wave mechanics differs in yet another way from that found in 
the usual textbooks. 

Professor Pauli was fond of expressing himself with formulas 
and only a few words. Once he said, “One should not write so 
much.” In working through the notes we have taken pains to 
preserve the originality and the special style of Professor 
Pauli’s lectures. We have also permitted ourselves to smuggle 


some of his characteristic remarks into the text. 
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xiv - PREFACE BY THE STUDENTS 


We are especially grateful to Dr. C. Enz for correcting the 


second part. 


F. Herlach 
H. E. Knoepfel 


Ziirich 
December 1958 


Introduction 


The discovery of the quantum of action by Planck in 1900 
initiated the development of wave mechanics. In the begin- 
ning there were many problems associated with the theory. It 
took all of 25 years until wave mechanics was formulated as a 
self-consistent theory. Of course, this could only be achieved by 
giving up a certain concreteness; in particular, a system of sev- 
eral particles can no longer be described with concrete waves. 
Very soon after the first papers by de Broglie, Heisenberg, and 
Schrodinger in 1927, the fundamental principles of the theory 
were logically complete. Since then, the theory has proved use- 
ful in numerous fields of physics, and it has been corroborated 
repeatedly by experiment. 
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Chapter 1. Wave Functions 
of Foree-Free Particles! 


1. ASSOCIATION OF WAVES WITH PARTICLES 


A particle with energy H and momentum p can be asso- 
ciated with a wave, Aexp[i(k-x—wt)], as follows (k= 
(22/A)n is the wave vector and n is the wave normal). 

For light quanta there are the relations 


E=ho, p= hike, pt] 


which are relativistically invariant. In addition, there are 
the relations 


6) E 
| ole ae a as : 
le 
ot ws [1.2] 
gS ra 


EB ——————s 
= = Vp? + mc (m = rest mass), [1.3] 


mec; = Mv 
~ V1 — o/c o/c?” a V1— 22/02’ 


where [1.3] follows from [1.4]. Because of the general for- 


[1.4] 


1 Remarks: 
1. In these lectures we use the symbol h to denote the quantity 1.05 x 10-* 
joule-sec. In the older literature this quantity was usually denoted by ’. 
2. Limits of integration —co and + co will usually be omitted. 
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mula from mechanics, 
ny 
dieu dp , or (in components) Cis ig [1.5] 


we can also derive [1.4] from [1.3]: 


0H Pt Pi 
i —iG =C¢ 
‘pi V p? + me? EH 
v2 C22 v Cy m?2e4 
oom? lati ge aoe [6] 


Equations [1.5] and [1.6] can also be combined to give the 
relation 


(1/c?)E dH = p-dp. 


The idea of de Broglie was that [1.1] should also be valid 
for a material particle, in which case [1.2] must be re- 
placed by 


a) i= w? mc? 
Wag > = oe (1.7] 


h 


which follows from [1.1] and [1.3]. For light (m=0) we 
again have [1.2]. 
Substituting [1.1] into [1.5] yields 


= ve ; {1.8] 


that is, particle velocity = group velocity of the associated 
wave. From [1.1] and [1.6] follows |v| =¢%(k/w); accord- 
ingly, for the phase velocity u, 

[1.9] 


@ Gc? 
= —— 
k ov 


Since v<e, therefore u>e. 
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2. THE WAVE FUNCTION AND WAVE EQUATION 
a. Superposition of plane waves. Wave packets 


The most general wave packet is of the form? 


w(x, o=|ffaus, ky, ks) exp [i(kex—wt)]dk,dk,dk, [2.1] 


where w is now given by [1.7]. 
This wave function yw satisfies the relativistic scalar wave 
equation, 


1 m2e 
(v—4 >? sa) we j= Nee. t) [2.2] . 


as can be seen by substituting [2.1] into [2.2]. Because 
of [1.7], the wave equation is then satisfied identically. 
Therefore, we can also write 


0 Pi) 
ae on STAY or 
— iki, 7 iw 2.3 


These correspondences, together with [1.1], yield the im- 
portant connections between the differential operators of 
time and space and the classical quantities p and EH: 


> 3 
agi ae, Via) OF 9A 
vm — aa 2] 


These form a translational key between the quantities of 
classical mechanics and the operators of wave mechanics. 


b. Transition to the nonrelativistic approximation 


In mechanics, for v<ce (coc) and p<me, we have 


2 mc? 


= V pr mii me(1+ 5 P ‘-. : 
es 
= ={me ia e+... ; [2.5] 


* See W. Pautt, Lectures in Physics: Optics and the Theory of Electrons (M.1.T. 
Press, Cambridge, Mass., 1972). 
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From [1.7] we also obtain 


EH me? h 
= — = + —_}24 .., 2.6 
=o h T 2m Bex: 2-3 
(where E=me?+E#,,,, Ey, =p?/2m). We define 
h 
wo’ = oe [2.7] 
in order that 
me? 
CS +o’, [2.8] 


and 


wp! (x, t) =| J fave exp [i(kex — w't)] dk, [2.9] 
with which 


2 
p(x, t) = exp|— —— tlp'(x,t) [2.10] 
Substituted into [2.2], this gives 
mo? im dp 1 dy’ mc? 
ae] f aman Phen, eas t 
Vegi ise + ee at ene 


and the nonrelativistic wave equation, 


,, 2m dp’ 1 dty’ 

NO ese ot ain (2.11] 
follows. Aside from the imaginary coefficient, it corresponds 
to the equation for heat conduction. The imaginary coef- 
ficient assures that there is no special direction in time; 
[2.11] is invariant under the transformation t>—t, yp’ y’*, 
whereby y*y remains unchanged.? 

From now on we shall always calculate with the primed 
quantities which we have introduced here; however, for 
the sake of simplicity, the primes will be left off. The 
quantities w and w’ only differ by a constant; however, 
this is not an essential difference since only frequency dif- 
ferences are ever of importance in wave mechanics. 


* We shall see later that the physically measurable quantity is not the wave 
fuction ¥, but only the probability density v*y. 
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3. THE UNCERTAINTY PRINCIPLE 


The kinematics of waves does not allow the simultaneous 
specification of the location and the exact wavelength of 
a wave. Indeed, one can only speak of the location of a 
wave in the case of a spatially localized wave packet. The 
number of different wavelengths contained in the Fourier 
spectrum increases as the wave packet becomes more local- 
ized. A relation of the form Ak;Ax;> constant seems rea- 
sonable, and we now want to derive this relation quanti- 
tatively. For the sake of simplicity, we carry out the 
calculation in only one dimension; the generalization to 
three dimensions is immediate. 

We consider a wave packet, [2.1], at a fixed time t, which 
we take to be t=0. As is required by the Fourier inte- 
gral transformation, [2.1] can be written symmetrically 
in # and k: 


all : : 
p(x) = am [40 exp [tka]dk , [3.1] 
ey Wo w(x) exp[—ika]daz. [3.2] 


Because of the symmetry of these formulas, all equations 
remain correct when the following substitutions are made: 


yp A @& Kk + 
: [3.3] 
A ee ee. 
Further, the famous formula of Parseval, 
v= [yrtoryle)ar=|a*aya@ar, [3.4] 


is valid. 
a. Average values of functions and operators. Normalization 


For a normalized wave packet the normalization integral 
N equals unity by definition. A plane wave of infinite 
extent gives Noo and, therefore, it cannot be normalized. 
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We define the average value of a function F to equal the 
following quantities: 


= F(ax)p*(x)p(a) dx 


7 J ¥(k) A*(k) A(k) dk 
P= eeestayda 


E(k) = " [A(R A(K) dk ov«£ 


The quantities y*y and A*A in these formulas have the 
meaning of a density. Later on, this interpretation will be 
better justified. 

From now on we always assume that p(x) and A(k) are 
normalized: 


F(x) = | F(a) y*(w) (a) de , [3.5] 


F(k)= | F(k) A*(k) A(k) dk. [3.6] 


By means of the translational key [3.3], we obtain the 
operator F(—i(d/d«#)) from the function F(k). How to 
form the average value of such an operator, that is, on 
what functions to allow it to operate when the average 
value is calculated, can only be decided by more detailed 
considerations. The result is 


F(k) = [re ) Le (- 2): (a) 


F(a y= [ae B(+i3)-a00| dk. [3.8] 


Because of [3.3], [3.8] is a result of [3.7[, and vice versa. 
We prove [3.8] for the case F=polynomial: 
lL. Pee 


da , [3.7] 


i! 
pt oe ce 2 : ( + Te 
a | ay* da Jaw exp [ikx]dk 


Pak ; ae P 

= seve faw (-i ex [tka }) dk 
A 

= sale i dk: fv exp [tka ]dx 


=| A*(k) (i a 7) dk 
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2. F(x)=a". The proof here is analogous but involves n 
partial integrations. Thus, formula [3.8] is proved for a 
polynomial. It can also be proved without difficulty for 
an entire function (using the Fourier integral theorem). 


b. The uncertainty relation 


We define 
(Sx)? = (a —%)2, (3k)? = (k—k)?. [3.9] 
For the sake of simplicity, let 
7— 0 and k=0, 


which indeed can be achieved by means of a simple co- 
ordinate translation. Using [3.7] and [3.8], along with 
partial integration, we obtain: 


———; yA* A 

a= [arth \(- Sa) At) dh +) ee yaa [3.10] 
it ) yy* d 

a= [ora \(- a) 9) do— + [ == Sede. [3.11] 


We now come to a quantitative determination of the un- 
certainty relation. In order to calculate 2?-k?, we start 
from the inequality 


Dal ss ve)+ | > 
ah et wt 5 (Ss ee 


from which, by using [3.11], [3.5], and y(co)—>0, it fol- 
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lows that 
1 aa 
_— k2 . 
0 <| D(x) dx i 


Thus, we have k?-x?>} or, more generally, (dk)? (8#)?>4. 

With Shes all 
Ab= + VisKe, An=+V (82), Ap =1(8p)?, 

the uncertainty relation follows as a purely wave-kine- 

matical law: 


h 
Ak-Av >}, Ap: An>>. (3.12] 
The equality signs in [3.12] hold only when D=0: 
oy = aaa 
ants Dt 


The solution of this differential equation is a Gaussian 
distribution which, when normalized, is 


p(x) = = exp [— az*] ; ¢——— [3.13] 


Thus, a wave packet with a Gaussian distribution is asso- 
ciated with the smallest uncertainty. Spectral analysis of 
such a packet gives a frequency spectrum which is again 
a Gaussian distribution: 


— 
1 2 
Ath) = Se V2 [expt ax? — ikx|dax (from [3.2]) 
_ 4/24 exp|/—ala+ wy, 
SS | ee 2a 
-4_ fe few _k 
 /in ¥ an Va xP gal? 


af 
lee 


c. Behavior of a wave packet in time (t4~0) 


2 
dx:exp |- 5 


2 


[3.14] 


We start with the time-dependent form of the Fourier 
integral representation (in [3.1], kw—>(kx—at), w given 
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by [2.7]): 


i 2 
p(x, t)= vm fA exp f (i = mn) A 


uu . 
hk 
A(k, t)= A(k)-exp |- Unie 1 


= om |v t) exp[— tka] da 4 [3.16] 


| A(k, t)|?=|A(k)|?. [3.17] 

Now we choose A(k) to be the Gaussian distribution [3.14]: 
1 ke 2ihat 
p(x, t) = yea flr (147 )+ ikea. 
With 
= a 
*Oisre + 2ihat/m 

and 


ak 


dk = jexp - ~ (k — 2ixa)? 


X exp [— «x?] = 2/an exp[— ax?], 


Kes 
few - re + tke 


it follows that 


47 
2a, 1 - 
_—— SS SS Sere | —— 8.18 
ues) am V1-+ 2ihat/m lg) [ ] 


(For ¢=0, we again obtain [3.13].) Hence 


en eu. 
| p(a, t) |?= x VIL (Qaht/m)? exp [— («+ a*)x?]; 


2a 


oo Seman [3.19] 


From probability theory we know the variance of a Gauss- 
ian distribution: 


B 1 
Wiz)= wor=/2 exp([—fx?] implies Saas 
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Using this, we obtain 


sora nfs 


m 


Similarly, it follows that (Ak)*=a, with |A(k)|? given by 
[3.14]. Therefore, 


1 h2(Ak)? 
A) 


(62) Sn ae ae [3.20] 


The variance of a moving wave packet increases quadrat- 
ically with time. Not only does this hold for a Gaussian 
distribution, it is also generally true. With A(k,t), we 
now have 


* 
5 ae (k, t). ate t) ‘ar=[]>4O9 u Fax, [3.21] 


in analogy with [3.10]. From [3.16] we see that 


dA(K,t) _ hk? dA(k) Ak 
7) = exp| vone |: ( oi —74 ar a(n). [3.22] 
Therefore, 
dA* DA 


ia 


m= [| 2 


= am ime. ms 


dk Ok 
ql oA* 
AwrA_ recast 

= fi ( iF yak a 


4. WAVE PACKETS AND THE MECHANICS OF POINT PARTICLES. 
PROBABILITY DENSITY 


ht, [3.23] 


We introduce the momentum p in place of the wave 
number k: p=hk. 
A(k, t) 
’ [4.1] 


Vp, t) = =e 
HP) =e ACh) 


[p(p, t) |? dp =| A(k, t)|? dk. [4.2] 
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Equations [3.2], [3.1], and [3.4] are replaced by 


1 . 
Q(p) = van | 0) exp | jpelae, [4.3] 
~ i 
p(x) = sen | exp ki polap, (a4 | 
and 
fier ap=fiy|rar=3, [4.5] 
respectively. 


We first remark that the relation between wave packets 
and the mechanics of point particles can be only of a 
statistical nature. A measuring arrangement defines a 
State. For a state at time t, the wave functions y(«#) and 
y(p) are given. However, these are not physically mea- 
surable quantities; we can only measure the probability 
W(a)da of finding a particle between w and #+dz. We 
call W(x) the probability density; in terms of it we can 
formulate the basic assumption of wave mechanics: 


W(x) = | p(x) |?, 
the probability for the interval between x and x-+ dw is 


W (a) dax ; [4.6] 

W(p) = |9(p)|', 
the probability for the interval between p and p+ dp is 
W(p)dp. [4.7] 


With these functions we can write the expectation values 
as (compare with [3.5] and [3.6]): 


f(a) = Sf(a)>a = | f(a) W(x) de , [4.8] 


g(p) = <9(P)>a = [arr wap . [4.9] 
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The reason for introducing nonmeasurable quantities y 
and » is found in the fact that these quantities obey simple 
mathematical laws, particularly the linear superposition 
principle: If y,(#) and y,(x) denote possible states, then 
the state described by C,y,(7)+C.y.(#) is also possible. 
On the other hand, the probabilities—just as optical inten- 
sities—are not additive (cross terms). Rather, they show 
the well-known interference effects on the observation of 
which wave mechanics is actually based. 

Of especial note is the complete symmetry of the theory 
in W(a#) and W(p); that is, all formulas remain correct 
when the following substitutions are made: 


p x a 
| as [4.10] 


The completeness of the symmetry comes from the fact 
that the probability, as a real quantity, does not change 
aS a result of complex conjugation. 

“‘Any speculation which destroys this symmetry by intro- 
ducing intuitive pictures is not to be taken seriously.” 
This opinion was also shared by Einstein. He believed this 
statistical description to be indeed correct but not com- 
plete. Howeyer, up to the present day, possibilities for 
extending the theory have not been found, even though 
no proof of the impossibility of extending the theory has 
been given [A-1].4 


5. MEASURING ARRANGEMENTS. DISCUSSION OF EXAMPLES 


If we consider a current of particles which passes through 
two small holes in a diaphragm, then we discover that the 
probability of finding a particle behind the diaphragm is 
a typical diffraction pattern, as in optics (Fig. 5.1). This 


“Comments [A-1]-[A-5] appear in the Appendix on pp. 193-197. 
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interference of the probability is independent of intensity, 
that is, independent of the density of the particle current. 
It depends only on the locations of the holes; the holes 
define for us a ‘“‘state” with its associated wave function. 


vm 


Figure 5.1 


By measuring the position and momentum of particles, we 
can determine only statistical properties of this state which 
was fixed by the experimental arrangement. This statis- 
tical description prevents contradictions between wave me- 
chanics and the mechanics of point particles; however, it 
leads to an uncertainty that is characteristic of wave me- 
chanics. If we consider, for example, the ‘‘state’’ repre- 
sented by a single, force-free atom, we find that this state 
is completely changed by every measurement of position 
or momentum. If we measure the position of the particle, 
for example, then an indeterminable momentum will be 
transferred to the particle during the process; this will 
make a precise prediction of the position of the particle 
at a later time impossible. Although the trajectory of a 
celestial body can be determined more and more accurately 
by making successive measurements to which the laws of 
classical mechanics are applied (Fig. 5.2), every measure- 


Figure 5.2 
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ment on an elementary particle throws the particle out of 
its trajectory (Fig. 5.3); that is, earlier measurements of 
position are useless for a further determination of the tra- 
jectory. 


Figure 5.3 


Try as one will to think of experimental arrangements 
and measurements, there is always an uncertainty present 
which is given by the Heisenberg uncertainty relation [3.12]: 


h? 


(Ap)?: (Ax)? > t° 


(In Section 3 this was derived purely from wave kinematics.) 
Following Bohr, we call two quantities complementary if, 
like p and x, they satisfy an uncertainty relation. For 
example, EH and ¢ also are complementary. 

We now want to show by means of three examples how 
the uncertainty relation actually works. At the beginning 
we can remark that the conservation laws for energy and 
momentum are presently considered to be firmly established 
experimentally and theoretically even for a single process. 
Likewise, it is firmly established that intensity has no in- 
fluence on the interference phenomena. 


a. First example: Position measurement in a microscope 


From the theory of the microscope we know that con- 
vergent light is necessary for making a position measure- 
ment in a microscope (Fig. 5.4). The limit of accuracy of 
such a position measurement is given by Abbé’s sine con- 
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dition, 


In this way we could in principle measure «# arbitrarily 
accurately by simply making A small. However, we must 
require that the microscope function classically. That is, 
we must require that we have a large number of light 


~~ 
Ax 


Figure 5.4 


quanta, so that at least one quantum is certainly scat- 
tered from the object; we can then observe this quantum 
with macroscopic means (eye, photographic plate, etc.). 
As a result of the scattering, an indeterminate momentum 
is transferred to the object because, unfortunately, without 
disrupting the microscope, we cannot determine which path 
in the microscope the light quantum took. The region 
available to the light quantum is included in the angle e. 
From this, using 
_hy_h 
|p | a as a ’ 
it follows that the uncertainty in momentum is 


Apy~ sine 


and it is really true that 
Ap, Az~wh. 
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b. Second example: Momentum measurement using the Doppler 
effect 

We allow a wave train of light of finite length Z to im- 
pinge on the particle whose velocity v, is to be measured 
(Fig. 5.5). Again, the wave train will contain so many 


Ck) [AREA A aa 
zx 
> <_ 
v hy 


x 


Figure 5.5 


light quanta that certainly one of them will be scattered 
from the particle. We suppose that this light quantum 
approaches along the negative «x direction and is scattered 
along the positive # direction, and we set up the energy 
and momentum balance for this situation before and after 
the collision: 


momen- hy hy’ ' i hy hy’ 
tum St Ps=— tis: Pa Ps — 
energy hy+H=hi'+ EE’, h’/'=hvy+H-—E'. [5.2] 


Using the general relations (valid even relativistically) 


py) Ok ' 
dp SENS aia ass [5.3] 


and differentiating [5.2] with respect to p,, we find 


dy’ yp, 
h =), Pe 


YP, yp. oe 


Since » is fixed, differentiation of [5.1] with respect to p, 
results in 
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substituting this into [5.4] gives 


or 


a ee ee 5 
oe 7 ee 2— Vee [5.6] 
Written in terms of uncertainties with 


Av’ = (dv’/dp.) Ap. 
we find 


hAv'= oo pen (050, Apes. [5.7] 

The wave train passes the particle in a finite time 7; an 
uncertainty in »’,5 

Ay’ ate ; 

arises from this temporal limitation on the interaction. 


Thus, 


h 
Ap.~ Gear [5.8] 


An uncertainty in position comes about because, without 
ruining the arrangement for measuring momentum, we 
cannot determine when the particle changed its velocity 
during the time interval T: 


‘ 
Ag~ (v,—,)2. [5.9] 
Taken together with [5.8], this gives 
Auwv:Ap,~wh. 


which we had already obtained. 
5 This formula can be obtained, for example, by Fourier analysis of the 


wave packet under consideration. Moreover, it is identical to the uncer- 
tainty relation for the energy, AE: Aiwh, 
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c. Third example: Coherence property of light 


We consider a diffraction experiment of light falling on 
a diaphragm with two holes (Fig. 5.6). For interference 


Screen 


Diaphragm 
———_— 


Figure 5.6 


to occur, classical wave theory requires coherence of the 
light. What is the situation when we consider the process 
quantum mechanically and take a single luminous atom 
as the light source? 

The interference pattern on the screen is a system of 
bright and dark bands which depend on the path difference 


p=" (41,— 41). [5.10] 
If we want to determine through which hole in the dia- 
phragm an emitted light quantum went, then we can do it 
by, for example, measuring the recoil momentum of the 
atom. To that end, we must know the momentum of the 
atom in the w direction to a precision 


h. #d 
Ape< 2 qs , [5.11] 


As a result of the wave nature of the atom, this implies 
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an uncertainty in position 


hk il h A A 


i Sr 
7" 9 Ap, 2 2hsind/2  4sind/2° [5.12] 


As far as we know, the atom could just as well be at posi- 
tion A’. There the path difference is 


yy) 
p= (AL, — AL), [5.13] 
and we have 
Ce <a 
D—D'= 7 2 sing Ar>x ‘ [5.14] 


This means that the interference pattern is washed out as 
soon as we have determined through which hole the light 
quantum passed. As always, the uncertainty relation pre- 
vents a contradiction between the wave and corpuscle 
descriptions. 


6. CLASSICAL STATISTICS AND QUANTUM STATISTICS 
Classically, we can say that there exists a probability 
density W(p, 2) from which, by integrating, we obtain 
won) =|We, ade, W(x) =|we, a)dp, [6.1] 


where, however, W(p) and W(x) are different functions. 
The normalization is the usual one: 


[ [we eapar—1, or [Werae=1, Jroap=s. 


In this context, a measurement means a reduction of the 
probability, that is, a decomposition of W(p, x) into parts 
which represent subsystems for which the probabilities are 
additive: 


W(p, 2) = 9,W,(p, x) + g.W.(p, £) (0< g<1), 
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and also, of course, 


W(p) = Wi (p) + 92.W2l(p), 
W(x) = g.W,(2) + g2W.(2). 


Classically, we can continue this decomposition until x 
and p are known to lie in the intervals (#, +A) and 
(p, p+Ap), respectively, in which case we obtain the sim- 
plest distribution: 


W(a) =0 outside of (x, Aw + 2) 
W(p)=0 outside of (p,p+Ap). 


Here there is absolutely no connection between Ap and 
Aw (no uncertainty relation: h=0). However, we can de- 
rive a formula for the change of the probability in the 
course of time, for example, for a force-free mass point 
with statistically distributed initial conditions: 


a= a+ vt=a+ Pi (r= 2). 


With the simplifying assumptions 


3m=0, da,dp=0, 
we obtain from 


that 


This is exactly formula [3.20] which was derived on the 
basis of quantum theory. The only difference is that there 
is no connection between Ap and Az in classical mechanics. 

The specific differences between classical mechanics and 
quantum mechanics are 


1. the uncertainty relation, and 
2. the interference of probabilities. 
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It is these two points which lead to complications when a 
reduction of the probability is to be made. A state char- 
acterized by a py function cannot be decomposed into sub- 
systems for which the probability is additive. Whereas in 
classical mechanics a measurement signifies decomposition 
of the probability into parts which represent subsystems, 
in quantum mechanics each measurement gives a new 
state instead of a selection of substates. 

In quantum mechanics we can distinguish between two 
CASES: 

Pure case (Fig. 6.1): The probabilities are quadratic forms 
of py functions. This is true, for example, in the case of 
the diaphragm with two holes, both of which are definitely 
open. Then it is indeterminate through which of the holes 
in the diaphragm a particle went. 

Now we introduce a shutter which always covers one 
of the holes, but we assume that we do not know which 
of the holes is being covered at any given instant. In this 
case, we say that it is not known through which hole the 
particle went. We call this case a mixture. 

Mixture (Fig. 6.2): The probabilities are additive; the state 
cannot be described by a y function. Therefore, a decom- 
position into subsystems is possible until the state repre- 


Pure Case Mixture 
=D =p 
| ah Ye 


: 

br ety | 

ee ‘ 
=—— 


v(x) = x1 (x) + Cays(x) W(x) ant | vi(x)] 24 94| w(x) |? 
P(D) = Cpi(P) + Cr¥2(D) W(p) = ol i(p) |? + ga] os(p)|* 
lex|* + |ea|*=1 0<g<1 


Figure 6.1 Figure 6.2 
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sented by each subsystem is a pure case; that is, in the 
case of a mixture, the probability is a sum of probabilities 
for pure cases. A mixture is always obtained if one aver- 
ages over the phases, or if the relationship among the 
phases is destroyed in some other way. 

In principle, of course, it is also possible to distinguish 
between a pure case and a mixture in classical mechanics, 
except that the pure case is then trivial in the sense that 
p and # are known exactly. 


Chapter 2. Description of a Particle in a Box 
and in Free Space 


7. ONE PARTICLE IN A BOX. THE EQUATION OF CONTINUITY 


We first restrict our attention to a one dimensional case 
(Fig. 7.1). 


— >  - a 


Figure 7.1 


We find the corresponding wave function by solving the 
wave equation 


dp th wp 
ot Bm da?” ae 


with boundary conditions 
p(0, 1) = y(L, t)=0. [7.2] 


First, using standing waves as a starting point, 


p(x, t) = u(@) exp — 7B, b=, [7.3] 


we determine the solution for a stationary state. Then the 
probability 
W(x) = | p(a, t) |? =| u(x) |? 


23 
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is independent of time. The most general solution of 


h? d?u h? ” 
Ss ear or OT 
is 
ua) = A exp |i) + Bexp — iF. 


Satisfaction of the boundary conditions: 
y(0,t)=0: B=—A, or v= csinE 


y(L, Dia: Bae 


sas x 
=n, OF U,—= Cin TH), 


n= 0,1, 2,228 
Normalization: 


fingtae=|ol-de—s, \ol=|/5 


a x i i Ad 
ye, t) = yz sin (x Fn) exp|— + Eat ee Ee [7.4] 
Probability density in configuration space: 
t ae x 
W(a)= zs (xn) ‘ [7.5] 


Probability density in momentum space: The quantity 
\pnl| = nah/L is fixed; therefore, 


Wom= w(- i) = ee eles 


Because of this one cannot speak of the motion of a par- 
ticle in a stationary state. Motion is only possible if there 
is a linear superposition of stationary states (wave packet): 


2 x tha 
(a, t) = Vz 2 Cy, SiN (x i") exp |- mL? nt = et) 
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As can be easily verified, the orthogonality relation 


L 


fudtoyua(a) de = bn [7.8] 
0 
holds. 
In a three-dimensional region we have an analogous 
eigenvalue problem (Fig. 7.2): 


YS Uy, y=0 onS 
y= ze exp|— 7 Eat Un(x) : [7.9] 
Veu, + A,uwe=0, A= = 
8 
Figure 7.2 


At this point we want to show that the equation of con- 
tinuity, in the form known to us from electrodynamics, 
follows from the wave equation.! Here the probability 
current density appears in place of the electric current 
density ,? 

i= ail. (p* grad yp — yp grad y*) , [i100] 


2m 


1 Of course, here the equation of continuity is more generally valid. It 
includes the equation of continuity ef electrodynamics as a special case. 

2 Since we can only specify the probability of finding a particle in a given 
volume element, a “current” of this probability appears in place of a par- 


ticle current, 
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and the equation of continuity is 


m, yay + [tase 0 


on (7.20) 


Cire ce 
AC y) + divi 0 


If we multiply the equation [7.9] for wy by y* and the equa- 
tion for y* by y, the sum of the resulting equations is 


(yt V+ seal p* Viv — p Viy*) = 0 


The equation of continuity follows from this equation by 
way of the formula due to Green and Gauss, 


div (a gradb) = grada-gradb+ a Vb. 


In our boundary value problem, i=0 on S; therefore, it 
follows that 


| y*y dV = constant. 


Vv 


For the eigenfunctions we have the orthogonality relation 


| wiucdV=0, d,s, - [7.12] 


Vv 
Proof: If we multiply the equation [7.9] for u, by ux, and 
the equation for uw, by —u,, the sum of the resulting 
equations is 


(um, V2, — Uy V2um) + (An — An )utu,=0. 


Using Green’s formula again and applying the boundary 
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conditions, we obtain 


(A, — 42) futuav=o 


Vv 


n=™M: A, =A. OF Awrealy 


y ay Jotuav=o, 
V 


An=hn, %#m: This is called a degeneracy of the states 
nm and m. 


The degree of degeneracy is defined as the number of 
independent solutions belonging to the eigenvalue A,,. These 
solutions are not necessarily orthogonal. However, by an 
appropriate choice of basis for the eigenfunctions, they can 
be made orthogonal. 


8. NORMALIZATION IN THE CONTINUUM. THE DIRAC 6-FUNCTION 


If there are no walls or force fields present, then there 
is no longer a restriction on the wave function; that is, 
p varies continuously. We can represent y(a, t) by a Fou- 
rier integral (see [4.1] and [3.15], with p=hk): 


‘ : 
w(x, d= aq [ow t) exp a dp 
7 _ 
<5 fow exp p+ P| ap. (8.1) 


If we compare this to the linear superposition of eigen- 
states (see [7.7]), 


p(x, t) = > ¢,u,(@) exp \-; Ee, ; 


then it is a simple generalization to introduce a function 
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u(p, ©) in analogy to the u,(#), in which case we have 


4 p* 


p(x, t)= | p(p)u(p, ©) exp - On t) dp 


=|¢9(p, t)u(p, v)dp +. [8.2] 


al 2 
u(p, %) = V2 k pa| 


We write the orthogonaliiy relation formally as 


| ita Uen pal oan [8.3] 


This is a convenient symbolic notation; the symbol 6(p—p’) 
(Dirac 6-function) has the following meaning: 


Dy 
[ro —pap 
D1 
ao. if p’ is in the interval (p,, pz) 8.4] 
0, if p’ is not in the interval (p,, p,) . oo 
The “function” 6(@) does not exist in the proper sense [A-2], 
because it would have to have the properties 


_f9, for x40, 
ate) =| , for «=0. 


Rather, it is the symbol for a limiting process which is 
carried out in the following way. We consider a function 
6,() such that 


festeyan—a A 


Sec. 8 | NORMALIZATION IN THE CONTINUUM 29 


For example: 


ile 6,(2) = n exp[— 2n22?] 
0 
Figure 8.1 
sin nx 
2 6,(£) = 
ey 
0 
Figure 8.2 


@, for |v#|>1/2n, 


: onl n, for |x|<1/2n. 


Figure 8.3 


Then we have 


lim |f(p)6,(p — p’') dp 


P1 __ | f(p'), if p’ is in the interval (p,, pz), 
“10, if p’ isnot in the interval (p,, p,). 
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We can also say that 6(#) can be approximated by means 
of proper functions 6,(x); for example, we can allow the 
rectangle to become taller and narrower indefinitely, in 
such a way that its area remains constant. More generally, 


Oe) se | Auto) exp[ipa|]dp with Jontoiae— jie 
then 
Ann) = te) exp[—ipx]da, A,(0)=1, lim 4,(p)=1, 
thus 


+0 


uF ; 
On) ae exp[ipx|dp . 


—- oO 


This means that one symbolically writes 


(0) = | f(0)8(«) aa 


in place of the exact Fourier formula 


a 
0) = 5 fap [f= exptipalae- 
The generalization of the 6-function to three dimensions 
is simple: 
6 (x — x’) = 8 (a, — 7) +8 (a0. — 4) +5 (23 — 5) . 


Integration of [8.3] over p yields 


De 
fafeinararwa'.9 
Py 


- 1, for p’ in the interval (p,, p.), 


0, for p' not in the interval (p,, p.). Sh 


This expression has the same meaning as [8.3] and [8.4]. 
Now, we can also demonstrate the correctness of [8.3] 
because the u(p, #) given in [8.2] satisfies [8.5] and, there- 


See. 8 | NORMALIZATION IN THE CONTINUUM 31 


fore, the orthogonality relation identically. That is, with 
Dp 


if. i 
ah exp|—> px 


aa 


dp 


"Vee | Ae 


it follows from [8.2] and [8.5] that 


a 
— eXD |— 7 Paw 


) 


Ee | a, XP [(t/h)(p'— p,)#] — exp[(t/h)(p'— p,)z] 
27 1x 
__ 1 [sim {(1/h)(p'— ps)2} — sin{(1/h)(v'— pada} 
27 2 
_f1, for p’ in the interval (p,, p,), 
~ |0, for p’ not in the interval (p,, p.) . 


This is correct because 


a Ce for a@>0., 


a x —1, for a<0, 
A | = 

a cos ba 4g. 

~ Lv 


9. THE COMPLETENESS RELATION. EXPANSION THEOREM 


Let f(#) be a function such that Jif(a) pdx exists. If 
we have a complete orthonormal set of functions {u,(2)}, 


that is, 
frstoynate) da = dam » [9.1] 
then we can make a series expansion 
f(x) = J Antin(2) , [9.2] 
in which 
A, = [fot dx. [9.3] 
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We require convergence in the mean for ihe series in [9.2]. 
This is a weakened convergence requirement because pos- 
sible ‘‘spikes’’ are integrated over: 
N 
lim i Ifa) — Santa(a) [2dr = 0, [9.4] 
el 
or 


lim [|2(2)|*2— 0. [9.5] 
Because of [9.1] and [9.3] we have 
foxitey— 3 anuatoy|*= facie) Sa, [meyud(e) de 
— Sar heme den +3 ata, = / ifltde— aka, [9.6] 


Because the left side of [9.6] is certainly not negative, 
Bessel’s inequality follows: 


N o 
> aan <]lfj?dz, or > |a,|P] |fltde. [9.7] 
n=1 n=1 


Bessel’s inequality implies that [9.4] means the same as 


Slett= inter. (9.8) 


This equation is known as the completeness relation be- 
cause it guarantees us that none of the u, were left out, 
that is, that the set of uw, is complete.® 

The generalization to two or more functions is simple 
with the substitution 


f(x) —> ef(%) + e.9(x), 
An —> 0,0, + O2dn, b, = foceyutto dex . 


* Since we are dealing here with sets of infinitely many functions, we 
cannot establish completeness by counting, for example. 
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Substitution into [9.8] and comparison of the coefficients 
of cic, yields, for example, 


> dnb, = i f*(x) g(a) da . [9.9] 


By means of the 6-function, the completeness relation 
can also be written symbolically as 
> uh(x) up(@') = 0(@ — 2’). [9.10] 


This can be verified in the following way. If we multiply 
the above equation by f*(x’) and integrate, then we obtain 
with [9.3] 

> vfta) nto (x')da! = f*(a) . [9.11] 
Multiplying this equation by g(x) and integrating once 
again yields 


> a(e)ud(o) def ua(a')f*(a") da! =[*e)a(a) ae, [9.12] 


which is [9.9]. 

Remark: Equation [9.11] is in fact correct in practically 
all cases which arise in physics. However, since we have 
only required convergence in the mean for the series, we 
must integrate again for the sake of exactness. That is, 
we must go over to Eq. [9.12]. 


10. INITIAL-VALUE PROBLEM AND THE FUNDAMENTAL SOLUTION 


We call that solution of the wave equation which de- 
scribes a particle that was definitely at position w’ at time 
c—0, 

y(z,0)=d(4— 2’), [10.1] 
the fundamental solution K(a, 2’, t). 
We first consider a particle in a box. The general solution 
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for that case is (see Sec. 7) 
w(x, t) = 2 on n(@) exp +5 i. [10.2] 


By using the completeness relation [9.10] we can also write 
down the fundamental solution immediately : 


y(t) = K(a, a! t) = ¥ walaryun a’ yexp|— 7 Bat], [10.3] 
or, for a three-dimensional region with arbitrary boundaries, 
K(x, x,t) = 2 tale) ur(x’) exp |— 


; t [10.4] 


Whereas only a linear superposition of stationary states is 
possible as the wave function of a particle in a box (the 
boundary conditions must certainly be satisfied at all times), 
there is no such restriction for a particle in free space. 
In that case we have the continuous variable p instead of 
the discrete index n (see Sec. 8), and the wave function is 


a ; ao 
y(@, t) = Joon, %) EXP re J dp 
ee Beha 10.5 
am Pi; \P?— a p [10.5] 
with the normalization 
[urw, x)u(p', x) dx = d(p— p’). [10.6] 


The completeness relation [9.10], which is now written in 
the form 


i u(p, ©)u*(p, 2!) dp = 6(@— a") [10.7] 


(note the symmetry between w and p!), again allows us 
to write down the fundamental solution immediately: 


w(a, t) = Jute our, @')exp|— 7 Fe tla, [10.8] 
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or, by comparison with [10.5], 


y(2, t) = K(x — 2", 1) 


_ : eae 
mom? ip plo—a)— F-| dp. [10.9] 


We can perform the integration: 


z } a 1m (w—«a’)2 
Kla—a',t) = exp |e OEP oe] 


4 (x—ax')m\? ¢t 
x ye —— Sa | 5 
Jexp h (p t Vea ep. 
The substitution 
(2 —a')m 2mh 
es | ae 


leads to a Fresnel integral: 


«at 
Fe 


exp [—iu?]du = /%-exp |- it 


Jexe [iu?]du = 1/%-exp 


afi 
=; (ieee 


The solution K(#—z2’', t) can now be written as 


4 m(a— 2x’)? 
f=’, t) = —_—— /moex ah sae | [10.11 
Te 2t 
Using [10.10] it also follows immediately that 
[xe- © ,tjag—1. [10.12] 


Now, without using the Dirac function, we want to show 
that the K(«—za’', t) defined by [10.11] has the properties 
of the fundamental solution. Let =a—z2'. We shall cal- 


culate 
és 
lim i K(é, 1) dé 
é 
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With the substitution 


we obtain 


Se 


fac dé— =a My ie f exp[iu®]du. 
é. 


1 


There are three cases to be considered: 


a. Se 0, &,>0: “7 0 U, > 0 
t>0: w4-+0e 4&->+co : 
ot ; lim | Kdé=0 
by €)=0, 6, = 0: u,<0 Uz< 0 t—>0 
t+0: uw2--cC Uu—>—co ‘ 
&, 
@. €:— 0, 6,20: U,< 0 Us, > 90 im | rae 
t>0: U>—CO %—>+0o} t>0 


The formulas in [10.10] were used in evaluating the last 
case. The three cases a, b, and c taken together are equi- 
valent to the equation 


K(é, 0) = d(€), 


and it remains to be shown that K(&, t) satisfies the wave 
equation 


4 


MK Oth BK 


‘Ut 2m de * 


OK Vm Jaen pie | iL @ Fa f m2} 


at WV Oahi\ 2ivt v4 h 20 h ot 
0K _y/ m1 im limes 
ME VV Onnt yy te ™P Roe |? 


h 2t 


eK ,/m 1 fim pee a m& _ 2m ok 
2G? Onhi gli we P 
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It is worth noting that the fundamental solution [10.11] 
is well known from the theory of heat conduction. If we 
replace the mass m in the nonrelativistic wave equation by 
a purely imaginary quantity, then indeed we obtain 
the heat conduction equation, and the Gaussian distri- 
bution [10.11] is the familiar solution for a 6-function heat 
pulse at time ¢=0. 


Chapter 3. Particle in a Field of Force 


11. THE HAMILTONIAN OPERATOR 


By means of the Hamiltonian function 


3 | 
Hay Viens fi 


baa 2M 
the equations of motion of classical mechanics can be 


written in the canonical form 


dp, oH OV 


dt soba, ae 
dt |p, m 


We obtain the corresponding equations for a particle of 
charge ¢ in a magnetic field H by introducing the vector 
potential A (H=Vx A): 


el é 2 
H= = on) {Ps ae x(a) “OV (eee) [11.39 
dp,, 0H OV € e 0A, 
dé dk, ey ta Pi, 5 | 


dm, H_1(, ey — 
dt 3p eee Pr . 


Again, we use the operator formalism given in [2.4] as 
a translational key between classical mechanics and wave 
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mechanics: 


ed) d 
— ith — E->th—. é 
Dp—>—i ne cs [11.5] 
As a result of introducing the operators, the Hamiltonian 
function becomes the so-called Hamiltonian operator: 


Pa) 
H(p:, %:) >H (=. », [11.6] 


Generally, an operator A is understood to be a calcula- 
tional prescription which associates a function uw with an- 
other function, called Avw.? 

In the force-free case (V=0), we can also write the wave 
equation [2.11] with the operators [11.5] and [11.6]: 


Hy= Ey. [11.7] 


This equation also gives the correct result for a particle 
in a field of force (V40). Written out, the equation then 
reads 
h? .. oy 
—_V? = ih—; 11.8 
it was first suggested in this form by Schrodinger.? 

It so happens that Eq. [11.7] is also correct in more 
complicated cases, for example, in the case of a charged 
particle in a magnetic field. Nevertheless, an ambiguity 
sometimes arises when the classical Hamiltonian function 
is replaced by the Hamiltonian operator. Whereas the or- 
dering of factors in a product is of no importance for clas- 
sical quantities, 

Di®Xe = LeDi » 


with operators the ordering is very important. For ex- 
11t is not necessary that the prescription be a differentiation. For ex- 


ample, the operator x is defined to mean multiplication by =. 
2. SCHRODINGER, Ann. Physik 81, 109 (1926). 
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ample, later we shall find that 


Thus, it is possible for several different operators to be 
derived from a given Hamiltonian function. Which of 
these different possibilities is correct can be decided only 
by experiment; that is, the expectation values derived 
from the solution of [11.7] must follow the classical tra- 
jectories. 


12. HERMITIAN OPERATORS 


The solutions of the wave equation [11.7] must satisfy 
the important condition that their normalization be con- 
stant: 


> fey a=). [12.1] 


(This condition can easily be understood physically by 
remembering the interpretation of y*y as a probability 
density.) ; 

Carrying out the differentiation in [12.1] and using the 
equations : 


h dp | bh dy* 
nny ff =e = * 
i Oe = and i (Hy) ’ [12.2] 
we obtain 
i [y*(Hy) — y(Hy)*]d%x = 0. [12.3] 


This is a condition on H. Any operator which satisfies this 
relation is called Hermitian. 
Since the superposition principle should be valid for the 


* Cross terms, such a8 p;%,, which lead to ambiguities occur with non- 
Cartesian coordinates, for example, 
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waves, we must require 


H(ey) = eHyp | 
[12.4] 
A(y, =f Ws) = Hy, ae Hy, { 
An operator with this property is called linear. 
If we substitute 
YHYit Ye 
in Eq. [12.3], we obtain 
Jivtatys — p,(Hy,)*]d*x = 0 [12.5] 


for a linear Hermitian operator H. If F and G are linear 
Hermitian operators and a and b are real numbers, then 


we have 
aF + bG& Hermitian , [12.6] 
FG+GF Hermitian, eee 
i(FG—G@F) Hermitian, [12.8] 


but F@ in general is not Hermitian. 


Proof of [12.7]: From [12.5] we have 
| {(H,y:)*Heyp. — y( ApH y,)*} da = 0, [a] 
[overtop — vi Hepy ae = 0, 
i (Hy tH, — tH Hy sds =0. — [b) 
[a]—[d): 
fivtttly. - vB a=. [el 


Substitute H, =F, H,=G in [ce]; then substitute iG, 


42 PARTICLE IN A FIELD OF FORCE | Chap. 3 


H, =F and add: 


| (yt FO + GF)y, —ylG@F+ F@)y)}de=0. QED. 


Examples of Hermitian operators: The operator 
p =—th(d/dx) 


is Hermitian, as can be seen by partial integration: 


h dp h (ryt h >p\* 
ES a a, =—_- ——_ = — — ; 
Jv er al aie {i 7 ee 


On the other hand, ip is not Hermitian. The operator 
p? =—h*(d?/dx") is also Hermitian because 


iG ze — — v) dz 0 (Green’s formula) . 
Any operator which prescribes multiplication by a real 
function of x is clearly Hermitian. Thus, the potential 
V(x) is Hermitian, and this means that our Hamiltonian 
operator, 
2 


h2 


is also Hermitian. 

H remains Hermitian even in the case of a magnetic 
field. When we write the Hamiltonian operator corre- 
sponding to [11.3] we need only pay attention to the or- 
dering of the factors: instead of writing 


(2 7 e€ e 
we must write 
é td e  @ e? 
(p.— Ay) = Di — ¢ De An — GAeBe + a Aj. 


In order to prove the Hermiticity of the Hamiltonian oper- 
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ator that is derived in this manner from [11.3], all we 
need show is that 


PuAn + Ax pe 


is Hermitian (special case of formula [12.7]). This follows 
by partial integration: 


13. EXPECTATION VALUES AND THE CLASSICAL EQUATION OF 
MOTION. COMMUTATION RELATIONS (COMMUTATORS) 


So far we have learned to calculate the probability dis- 
tribution for observable quantities by using wave functions. 
If classical mechanics is to be contained in wave mechanics 
as a limiting case, then the expectation values of these 
quantities must obey the classical equations of motion. 
We now demonstrate that this is the case. 

First we calculate the time derivative of <x,> (see (ideD |, 
as well as [12.2] and [12.5]): 


C2,» = | p*a,yp de ; 


(0) =* (Ha, — 2,2). [13.1] 
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In a similar way, we can derive 


d i , 
di = ;, He — p,H>. [13.2] 


Now, we must evaluate the right-hand sides of [13.1] and 
[13.2]. In general, an expression of the form 
ER be Lb (es! [13.3] 
is called a commutator. The relations 
[F, Pay Ps) = Falls, fs) 12s, F] F, [13.4] 
and 
[¢,F, + oF, Fs|=e¢,{ Fi, Fs|+ of f., Fs] [13.5] 
follow immediately from the definition [13.3]. 
In agreement with the formulas for Fourier transforma- 


tion, [4.3] and [4.4], we now define the operators analogous 
to [12.9] associated with the momenta and coordinates: 


h 9d 
Pyy(*) = ore L,y(*) = Tey , [13.6] 
he 
P_P(P) = Dr@ 5 x,.9(p) = — pen ae [13.7] 
Px 


4 


With these definitions we can immediately derive the basic 
commutation relations: 


PxPi1— PiPe= 9 : [13.8] 


For example: 
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Using [13.4] and [13.8] we obtain 


a (Mis — Xr A = Dd {Pi Pike <_< XL, P;) aie (Pi X, cay XL, Pi)Pi} = 


Here we have introduced the new operation of differen- 
tiating a function of operators with respect to an operator. 
We define that this differentiation is to be performed in 
exactly the usual way in which a function is differentiated 
with respect to an ordinary variable. 

By repeated application of this procedure we obtain 


hoor 


= 13.9 
ond [13.9] 


F(p)*, — & F(p) = 
for an operator # which is given as a polynomial in p,, 
or, in the limiting case, as a convergent power series in Pi 
(that is, as an analytic function of p,). Of course, the 
relation 


F(p)P.— Px F(p) = 0 [13.10] 


also holds. In an analogous way, using [4.3] and [4.4], 
we obtain 


G(x) P, — Py U(x) = —=—, [13.11] 


G(x)a, — 0,G(x) = 0. [13.12] 


Summation of the last four equations yields 


h dH 
— = —_— ot 13.13 
H p,— Py i 3a,’ [13.13] 
h oH 
— —? ——— 13.14 
He, 0, H 4 Pe [ ] 


for an operator 


H(p, x) = F(p) + Ge). [13.15] 
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Thus, the result applies to the Hamiltonian operator with 
a potential. It is easy to show that [13.13] and [13.14] 
are also valid for an operator 


H(p, x) = P(p) + Gx) + 3 (Ante x)Px4-PrAr(x)} [13.16] 


(Hamiltonian operator with a magnetic field). If we sub- 
stitute [13.13] and [13.14] into [13.1] and [13.2], we obtain 
the classical canonical equations of motion for the average 
values of x, and p,: 


@ oe, = 7 OB d ony = 
at C#y> =a (22) i) dt Py = -G » {13.0 


(That is not purely by chance; rather, the theory was so 
constructed that [13.17] would result. It will always be 
required that a reasonable theory contain classical me- 
chanics as a limiting case.) 

As a generalization we proceed to consider the case in 
which an operator F depends on time explicitly (for ex- 
ample, time-dependent forces). From 


(FP) = | p*(Ey) dea 


’ 


and the wave equations [12.2] follows 


a oF 
GE = 5 [ey By) — peta + [ye Ey aze 
; a: 
=|o" areata, (see (12.5) 
i ; 
Kee -G (H, FI +=). [13.18] 


For F=H we obtain 


d oH 
qe = (2 [13.19] 
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If )H/dst=0, it follows that <H>=constant; this is a 
statement of the law of conservation of energy. 

As an example, let us specialize to the case of a mag- 
netic field: 


ii e 
a. (vt a (DP, Ay + Ax Px) + = At) + ee: [13.20] 


Using [13.17] we obtain 
d 1 e 

dt (%,)> = = ( = <A.) [13.21] 
d , 1 ¢ dA; vA; 
= (on SOE Get He. 


A; 41; 
35 (= ye. Ae * =) 
A TT Oa; OD, 


0A, e ov 
—|p,—-A;|\t-—-— ). 13.22 
Ye (2 a= I =) once! 
Note the symmetrization, which is necessary for preserving 


Hermiticity. 
By not using the canonical form we can also write the 
equation of motion more simply: 


<mx,> = (bs = - Ay )), ite — i, [13.23] 


where the definition 


is used and 


av ; ‘ 
fo ee (R24 ey), [13.241 
OL, 2c t 
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with a magnetic field (written as an antisymmetric tensor) 


_ ee [13.25] 
aa OL, Ov; 


and with an induced electric field 

£,==-22 2. [13.26] 
Equation [13.22] can be easily derived from [13.23] and 
[13.24]. 


We also have a continuity equation if there is a mag- 
netic field (see [7.11]): 


< (yty) + divi=0, [13.27] 


where the density i of the probability current is given by 


; 1 e € 
y= oe ty" (p.— = ds] v—v(2e+ <4.) vt [13.28] 


or 


. h oy ap* € 
= * = a * 
. ami? Dery vie) per y [13.29] 


Gauge invariance 


As we know from electrodynamics,‘ a gradient can be 


added to the vector potential A without changing the 
magnetic field: 


t 
A,~> Ay+ ate?) [13.30] 
k 


However, because of [13.26], an extra term must then be 
added to the potential, 


“See W. Pau, Lectures in Physics: Electrodynamics (M.I.T, Press, Cam- 
bridge, Mass., 1972). 
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this means that the form of the Hamiltonian operator is 
changed. However, from the wave equation [12.2], it fol- 
lows that the only change in y is that it is multiplied by 
a phase factor: 


yopexp 


ae 


The group of transformations defined by [13.30], [13.31], 
and [13.32] is known as the gauge group. Correspondingly, 
those quantities which do not change under these trans- 
formations are said to be gauge invariant. Examples 
of such quantities are the probability density y*y and 
the probability current i. 

As a further example we consider the virial theorem.5 
The quantum theoretical expression 


[{m d? © )) 


2 ae | 
=> S (K 2, “fb ok ,) + > (@, a, sie iid) [13.33] 


follows from the relations 


& (er) = ©, + Wy [13.34] 
and 

a2 

aa (2,4,) = Bay, + 2%, + 00, + Uap; [13.35 


for i=k, the expression reduces to 


m d? : 
G ap (7 » = (Kya, + map» . [13.36] 

Formulas [13.33] and [13.36] are also valid classically ; 
but the symmetrized form is required by quantum mechanics 
for K, (Eq. [13.24}). 


‘See W. Pavuti, Lectures in Physics: Thermodynamics and Statistical 
Mechanics (M.1.T. Press, Cambridge, Mass., 1972). 


50 PARTICLE IN A FIELD OF FORCE | Chap. 3 


From 
0”, = ©; — <4; [13.37] 
it follows that 
<(82,)*> = <ai> — (<a>)? [13.38] 
and 
(8a, 3a,)> = (00, — (UD e) ’ [13.39] 


as we know. Using these relations, we can derive the 
more general expressions 


d?2 
(2o a (82, bo) = (S (5 (Ei, Ban+ Bod) 
+m (86,86, 4 84, bis) ) [13.40] 
and 


m a? 


oF an! (32;) y")= CK, 8x, + m(8z;)?> . [13.41] 


In the force-free case (A, = V= 0, ma@, = p,), we again have 
the formula analogous to [3.20]: 


m? d2 
e oF aa (8n4)* ) = <(Sp,)2>. [13.42] 


Chapter 4. More than One Particle 


14. MORE THAN ONE PARTICLE 


When there is more than one particle, a new feature of 
the problem is the interaction between particles. According 
to the theory of relativity, one should take into account 
the finite speed of propagation ¢ of the action of these for- 
ces. In relativistic quantum theory the finite speed of pro- 
pagation will be taken into account. In our nonrelativistic 
approximation we set ¢= co for the sake of simplicity; that 
is, the action of a force propagates instantaneously to all 
particles. For this reason, we are able to use a single time 
coordinate t, which means that time is not considered on 
the same footing as the 3N space coordinates. We shall 
assume the number of particles, NV, to be constant; thus, 
no radiative dissociation processes will be considered.} 

From here on we designate the space coordinates by 


Qiy very Or (f= 3): 
a, es 5", OPP, ..., 03 } hy, GOs, Fy Is, -.-9107} - 
Similarly, we designate the momentum coordinates by 


{Diy +++» Prt 
and define 


d/g = dq, 'dq.°...° dg; , d’p = dp, dp,°...-dp;. 


1 See W. Pauwt, Lectures in Physics: Selected Topics in Field Quantization 
(M.I.T. Press, Cambridge, Mass., 1972). 
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Further, we introduce two generalized wave functions 
(diy +--+» UY, t) and (pi, ..., py, t) which are defined such 
that 
W(di, --95 Gs, dG = p* pd’ , [14.1] 
W (Dry +++) Ds t)d’p = p*pd’p . [14.2] 


These are the probabilities for finding the coordinates 
and momenta of the nth particle between g, and q;+dq,, 
and p, and p,-+dp,, respectively (¢=3n—2, 3n—1, 3n; 
i=, ...7ee)e 

In analogy to [4.3] and [4.4], we have the following rela- 
tions: 


al! 
OC neey Usy )= Geaym [ole Pr t) 


xX exp +5 (P19: ide... pai) d’p, [14.3] 


1 
(Pry -++9 Do, t) = aan [ves wey Vy, t) 


i 
X XP |— 5 (Pada + Pade + +++ + par| d’g. [14.4] 
The wave function y is to satisfy the wave equation, 
4 h ay 
—= =p, [14.5] 
where again H is a linear, Hermitian operator: 
[vityaa= p(Hy,)* aq. [14.6] 
Generalized to the 3N coordinates, the relations 
h 
Ape An >>, [14.7] 
h 9 h 9d 
aah dr een and wade 7 [14.8] 


are also valid in this case. 
Now we consider the special case in which there is no 
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interaction between the particles, which means that the 
Hamiltonian operator reduces to a sum of independent 
terms: 

H=Ho+ HO 41 ...1 H™, [14.9] 


Here H® is only to operate on the ith particle. If 


P(N» Tay Ia)y -oey yo '(In-2s Gn-1) Iw) [14.10] 


are solutions of the wave equations 


=Hyo (4=1,2,...,N) [14.11] 


of the isolated systems, then 


yp = py. mae [14.12] 
is a particular solution of 


—~— = Ay=(HY+ Ho+...+ Hy. [14.13] 


The general solution is a linear combination of products 
[14.12]. 

Thus, decomposition of the Hamiltonian operator into 
independent terms according to [14.9] corresponds to a 
reduction of the wave function to a product [14.12] of 
independent factors, in agreement with the fact that the 
probability reduces to a product when the particles are 
statistically independent: 


W(Q15 ++ Wry t) = Wi (x, ft)... Wx(x™, t) 
W (Pr, +++) Pry t) = Wi(p™, t)... W,(p™, t) 
Therefore, we know the Hamiltonian operator for un- 


coupled particles in external fields of force (see [11.3] 
and [13.20]): 


[14.14] 


N 
Hy => H®; [14.15] 


a=1 


(a) h? a 3 (a)/ (a) : (a) / (a) 14.16 
H = Se a er: nae) Ve) [14. J 
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If the forces between the particles can be derived from 
a potential V(q), for example, in the case of a Coulomb 
interaction, 

Cal 


Vi sec qd) = > 


ad Tab 


: [14.17] 


where 7,,= |x — x” |, then it is appropriate to write 
A= H+ V(Qi5 ---> MY) : [14.18] 


The commutation relations [13.8], as well as formulas 
[13.13] and [13.14], are also valid in the case of more than 
one particle; the only modification necessary is that the 
indices ¢ and k run from 1 through f. 


Chapter 5. Eigenvalue Problems. 
Functions of Mathematical Physics 


We have already solved a particular eigenvalue problem 
in Section 7, the problem of a particle in a box. Now 
we have the necessary tools available for treating the more 
general eigenvalue problems of particles in more compli- 
cated fields of force. 

Since we seek stationary solutions (standing waves) in 
eigenvalue problems, we use 


p(x, t) = u(x) exp 5m 
as a trial solution, in which case the wave equation [11.7] or 
[11.8] leads to the time-independent Schrédinger equation ' 
Hu= Eu. [T] 


We now solve this equation in some simple cases. 


15. THE LINEAR HARMONIC OSCILLATOR. HERMITE POLYNOMIALS 


For the linear harmonic oscillator we have the well- 
known relations 


ee eae eee 
a ome 9 ods =k, [15.1] 
aH : 
a ae: 
- ames >Gtopq=0. [15.2] 
— pm 


1 ¥rom a historical standpoint, it should be noted that the time-independent 
Schrédinger equation was found before the time-dependent one: EK, SCHRODIN- 
GER, Ann. Physik 79, 361 (1926). 


55 


56 EIGENVALUE PROBLEMS | Chap. 5 
From [15.1], corresponding to [I], we obtain the time- 


independent Schrédinger equation for the one-dimensional 
harmonic oscillator: 


h2 ad2u 


2 
eer: de ats — Wo92u = Eu. [15.3] 


With the coordinate transformation 


em |re-a : a a ee = < y ee 
and with 
A= i. ; | [15.5] 
[15.1] and [15.3] simplify to 
H =}3(p2+ 2°) hop, [15.6] 
—o8 youn, [15.7] 


and the commutation relations [13.8] become simply 


f il 
Pa 5" [15.8] 
With the substitution 
a? 
u= yexp|— a [15.9] 


[15.7] becomes the differential equation 
y"— 2ay'+ (A—1ljyy=0, (15.10] 


which is solved by the Hermite polynomials. The Hermite 
polynomial of degree n is given by 


H,(%) = (— 1)" -exp[ 


a). [15.11] 
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By defining 


n 


d 
4 = Z,0 OxPl— x], [15.12] 


we obtain with the polynomial rule of differentiation 


a Ane jan 22 exp[— x?]) 
eat 
= Aaa a exp[— 27], {15.13] 
oy — Qa = 2(n-+1)yx [15.14] 
yo + 2aey’+ 2(n+1)y=0. [15.15] 
From 
y=exp[@’]y, xz=exp[—z* ly, 
x = exp[— 27] (y’— 2ay) , 
z'= exp[—2°}{y"— 4ay'+ (40? —2)y} [15.16] 
follows 


y"— 2xy'’+ 2Qny=0, [15.17] 
and this equation is identical with [15.10] if we choose 
A=2n+1. [15.18] 


With this result, the solution of Eq. [15.7] can be written as 
x? 
h,(@) = ¢, exp |- a GE) [15.19] 


with c, a constant and the eigenvalues (see [15.5]) 

E,, = ha,(n+ $) . [15.20] 
Even in the case n = 0 there is a nonzero energy HE, = ha,/2 
(zero-point energy). Now, we show that the eigenfunc- 


tions h,(x) are orthogonal; that is, they satisfy [7.12]. 
According to [15.7] and [15.18], we have 


h, + (2n + 1—2*)k, = 0, 
hy + (2m+1—2°)h, =0. 
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Multiply the first equation by h, and the second by —h,. 
By adding the resulting equations and integrating, and 
using ' 

hnhn oa hike (limln a POE 


we obtain 
+0 


2(m —n) | hahndt = (mh, —hyh,)|=0, [15.21] 


which for nm, shows the validity of the orthogonality 
relations. 

Before we can normalize the solution, we must derive a 
few more properties of the Hermite polynomials. The Her- 
mite polynomials [15.11] are alternately even and odd: 


H,'x) = 1, H,(x) = 2a 

H, (x2) = 4a? — 2, H,(#) = 8a° — 122 [15.22] 
HA,(a) = 16a4 — 48a? 4-12, 

Their generating function is 


i)? a 


f(a, t) = 5 <P n= exp[ot] $ 


n=0 n! 


= exp [x?] exp|— (w —t)?]= exp i t24 2tx]. [15.23] 


pe eePl—- 2") 


From this, 


= 2tf [15.24] 


follows immediately, and, by comparing coefficients, we 
obtain 
Hy (x) = 2nH,-4(a) . [15.25] 
From 
d 
a + 2(t— a)f=0 [15.26] 
we obtain 


wH,(%) = $° Hy+1(%) + n° Hy4(2) [15.27] 
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analogously. Of course, these relations, [15.25] and [15.27], 
also could have been derived directly from the definition 
[15.11]. From [15.25] follows 


ar 


age Hal@) = 2n-2(n —1)+...+2+1-Ho(w) = 27! [15.28] 


for the nth derivative. Now we can determine the con- 
stant c, in [15.19] by normalizing 


fae da = et fexp—a"]E(2)a2— ae [15.29] 
With [15.11] and [15.28] we find 


Gt = (—1)"|{ Hala) exp[— 2*]do 


-/(3 c exp[— a ]daw=2'n!/m. [15.30] 


The differential operator was moved to the left by partially 
integrating n times. Finally, with c, computed from [15.29] 
and [15.30], the solution [15.19] reads 


1 ie ; 
oe 9 | ae): [15.31] 


a phase factor for this solution can still be chosen arbi- 
trarily. From [15.27] and [15.31], we immediately obtain 


Dh, (#) = bert + An—1(£) \" ; [15.32] 


and, from [15.25], [15.31], and [15.32], we obtain 


hal = V2nhn- 1(@) — hy, (x) 
a ou hn+i(®) a : + h,- (2 aE : [15.33] 
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The completeness of the Hermite polynomials 


Now we want to show that the set of functions h,(x) is 
complete; that is, we want to show that Eq. [15.10], 


y"— Qay'+ (A—1ljyy=0, [15.34] 


has no regular solutions other than the h,(#) (that is, no 
regular solutions except for A=2n-+1 and n integral). To 
this end we look for the general solution of [15.34] and 
show that it is regular only for A=2n+1. We take a 
power series, 


=> ae, [15.35] 


as a trial solution. Since there is no singularity at 7=—0, 
there are no negative exponents in [15.35]. 

At this point we want to formulate a general theorem 
which is useful for the wave equation Hy=Ey: If H is 
an even function of the p’s and 2’s, as is normally the 
case when the interaction potential is even (V(x) = V(—2)), 
then every solution of the wave equation is the sum of an 
even solution and an odd solution (as is easily seen with 
the aid of the wave equation). 

In particular, our Hamiltonian function [15.6] is even, 
and, since u—e*"*y, our theorem also applies to the so- 
lution of [15.34]; that is, [15.35] can be written as the 
sum of two power series solutions with 


P00 2,4): [15.36] 
and 
? = Taos [15.37] 


By substituting [15.35] into [15.34] and comparing the 
coefficients of 2°, we obtain the two-term recursion for- 
mula 


a,*p(p —1)+ a,-.{— 2(p — 2) + A—1}=0. [15.38] 
Let po be the smallest exponent; then, according to [15.38], 


Sec. 15 | THE LINEAR HARMONIC OSCILLATOR 61 


we have 


Po(Po—1) = 9: P=0 or po=1. 

even solution odd solution 
This shows that the initial values in [15.36] and [15.37] 
are correct. Now we want to investigate if and when the 
power series [15.35] terminates. The value p,,, is to be 

the largest exponent occurring in the power series: 


a,—0 for Di tne 1 5 @,7-0. [15.39] 
Then, with p=—n-+-2 in [15.38], 

0=—2n+A-1 [15.40] 

follows; that is, the power series only terminates for 
— 2 te [15.41] 
For a given odd value of A, either the even solution ter- 
minates or the odd solution terminates, depending on 
whether » is even or odd. Thus, there is exactly one poly- 
nomial; for any other A neither solution terminates. Now, 


we want to clarify the convergence properties of the non- 
terminating series. From [15.38], 


y= Oya er : [15.42] 
follows 
tg ey IO PEE (15.43) 
for p even, and 
a= a, BATA) (2p —3— A) wis, da 


p! 


for » odd. For sufficiently large p,, we certainly have 


pe) 
Du > ; [15.45] 


thus, the coefficients [15.43] and [15.44] no longer change 
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sign for p>p, and we may choose a4,>0. Further, 


lim p-—2 =2 [15.46] 


po Ay-2 
follows from [15.42]. As a simplification we introduce 
p= 29 sot 7 — 29-1 (¢@=0,1,...); [1547] 
then the series [15.35] becomes 
Y= Daye OF Y= HY Aa B*. [15.48] 
q 


q 
We shall not be concerned further with the factor 7 and 
with the difference between a,,;, and a,, in the right-hand 
formula. From [15.46], 


lim ie 


we have 


Aog — Aog—2 . fo) 


>6 or Ang S |g [15.49] 


1 
for given 6<1 and sufficiently large g. By repeated ap- 
plication of these expressions, we obtain 


Cc: 62 
; Ag > a : [15.50] 


using this in [15.48] we can write 


O° £\@ 
y=) onsen: SY ( rr =c-exp[+da?] [15.51] 
q @ = 
and 
x a 
u = exp |— a > c-exp (0 — 5) a [15.52] 


Since we can choose 6 arbitrarily close to 1 and thus cer- 
tainly greater than 4, the normalization integral of wu does 
not exist; that is, the nonterminating power series do not 
provide normalizable solutions of [15.7]. The h,(a) are the 
only regular solutions of Eq. [15.7], and this concludes the 
proof of their completeness. 
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16. MATRIX CALCULUS ILLUSTRATED WITH THE LINEAR HARMONIC 
OSCILLATOR 


If we have a complete, orthonormal set of functions {u,}, 
futuae =o, [16.1] 


then we define the matrix element of the operator / with 
respect to the uw, as 


(k|F|n)= Jui Fu,da. [16.2] 


The quantity (k|F|n) is the element in the kth row and 
the nth column of the matrix of F; this element is also 
designated as F,,,. In this way we have associated the 
operator Ff with a matrix which we also designate by F. 
This notation will be justified in the following by the 
equivalence of matrices and operators. 

By means of [16.2] the condition [12.5] for Hermiticity 
of an operator F' can be written as 


(| F|n) = (n|F[k)*. [16.3] 


By analogy, we call a matrix satisfying this condition Her- 
mitian. 

Because of [16.1] and [16.2] we can make the following 
expansion: 


PFu,(#) = > u(e)(k|E\n) . [16.4] 
k 
For example, for the operator 2 we have 


Lunt) = Yu(a(k|a|n), (kla|n)=(n|a|k)*; [16.5] 


and for the operator p, we have 


du, il 
dz 


oat Dui(w)(k|pelr),  (k| Palm) = (n| Pel k)*. [16.6] 
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The multiplication of operators corresponds to the mul- 
tiplication of the associated matrices; an example is 

du, 

Pa{run(ax)} — » —+4 7 


k 


= Smla) SC) pal) (el e|m) = Yule) (U|pew|n). [16.7] 


This agrees exactly with the usual definition of a matrix 
product: 
(1| AB|n) = > (| A|k)(k| Bln). [16.8] 


k 


Of course, the relation 
x(Paun(x)) = >) ur(x) (U| x p2|n) [16.9] 
t 


also holds, which means that the commutation relation 
[15.8] takes the form 
Pot — Pe= (—1)1 [16.10] 


in the matrix representation (/ is the unit matrix, (k|Z|l)=6,,). 
By comparing [15.32] with the expansion [16.5], we see 
that only two of the matrix elements of the operator x are 


nonzero if we expand in terms of the eigenfunctions 
h,(): 


(n+1|e|n) = = [16.11] 
(n—1|"|n) = Ve [16.12] 


Because of the Hermiticity of #, [16.12] follows from [16.11] 
and vice versa: 


n-+1 


(n|x|m+1) = (n+1|a|n)*= ; 


[16.13] 


The relation [16.12] follows from this by means of the sub- 
stitution n->n—1. We obtain the matrix element of p, 
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from [15.33]: 


P2h(L) aaa (= 1) 


dh,(z) i +1 


dz 9 has: (x) 


= yt Pane) — ,» h,(a)(k | pe|n); [16.14] 


(n+1]p,|) = iW nel [16.15] 
(n—1|p,|n) = aye. [16.16] 


Written out, the matrices are 


thus, 


0 v1 oO 0 
‘ VI 0 v2 0 
e=—=-| 0 V2 0 V3 
0 0 V3 0 
~ a eran ee ee . 116.17] 
VJ1 0 0 77 
—vV/1 0 V2 0 
aes 6 =—-v2 0 v3 
= iV 2 
0 0 —v3 0 | 
Note that the indices begin with zero. With 
i 0 v2 0 0 
0 3 0 V6 0 
—_ V2 0 5 0 v12 
mo 210 /6 0 7 0 
0 0 vVi2 0O 9 
1 (ear One ile ae 
0 3 0 —V/6 0 
i 1 =—~/72, 0 : 5 0 =a le 
a. | = 6 0 “i 0 
0 0 —/12 0 9 
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we obtain 


(hoy) 1H = 3 ("+ p2) =" [16.19] 


o Oo S&S ft 
oe 8 ww © 


for the Hamiltonian matrix. As a result of expansion with 
respect to the eigenfunctions, the Hamiltonian matrix is 
diagonal, and the diagonal elements are the energy eigen- 
values. This is true not only in the above example, but, 
rather, it holds more generally, as we shall see later. 

For some applications it is also useful to introduce the 
following matrices: 


0: 44aO 0) 
1 d 0 V2 0. 
A= > = 75 aah 
| eae, re 
0 0 0 0°! 
: . qi a V1 0 0 0 
0 0 V3 0 


Then the relation 
(k| A*|n) = (n| A|k)* 
holds, and the relation 
[A, A*]=AA*— A*A= 1 
follows from the commutation relations [16.10]. 


a. Equivalence of the Schrodinger equation to a system of equa- 
tions in Hilbert space 


According to [9.2], we can expand an arbitrary quadrat- 
ically integrable function y(x) in terms of a complete 
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orthonormal set of functions u,(a 


= > uae) , mie a)dv. [16.20] 


In this way, a vector with components y, in an infinite- 
dimensional Hilbert space is associated with each func- 
tion (x). In this case the completeness relation is 


il p(x) |?da = 2 |v [2, [16.21] 


from which follows the existence of > |y,|?. Therefore, 
k 


instead of calculating with wave functions, one can just 
as well calculate with vectors in Hilbert space, applying 
the rules of matrix multiplication. 

If we substitute the expansion [16.20] into the Schrédin- 
ger equation, multiply by U;,(@), and integrate, we then 
obtain 


feeutte Dy (2) “pa f de¥(e D> Lune) ya. [16022] 
Combined with [16.2] and [16.4], this yields 


> (k|H|n) — £,(k|1|n)}y,=0. [16.23] 


This linear, homogeneous system of infinitely many equa- 
tions (k=0, 1, 2,...) in the unknowns y, is exactly equiv- 
alent to the Schrédinger equation and gives the same results. 


b. Example: Linear oscillator with an additional potential 


With the following considerations we want to prepare 
the perturbation theory. If an additional potential V(x) 
is added to the Hamiltonian function H, of the linear har- 
monic oscillator (we put ba,=1 for the rest of this sec- 
tion), 


H=}(a? + pz) + V(a) = Ho+ Via), [16.24] 


68 : EIGENVALUE PROBLEMS Chap. 5 


then the corresponding Schrédinger equation, 


Hy=E£y, 
or c 
(H—#-1lhy=0, 


in matrix form, is usually no longer exactly solvable. If the 
additional potential is small (a ‘‘perturbation’’), then an 
approximate solution can be found which is especially easy 
to calculate in the matrix representation. 

Since the h,(x) are complete and orthonormal, we can 
make the expansions 


y(a) = 2 hrla)n » po=[ytorn a)da, [16.25] 


a)h,(@) = > h(x) (k| Vn) 


(k|V|n)=[ hee) V(ayh,(o) da]? ‘88-281 
and there is the relation 
fisey(a) ae = (k|1|n).- [16.27] 


As we did with [16.22], we substitute this expansion into 
the Schrédinger equation, multiply by h(x), and integrate: 


foonzte ) XK Ao — E)h,(&%) + V(a)h,(@)}yn=0. 
We then obtain 
Di(n+3—E\(kl1|n)+(k|V[n)}yp,=0. [16.28] 


This is again a system of infinitely many equations in the 
unknowns y, which is equivalent to the Schrédinger equa- 
tion. We shall solve the system later when we consider 
perturbation theory. 
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c. Determination of the eigenvalues of the linear harmonic 
oscillator with the matrix method 


We have just derived the matrix representation of quan- 
tum mechanics from the wave equation, and we have shown 
that the two types of representation are completely equiv- 
alent. Nevertheless, historically, matrix mechanics was de- 
veloped before wave mechanics,? and at first the identity 
of the two theories was not at all noticed. 

As an example of a matrix mechanical calculation without 
recourse to the wave equation, we treat the linear har- 
monic oscillator again. Thus, in the Hamiltonian matrix 
for the linear harmonic oscillator, 


H=}(pi+#), [16.29] 


where, for example, 2? stands for the matrix 


(n|at|n’) = ¥(n| a/R) (el ln’), 


we must substitute Hermitian matrices for # and P2 Which 
make H diagonal. In addition, 2 and p, must obey the 
commutation relations [13.13] and [13.14], which, of course, 
can be derived within the matrix mechanical framework. 
In our special coordinates [15.4], the commutation rela- 
tions are 


H 
Ha — 2H =~ id= = — De. [16.31] 


2 W. HEISENBERG (Z. Physik 33, 879 (1925)) introduced the matrix elements 
as the quantum mechanical analogues of the Fourier amplitudes of classical 
mechanics. Just as a classical quantity is determined by its Fourier ampli- 
tudes, so the associated quantum mechanical quantity is to be given by the 
set of associated matrix elements (although with the matrix elements a form 
corresponding to the Fourier integral is not possible). However, Heisenberg 
did not yet use the expression ‘‘matrix element” in his first article. Born and 
Jordan (Z. Physik 34, 858 (1925)) were the first to recognize that the law 
of multiplication given by Heisenberg for the quantum mechanical quantities 
was identical to the rule for matrix multiplication. By means of the matrix 
calculus the entire theory could then be still better founded and worked out 
(M. Born, W. HEISENBERG and P, JoRDAN, Z. Physik 35, 557 (1926)). 
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If we now assume that H—E is already diagonal, as a 
result of the correct choice of 2 and p,, then such a com- 
mutator becomes simply 


(n'|HF—FH|\n") = (Ly — E,+)(n'| BF \n") 
(Pauli: ‘‘... just calculate it!”), which, with [16.30] and 
[16.31], yields 
(Ey — Ey*)(n'|pe|n") = i(n'|2|n") , 
(EZ, — Hy) (n'|a|n") = —i(n'|p,|n") . 
Therefore, (n'|p,|n") and (n‘|z|n") are either both zero or 


they are both different from zero; in the latter case, as 
a result of 


(Ey — Ey»)? (n'|@|n") = (n'|a|n"), 
we have 
(By — By)?=1, or Ey —#Hyv=+1. [16.32] 
Since the eigenvalues of the linear harmonic oscillator are 


nondegenerate, as can also be easily shown within the 
framework of matrix mechanics, the selection rule * 


¢ n'—n=+1 [16.33] 


follows from [16.32]; simultaneously we have fixed the 
order in which the eigenvalues are to be numbered: 


HE, =n + constant . [16.34] 
With this numbering, only the matrix elements 
(n|2|n+1), (n|aln—1), (n|pz|n+1), (n|pz|n —1) 


are different from zero; all others equal zero. From the 


’ Only those matrix elements whose indices obey the ‘“‘selection rule’”’ are 
different from zero. We shall see later that the intensity of light emission 
is proportional to the square of the matrix element associated with the tran- 
sition responsible for the emission. Thus, the selection rule specifies between 
which states transitions are possible. 
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above relation, 
(n|p,|n F1) = +i(n|2|n F1) [16.35] 


now follows. 

We have not yet been completely justified in writing 
Eq. [16.34]. We must first show that the sequence formed 
by the £, has no holes; that is, we must show that the 
next larger energy following each H, is E,4,=E,+1. For 
this purpose we consider the diagonal element of the com- 
mutator [16.10], 


(n | pax — xp»\n) = 2{(n|e\n+1)(n+1]2|n) 
—(n|a\n—1)(n—1|2|n)}=1; [16.36] 


in arriving at this expression for the commutator we have 
used [16.35]. As a result of the Hermiticity of 2, we obtain 


|(mj2|n+1)|?—|(n|2|n —1)|*=3 [16.37] 


from [16.36]. That is, we have 


|(n|aln—1)|2= 5+ constant , [16.38] 
since this implies 
1 
[(m[@|m+1)|?=[(n+1|2|n) |? = a + constant 


which then leads to [16.37]. Using formula [16.35] and the 
Hermiticity of 2 and p., we also obtain 


(n|p3|m) = (n| a2 |) 
=|(n|a|n—1)|*+|(n]aln-+1)|*, [16.39] 


with which we then obtain 
E,,=|(n|a|n—1)|?+|(n|a]n+1)|?=n-+ constant [16.40] 


from [16.29]. The left side of this equation is again posi- 
tive. Therefore, there must be a smallest n=, such that 
E,=0 for n<m. Since we have only decided about the 
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sequence of the numbering until now, we can take m =0. 
Then we have H_,=0, which means 


(—1|#|0) =0; 
together with [16.37], this yields 
E,=%3 and H,=n+4. [16.41] 


The E, is the famous zero-point energy which we now have 
derived correctly with the matrix method also (see [15.20]). 


17. THE HARMONIC OSCILLATOR IN A PLANE. DEGENERACY 


The isotropic harmonic oscillator in a plane has the po- 
tential 


V(Gi, 43) = 5 Oagi + q3), [17.1] 


which leads to the Schrédinger equation 
d’u  d@u 
da ag ee [17.2] 
with 
A=2E/ho, and 2,=Vma,/h-q; 
' (see [15.4] and [15.5]). 
Equation [17.2] is the sum of two equations in a, and a,. 
Therefore, the solution reduces to the product 
U = hy (Ly) "Rn (eo) , [17.3] 


and the eigenvalues are 


A=A,+A,=2n,+1+ 2n,+1 
= 2(m + M+1) = 2(n+1). [17.4] 
(When the solution is a product of this form, we say that 
the equation is separable.) Here we have the situation in 


which several different states are associated with the same 
eigenvalue; this is known as degeneracy. If n states belong 
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to the same eigenvalue, we speak of n-fold degeneracy. 
In our case the degeneracy is (n-++1)-fold: 

wt: 0, 1h Dyce tl 

No: n,n—1,n—2,...,0. 
We can remove this degeneracy if we appropriately modify 


the potential—for example, by going to the anisotropic 
harmonic oscillator: 
si 


V(qi, We) = 2 (39; + ©3493) - [17.5] 


The equation is still separable and the solution is 


“= ries as}-he i) ) [17.6] 


E= h{a,(n,+ $) + @,(n,+ })}. arg | 


with 


If w,/w, is irrational (w, and w, incommensurable), then 
the degeneracy is removed. 

By going to the limit w,/w,—1, we again obtain the iso- 
tropic harmonic oscillator. Quite generally, a degenerate 
system can be thought of as obtained, usually as a limiting 
case, from various separable nondegenerate systems; that 
is, the degenerate system will be separable in various coor- 
dinate systems. 


a. The solution of the plane harmonic oscillator in polar 
coordinates 


For example, the isotropic harmonic oscillator is also 
separable in polar coordinates: 


L,=7COsp, 2,=PaMyp , dz,dz,=rdrdp , [17.8] 
Da es 1 of ou 1 0?u 
(Sat ae) a ae) ap 


P) 
Ou at ou 1 du 


~ ore or | 7? dg?- [17.9] 
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The Schrédinger equation is then 


ou 1 du. il du : 
ee = =. elton 
ee ee [17.10] 


The solution can be separated as a product of the form 

= 0.(7) ee, (17.11] 
where m is a positive or negative integer; then the differ- 
ential equation for +,,(r) is 


Gea, | IL Gk ape 
ee a 4 —r)v, = 0. IL Ue 
a! pa a 


Since the energy values of a system are not changed by 
a coordinate transformation, the degeneracy must be the 
same as in [17.4]; on the other hand, the eigenfunctions 
will be different. Indeed, even in one coordinate system 
they can be made quite different by taking linear com- 
binations. We can also remove the degeneracy by going 
to the isotropic anharmonic oscillator: 


(hw,)-2V (r) = $r?+ e V(r). [17.13% 


The e is to be a small number which is not proportional 
to r?. The separability [17.11] is not altered and we 
obtain 

d?v,,  1dtv, mM? 


+——~—- — 


dr? ip (Ole r? 


Um + (A—r? — 2eV(r)) on, = 0 [17.14] 


in analogy to [17.12]. The eigenvalues of this equation 
are still two-fold degenerate for the values +m; in order 
to remove this degeneracy we would have to introduce a 
magnetic field. 

Now, we want to solve Eq. [17.12]. Since the physically 
allowed range of r is 0<r<oo, we can simplify [17.12] by 
the coordinate transformation 


ey 


an P=A/e r— C—; 
ue or a’ 


[17.15] 
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we then obtain 


1od/. du Py) ou ou dU 
(+ =) > ( » 1(ore 


A—« 
4 


GO = 0 - 
With dy/dz=y' and the substitutions 


—— golim|/2 e-alty : 


Me gala pinta} ay! 1 mii 
smemnsnlyeg (2, 


| ee ye a ted esl as 2 
- 27° 2 Qa?" or 


we can write Eq. [17.16] as 


7, (m+1 Pee ir) - 
y's ("EP ays 5 (F—m—1)y=0 
or 


wy”+ (m+1—ax)y'+ ky=0, 


where 


As a result of the relation [17.4], we also have 


n=2kim. 


75 


[17.16] 


fe lees | 


[17.18] 


[17.19] 


We treat the case of m and k integral and nonnegative. 

The treatment of the general differential equation [17.19] 

yields physically useful solutions only for k=0, 1, 2,.... 
In the following we want to discuss some functions which 


solve [17.19]. 
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b. The Laguerre polynomials 
The Laguerre polynomials are defined by 
Iya) = £5 (ore) = 3 (-1"(") 
xX k(k —i)-...:(n+1)-a", [17.20] 
and they satisfy the differential equation [17.19] with m =0, 
aly + (1—a#)Iy+kL,=90. [17244 
This can be shown, for example, with the generating 
function for the Laguerre polynomials, 


© Ix(0) ,  exp[—a/(1—2)] 
EN 2. k! = ae ’ 


which we derive generally in [17.35]. Differentiation with 
respect to # and z and comparison of coefficients yields 


Li ah. [17.22] 
Lo (2k+1— xv) Ly, — k? Ly) . [17.23] 


The differential equation [17.21] follows immediately from 
this (substitute Z,-, from [17.23] into [17.22], and elimi- 
nate L,;, with the relation L,,,—(k+1)L,=—(k+1)L,). 
With the definition 
we obtain 
a( Le)’ (m + 1 — @)(Le)'+ (k — m) Le = 0 [17.25] 


from [17.21] by differentiating m times. If we replace k 
in this equation by k+m, then we obtain exactly the dif- 
ferential equation [17.19], which proves that 


ve aa ee [17.26] 


satisfies Eq. [17.19] for m>0. 
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Therefore, the solution of Eq. [17.16] (and the corre- 
sponding Eq. [17.12] with which we began) is 


Vrm = Constant x w™/? €-2/2 D(a). [17.27] 


We shall normalize this solution in [17.38] and generalize 
it to the hypergeometric functions in [17.37]. 

The Ly, are a set of orthogonal functions; that is, they 
satisfy 


[tn lP ena — 0 for kk. [17.28] 


of which anyone can convince himself by performing the 
calculation. 
By applying the formula 


f(x) = Bel ae dt, [17.29] 


where k>0 and C is a circle about 2, to f(x) =.x*e*, and 
using [17.20], we obtain 


—ttk 
Ly (a) = ek! dae aia [17.30] 


— i) ae 
By substituting —t+ <2 for t, this can be written as 


! 
L,(2) = aoe $ ett — wrt Dat . [17.31] 


Cc 


Replacing k by k+-m and differentiating m times with re- 
spect to #, we also obtain an analogous expression for L”,,,: 


_ da” k {]2 
Rem = Fm Latm = (— ies —— 


a a a)*t-%+mtDdt, [17.32] 
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Using this practical integral representation, we want to 
derive the generating function for the Laguerre polynomials: 


~ Ex) ee 1 ty * dt 
ja 7 t= spe S| ; :) a [17283] 


We can assume |z|<1; in this way the convergence of the 
series is guaranteed for sufficiently large ¢: 


_— 2 
| a ‘lela. 


The calculation yields 
M 
iGae a ve dt ] det dt 
2nt (1 =F t -Qaet t—(t—2x)z 


= ‘spe ———- (iae 


The integrand has a simple pole at the point t = — (xz) /(1—z). 
Therefore, the integral equals 2m times the residue, or 


271 ex as 
> P T= sP 


é 


and the final result is 


{(@,2)= ) 2 ; [17.35] 


The generating function for L”,,, follows from this simply 


by differentiating m times with respect to a: 


x Lie(x) _ q\m 9m OXP[— w2/(1 — 2)] 
i (1—z)m4 


Replacing k by k+m then results in 


& TEiml) oy (_yyq SXPL— a2/(1—2)] 


foo (+ m)! (1— 2)mH [17.36] 
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¢. Normalization of the solution [17.27] 


Now, we want to normalize the solution of Eq. [17.16] 
which is given by [17.27], 


Vx,m(@) = constant x 7”? e-2/2 D(x), [17.37] 


and prove its orthogonality at the same time. We calculate 
the normalization integral, 


fonmts’ndar ene OL. ae), (e)\daw—wW,,,0q, [ns] 
G 0 

by the elegant method of E. Schrédinger; this method can 
be applied completely generally (as can the derivation just 


given of the generating function). To this end we use 
[17.36] in the form 


[17.39] 


(ne nn &xP[— at/(l1—1)] 
Pe we —e : 


and we integrate the product of [17.36], [17.39], and 2” e~* 
with respect to 2: 


@ © gkgpkl re ra o 
2 2 as eeen if x2” exp [— 2) Lis m(@) Ly +-m(X) dx 


oo 


7 ; per forex {201 an 012) }ae 


(1— 2)" 


ill Zl m! 
= aes [ee -alty = Gp 


0 


aes ( "| 2. (kK+m)! 


k=aQ 


k oe [17.40] 
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For this calculation we used the following relations: 


pa [ureray=mi, 
0 


a at)” 


_ (m+1)(m + 2)-...-(m +) (m+ kt [17.41] 
k! kim! 


By comparing coefficients in [17.40], we find that for- 
mula [17.38] is satisfied with 
[(% + m)t]§ 


Nim = mer ak [17.42] 


Indeed, the coefficients of all terms of the form 2*t* with 
kk’ on the right-hand side are zero; the orthogonality 
follows from that fact. 

Remark: We define 
(*) a(a—1)-....(e—n +1) 


= = [17.43] 


also for nonintegral values of a. Formula [17.41] is valid 
in this sense: 


( )ear = ee Ia +n) 


n n! ~— T(a)ent * [17.4g 


d. Some properties of the I" function 


Here we gather together some properties of the J’ func- 
tion which we shall need later. The /" function satisfies 
the well-known functional equations 


T(e+1)=2-Ie), [17.45] 


T(z): T(—2) [17.46] 
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liz=—n (n=0,1, 2,...), then the first equation with (1)=1 
reduces to 


I'in+1)=n!. PL7.A47| 

Iz) has simple poles at z=0, —1, —2,..., and it is reg- 

ular everywhere else; 1/J(z) is an analytic transcendental 

function. The J’ function has the Euler integral represen- 
tation 


Pe+1)=ferra for Re(z)>—1. [17.48] 
0 


In addition, there is the simple, useful, and beautiful Hankel 


relation, 


1 1 
—— =— fe't-* liz A 
We oni e* t-*dt for all z, [17.49] 
a 


where the path of integration C, as shown in Fig. 17.1, 
goes around the origin and converges toward the real axis 


C 
i 


Figure 17.1 


at —oo. Hither one of the integral representations can be 
transformed into the other by means of [17.46]. Because 
of [17.45], the relation 


EO es : [17.50] 


yey pl) (y+ n— 1) 
also holds. 
e. The confluent hypergeometric function 


The general differential equation (see [17.19], m+1—y, 
k—> —«a) 
xy" + (y—ax)y’—ay= 90, [17.51] 
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where a, y, and # can be arbitrary real or complex 
quantities, leads to the confluent hypergeometric function 
F(a, y, ©), which is a limiting case of the generai hyper- 
geometric function F(a, B, y, a): * 
: x 
F(a, y, %) = aa fF (>. B, 3) : 
First, we want to derive an integral representation of the 


hypergeometric function and then investigate its asymp- 
totic behavior. To solve [17.51] we use the series 


yf =) aa [17.52] 
n=0 
as a trial solution. By comparing the coefficients of 2" we 
obtain the recursion formula 


(n+ 1)Nan41+ y(n 4+ Lani — na, — aa, = 0, 
atn 1 


yin Pee . [17.53] 


On+1 = 
In order that the solution not vanish identically, we must 
assume 4,40. From [17.53] it is also evident that no so- 
lution of the form [17.52] exists for y<0 and integral. 
If we set a,—1, then the solution is 


a(a-+ 1) 
2!y(y +1) 
a(«+1)...(#+n—1) 
niy(y+1)...(y+n—1) 


Pla, y, a)=1+ Fo + Cree ee 


a+... [17.54] 


‘The general hypergeometric function is defined by the hypergeometric 
series 


ab 1 
F(a, b,c, 2) = 1+ — poe Dea 
lle 2!ce(e+ 1) 


in which all quantities can be complex. The series is absolutely convergent 
for |z|<1 and divergent for |z|>1; it is also absolutely convergent for 


|z| = 1 if Re(a+b—c) <0. The hypergeometric function satisfies the equation 


a’u du 
#(1—2) 5 + {e— (a+ b+ no es 
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For «=0, —1, —2,..., —k, ..., the series terminates, which 
means that F(—k, y, x) is a polynomial; as we shall imme- 
diately see, for y=m-+1 these are essentially the Laguerre 
polynomials. 

Using [17.44] and [17.49] we can transform this solution 
into an integral: 


Fla, y, 0) =I(y) 5 i. "| ee 


n=0 n y Ae n) 
al <. {—« — x\" 
— t —¥ | ——_ 
ronda 3( "eee 
Cc 
The requirement 
“ Et, [17.55] 


which must be made in order that the series converge, 
means that the path of integration C (see formula [17 .49]) 
must surround the point x as well as the origin (Fig. 17.2). 


Figure 17.2 


Using the binomial theorem, 
o /_¢ — g\" n\-% 
oat | ea 
DC ea ( A 


F(x, y, @) = ion | et t-"(t — w)-2 dt. [17.56] 


we obtain 


Qt 
0 
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Substituting 
y= mm +1 ; (+ oe k 


into this formula (in which case [17.51] reduces to 
Eq. [17.19]), we again obtain the Laguerre polynomials 
of [17.32]; the only difference is that they were normal- 
ized differently: 


_ (e+ m)i(- 1)" 1 


Dyna (2) k ! Qn1 


| et t-e+m41)(t — dt 


(4 


=n" ™\a+ m)!F(—k,m+1, 2). [17.57] 


We now want to check that the solution [17.56] which 
we derived actually satisfies the differential equation [17.51] 
for the confluent hypergeometric function. To this end 
we state the identity 


oF" + (y— 0) F'—aF 


=O ay Fier —aytyat, (17.58) 


oO 


which can be derived immediately with the aid of [17.56] 
and the relations 


et e-7(t— a)-* (fc a ot (y—a) a} 


= et t#-7(t — 4)-7-1(— a) {(-#-1) + (a+1ap)+th 


We know that the path of integration C (see Fig. 17.1) 
runs from —oo, around the origin, and back to —co. How- 
ever, at the end points —oco the integral vanishes, so that 
the right-hand side of the identity equals zero. This shows 
that the equation is satisfied. 


t 
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If the path of integration C is broken up as shown in 
_ Fig. 17.3, then we have 


(a 
+ + 4 
, FP=F,+F,, 
with 
1 = ee fe 1°-7(¢ — x)-*dt. (17.59) 


1 

Cs 
According to [17.58], since the end points of the integration 
are all at —oo, both F, and F, are solutions of the differ- 
ential equation [17.51]. Thus, we have split our solution F, 


Figure 17.3 


given by [17.56], into the two independent solutions F, 
and F,. If —«=k=0,1, 2,..., then there is at most one 


pole at the origin and 
Tie =O ; F => He 5 


for y=m+1>—k this again yields the Laguerre polyno- 
mials. On the other hand, if the origin is a regular point, 
which means «—y integral and nonnegative, then we have 


Po Or, ar ae 
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f. Asymptotic behavior of the confluent hypergeometric function 


With the aid of the integrals just introduced we can 
investigate the behavior of F(a, y, 2) for large |a]. 

First, we calculate F,. For large |x| we can use the 
expansion 


= (—2)-* (1+ att... [17.60] 
From [17.49] 
ill eee = iI 
a fe ilt= [17.61] 
and with [17.45], 
i iL 
Fy—e—1) Pyoay * 
we get 
_— Ty), fy elev ED) 
P= Ta! a) {1 a top. [17.62] 


However, at one point in this derivation we ‘‘swindled,”’ 
because the’ condition |t/z|<1 for the convergence of the 
expansion [17.60] is only satisfied over one part of the 
path of integration (see Fig. 17.3). Therefore, we have 
to deal here with a so-called asymptotic series.» We must 
not take too many terms in such a series; otherwise the 
approximation will become poor again. For y—« negative 
and integral, the result F,=0 also follows from the 
asymptotic formula [17.62]. 

The expansion is completely analogous for F,. We need 
only make the substitution 


t—“#=T 


* For more information on this subject see, for example, E. T. WHITTAKER 
and G. N. WATSON, 4 Course of Modern Analysis (Cambridge University Press, 
New York, 1962). 


| 
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F(a, y, 2) = = ofa +t)*-77-*dt. [17.63] 


2 


This integral is formally the same as the one in [17.59] if 

we substitute («—y) for —a and —w# for -+-# in the in- 

tegrand. With these substitutions we can immediately write 
down the solution from [17.62]: 


Pyle, 7,2) = peony {y 4 LOM) [17.64] 


Of course, this too is only an asymptotic series. We again 
find F,=—0 for the case of the Laguerre polynomials, 


-o—k=0,1,2,.... 
With the asymptotic formula 
Ty) = 
Bia, y; eg em ye) 
e ne eae 4..., [17.65] 
and with 
y=m+l, —a=k, [17.66] 


we can now write the solution [17.27] for the plane har- 
monic oscillator, 


Vim = aml? e~2!2 B(— k, m+ 1, @). [17.67] 
For large a, 
m! 
rc ptink gomatk enue 
eo Pmt 
1" _g-(mi2+h+1) etai2 , 17.68 
Ta ] 


In the solution [17.67], # is real (see [17.15]); from this it 
follows that the second term in [17.68] diverges for large «. 
It is for this reason that the only allowed eigenvalues for 
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a physically possible (that is, orthogonal and normalizable) 
solution are k=0, 1, 2,.... The asymptotic solution of the 
plane harmonic oscillator is then 


m! 
Dip jen ~ etink pn ak e-ze : 17.69 
Gs ini +k) 
18. THE HYDROGEN ATOM 


a. Separation of the wave equation in spherical coordinates 


We now consider the Schrédinger equation for a central 
force field with potential V(r), 


h2 
a, Sa 4p Rs 
5 V+ Vir) u=EH-u, [18.1] 


written in spherical coordinates 7, %, m, which are related 
to the Cartesian coordinates 2, y, 2 (Fig. 18.1) by the 


Figure 18.1 


well-known formulas 
x =rsind cosp 


y=rsind sing ¢. [18.2] 
z2=rPrcos? 


Then the Laplacian operator takes the form 


ar a Cx) ee | Do! ... Sw 1 d@u 
py ee =) a hears ae 
aS T 7 sind od (sing =) : sin?? a3 [18:3} 
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or 
Yum For +S. 
By assuming the form 
u = v(r)- XY (4, ¢) [18.4] 


for u, Eq. [18.1] splits into a radial part and an angle- 
dependent part: 


il Gl dv 2mr2 
a =e VATE 
(r =) + re {EH — V(r)} 


ee aa lege 
= ae 35 (sin 09 5) +a f. [18.5] 


sin? 3 sin?) dp 


Since the left side of this equation depends only on r and 
the right side depends only on # and 9, the two sides can 
be equated to a constant A and separated: 


1d ( 2dv 2m A 

ava (r a+ fe [Z— V(r)]— — v=0, [18.6] 
ee IU AN 

sind 30 ( = T Fined sin?d dy? eo ; Gl 


The only physically useful solutions of Eq. [18.7] are for 
A=ll+1), C0 ih, 2 yee 5 [18.8] 


they are the spherical harmonics Y,(0, py). Here, we con- 
tent ourselves with writing down the most important prop- 
erties of these functions, since we shall be concerned with 
them in more detail in an exercise (see Sec. 43). 


b. The spherical harmonics 


With condition [18.8] we now write the angle-dependent 
differential equation [18.7] as 


Oo) oY E28 
ee sing —— ——+4114+1)Y=0. [18.9 
sin? od (' ‘ rt aad dy? ae?) [18.9] 
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This equation follows from the requirement of a central 
force field; however, it is independent of the particular 
form of the potential V(r). We solve it with the trial 
solution 


Y(9, p) = 0(d) emo , [18.10] 


in which the requirement of being single-valued leads to 
integral values of m. In this way we obtain the differential 
equation 


jeed if 
a(t ind 55) + re ae 


m2 
ey \0= 0.  [ienay 


nd 


With the substitutions 


ee: sind dd = — dz 
= » lasek 
8(8) = y(a) » Whence in SE = — (1—at)y! ali ] 


the equation is transformed to 


(1— a? )y"— 2ay'+ {U-+1)— 7 by = 0. [18.13] 


The solution of this equation can be written in the form 


¥ i aitm 
Py (x2) = (1— wt Th i age ——_— (4? —1)', [18.14] 
where the condition 
me ae [18.15] 


on the integer m must be satisfied. The formula is nor- 
malized such that 


7) ie | [18.16] 


By applying formula [18.14] to positive as well as nega- 
tive m, we obtain the relation 


PP = CinP5”, im = (— 1)” oar [18.17] 
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Therefore, the solution of the original differential equa- 
tion [18.9] is 
Yi m(9, p) = PP (cosd) etme , {18.18] 


The functions Y,,, (or the P? alone) are called tesseral or 
spherical harmonics. The nodal lines of their real and 
imaginary parts divide the surface of a sphere into four- 
sided regions (tesserae) of different signs, separated by par- 
allels of latitude and meridians of longitude. The functions 
Y,,=P', written simply as Y, and P,, are called zonal 
harmonics, because they divide the sphere into latitudinal 
zones of different signs by their nodal lines. With the 
normalization [18.16], the P,(7) are the Legendre polyno- 
mials: 


a1, Pa P, = $0? —}4, P. 


3 
35 15 3 
a 3 ct — Zoe in On 


The spherical harmonic Y,(0,q) is 1/r’ times a homo- 
geneous polynomial of the Ith degree in a, y, z which sat- 
isfies Laplace’s equation: 

_ A(2, y, 2) 
=. 


Yb, ~) Vii, = 0. [18.19] 


The Y,(#, y) can be expanded in terms of the linearly inde- 
pendent and orthogonal functions Y,,,(0, 9): 


Y, (A, ~) = y Oe Se (8, Q) ° [18.20] 


m=—t 


The value of the normalization integral is 


m ie re (t+m)! 2 


Using this value we obtain 


P7'(cos#) exp [img] 
/N™ V 27 


for the normalized spherical harmonics. 


Vind, p) = [18.22] 
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c. Solution of the radial differential equation 


With the condition [18.8], Eq. [18.6] becomes 


[18.23] 


— (rv) + Fle 11) — EO | e= 0. 


2mr? 


The eigenfunctions v satisfy the orthogonality relation 
foro dn 0 for EAH, V=t. 
o 


We see from the differential equation that there is for- 
mally an additional potential 


(d+ 1)h?2/2mr2 


for the radial motion in a central field. This expression is 
the analogue of the potential of the classical centrifugal 
force, 


a8 P 
A 2mr? 


(P = angular momentum) 
if we set 


PP=W1+1)h. 


This gives us an important relation between the quantum 
number / and the angular momentum |P| of the particle. 

In order for the radial differential equation [18.23] to 
be solved, the potential V(r) must be given explicitly. We 
choose as an example the Coulomb potential: 


Vr) =— 5 [18.24] 
with 
¢,=+ Ze (Z =atomic number), 


€,=—e (¢ =magnitude of the elementary charge) , 
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and obtain 


cae (E+ =) (rw 0) — CED poy = 0. [18.25] 


Our goal is the wave mechanical description of the sys- 
tem consisting of an electron plus a nucleus. For this 
purpose we should use the corresponding Schrédinger equa- 
tion with the 6 coordinates %,, ¥,, 21, %2) Ye, % of the two 
particles. Since the potential depends only on the relative 
coordinates, V=V(a%,—22, ¥:— Yo, 2%: — 2), this Schrodinger 
equation (in analogy to classical mechanics) can be split 
into an equation which describes the motion of the center 
of mass plus an equation which determines the relative 
motion of the two particles. The latter equation is formally 
the same as the equation for the motion of a particle of 
mass m moving in a potential V, where m is the reduced 
mass, 

1 it ne a 


; [18.26] 
m Melectron Maacieus 


By substituting the reduced mass m (which is practically 
equal to M,.cto. in this case) into [18.25] and solving this 
equation, the energy eigenvalues and the eigenfunctions 
for the relative motion of the electron and nucleus can be 
determined. 

For this purpose we introduce the dimensionless quanti- 
ties @ and e: 


r=; EEN DE [18.27] 
he : 
% = Fim (Bohr radius) 
[18.28] 
R- Z?e*m h? (energy corre- 
= 


2h? ~ 2mat sponding to a) 


In terms of these quantities we obtain 


d? 2”) _ 0 18.29 
So (o0) + (0+ 0? ee ; [18.29] 
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For small 0, we can immediately verify the solutions 
v~o', oro ; [18.30] 


for large g, by dropping terms of orders 1/9 and 1/97, we 
obtain 
»~ exPltev= =] 
e 


[18.31] 


One way of solving the differential equation [18.29] is by 
using @ power series as a trial solution; however, we want 
to profit from our knowledge of the hypergeometric func- 
tions. Taking account of the partial solutions [18.30] and 
[18.31], we assume 


v = op’ exp[-+t ovV/— €]-a(o) [18.32] 


as a trial solution, where w(@) is to remain finite for 900. 
If we compare the differential equation for w, 


d?w 
do? 


+(4 a/=e |e 


de 
24+ 2(+1)vV—é 
@ 


ai =0 [18.33] 


(first dividing by —4e) with the equation for the hyper- 
geometric function F(«,y,x) (see [17.51]}), 


27 
¢ +(2-1)5F ro, [18.34] 


da? x da 


we see that they are identical if we put 


il 
wot Fett, [18.35] 


y=2(1+1), [18.36] 
w=F20/—e; [18.37] 
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thus, w is given by 


w = constant 


P be ee 
From the power series expansion [17.54], 
a 
F(a, y, Ele, Onis ; 


it follows that the solution [18.32] actually does begin 
with o’ for small 9, as we conjectured in [18.30]. 

In order to carry the discussion further, we must distin- 
guish between «<0 and e>0, that is, between ,/—é real 
and imaginary. 


d. Case 1: e<0. The discrete energy spectrum 


We take the lower sign for ,/—é throughout; in [18.32] 
this means that we take the minus sign. The upper sign 
does not lead to a physically acceptable solution, because 
the normalization integral diverges. 

Using [18.35], [18.36], [18.37], and the radial quantum 
number n,, 


1 
t= —t=———— 1-1, [18.39] 


then for large 9 we can write the asymptotic formula [17.65] 
for F(a, y, #) as 


i 1 —s 
he a, cae Ee) 
—— = ee) - (ye) -t-1 
rae) ED Gy —ea 20 2) = : [18.40] 


Because of the second term, this solution approaches infin- 
ity exponentially for 9->0co; such behavior is incompa- 
tible with the physics of the problem. We must demand 


96 EIGENVALUE PROBLEMS =| Chap. 5 
that this term vanish; this means that we must require 
RO; Ay Besos [18.41] 

We then obtain the asymptotic expression 
o~ exp[— oV/— €]- 01-9) 1 [18.42] 


for the solution [18.32]. According to formula [17.57], the 
condition [18.41] reduces the hypergeometric function 
F(a, y, 2) to the simpler associated Laguerre polynomial, 


L™ ,,.= constant x F(— k, m+1, 2). [18.43] 


Therefore, according to the trial solution [18.32], the eigen- 
functions of the radial differential equation [18.29] are 


a Q Wiehe 2 0 
pe ay Oe ae ee) 


The normalization can be carried out in a way analogous 
to that for the plane harmonic oscillator. 


v = constant X eo’ exp 


1. The quantum numbers. Energy eigenvalues and degen- 
eracy. The principal quantum number n is defined by 


alt 


Its values are determined by the condition [18.41]: 
n=141,142,.... [18.45] 


These relations also imply that for a given 1 the angular 
momentum quantum number | can only take the values 


LO; 1,..0.5%—1. [18.46] 
The energy eigenvalues follow from [18.44] and [18.28]: 


1 Z76*m 


Ree — 
n= Dh? 


[18.47] 
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With these results we can write the famous Balmer formula: 


224, 
hy = E'— BE", or y= R(— aE —= 
2h 


It is worth noting that the eigenvalues EZ, do not depend 
on the angular momentum quantum number / at all. For 
a fixed 1, we have seen that there exist 2141 different 
spherical harmonics or eigenfunctions of the Schrédinger 
equation [18.1]. To this (21+1)-fold degeneracy, which is 
present for a completely arbitrary central field V(r), there 
must be added the n-fold degeneracy [18.46], which is re- 
lated to the Coulomb field. Thus, for the total degree of 
degeneracy g, of a state with principal quantum num- 
ber n, we find 


ital 


In = > (2141) = ‘ST (41) = m?. [18.48] 


1-0 


We did not find this degree of degeneracy in any of the 
cases that we solved previously. 

We note that the grounds tate (n=1, 1=0) is nondegen- 
erate. Furthermore, since ¢«-=—1, the solution reduces in 
this case to 


u = constant x ee. [18.49] 


2. Examples in which the degeneracy is removed. The va- 
lence electron of an alkali atom moves in a central field 
which, however, is not quite the same as the Coulomb 
field. This has as a consequence that the n-fold degeneracy 
due to the Coulomb field is removed, and that the nth 
hydrogen-like level of the valence electron splits into 
levels. We have a further example if there is an external 
magnetic field. The magnetic field alters the character of 
the original central field completely, with the result that the 
mth level splits into n? levels. 
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e. Case 2: e>0. The continuous energy spectrum 


We come to the case which describes a “‘free’’ particle. 
With [18.38] and the assumption 


f/—e=+ive (positive root!) , [18.50] 


t 


we can write the solution [18.32] as 


v(e) = constant 
_ i 
U ante y —_. 
Xo exp[ iviolR (Ss 
where we again restrict ourselves to the “lower” sign of 
+./é. Introduction of the wave number k of the free 


particle, 


oe ae, Zio?) , [18.51] 


kh = V2mE, [18.52] 
leads to the relations (see [18.27]) 


- \/E 
ka, = 4/e = Vz 
[18.53] 


TE nw \[? 
\/£0 = a r= kr (oe 


with which the solution [18.51] can be written in the 
somewhat clearer form 


v(r) = constant 
x (ker)! ee F (= a OT 40, 2itr] . [18.54] 
0 


As a matter of fact, this solution is actually real. To prove 
this, we separate F into two parts, 


De 9 sie oa Ua [18.55] 
and use the integral representation [17.59], 
n= P(2l a 2) 
2 270 


x fertinas sa — Qikr) ~(skeae)-t-14¢ , [18.56] 


Oo: 
C3 
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for P\. If, for example, we substitute 
t—>t’— 2ikr [18.57] 


in the expression for F,, it follows that the path of inte- 
gration C, is the same as C,, from which 


6 — (ener, )* [18.58] 
follows immediately. Thus, we can write 
etkr Ff — e-thr Fr + (e-#r Fr, )* ; 


which proves the assertion. 

Now, we want to investigate the asymptotic behavior of 
the solution [18.54]. Using [18.58], the asymptotic formula 
for F,,[17.64], and the relations[18.35] to [18.37 ]and [18.53], 
the solution for large r can be written in the form 


v(r) = constant x (kr)! e-**" F ~ constant 


I(21 4+ 2) hie cai 
S—— SS kr)* etter (Qo ker) @/kae)—i-1 
aoa 
+ complex conjugate . [18.59] 
Using 
2ikr = exp hog AL iF 
and 


(Qikr)(a)-1-1 — 9-1-1 (kp)-"-1 exp = log (2kn) 
0 


aot et 
x exp I+ F(5-1-1)| ‘ 


we can also write this expression in the form 


(21-1)! 
Ti (é/ka,) + 1 1] 


1 a LS 
. i 227 | 
xX exp ier + er log (2kr) —4 5 (i+ } 


v(r) = constant x 


1 4 ; 
Apr =P - 5 = +complex conjugate. [18.60] 
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With the relation 


r( +141) =|Flexplia(t kao) ] [18.61] 
0 

(where o is the phase of the I’ function), we finally obtain 
for kr>1 


(21+1)! exp[—2/2ka,] 1 


v,(r) = constant x ot \I[(i/kay) + 1+ 1)| kr 


x COS jer ae {log (2kr) aa > (t al) all ka.) . [18.62] 
0 


Essentially, this yields a spherical wave for the solution, 
which, however, contains a logarithmic correction to the 
phase. The effect of the correction is to make the phase 
of the wave function for large r change more rapidly than kr. 
This correction, which is typical of the Coulomb field, can 
also be justified classically. 

For Z=0, that is, for 1/a)=0 (see [18.28]), the force- 
free solution follows from [18.62] (see also [21.10]). 

We shall show in an exercise (see Sec. 46) that normal- 
ization according to 


fousrint dr = 6(k’— k) [18.63] 


0 


leads to the normalized asymptotic formula 


: 1 
x sin {er + — log (2kr) — he ee o(l, kas) : [18.64] 
kay 2 


Remark: As is evident from the above formulas, the states 
of positive energy (e>0) are infinitely degenerate. This 
is because, for a given value of k, the quantum number 1 
can independently assume all values between 0 and co 


Sec. 18 | THE HYDROGEN ATOM 101 
f. Solution of the wave equation in parabolic coordinates 


The wave equation for the hydrogen atom can also be 
separated and solved in parabolic coordinates. We want 
to obtain this solution and show, as is to be expected on 
physical grounds, that the same energy eigenvalues and 
the same degeneracy result as with the spherical coordi- 
nates. Further, we shall also treat the scattering problem 
with this solution and derive a solution corresponding to 
formula [18.62] whose characteristics are more evident. 

The connection between the parabolic coordinates A,, A2, 


Figure 18.2 


g and the spherical polar coordinates (Fig. 18.2) is given by 
Ay=r+e=r(1+ cos) 
A, =r—2=r(1— cos#) +. [18.65] 
aa 
In addition, we have 
r = $(A,+ A) 


z Se A,) ‘ [18.66] 
o! = r2 — 22 = Avas 


® See A, SOMMERFELD, ‘‘Uber die Beugung und Bremsung von Elektronen,”’ 
Ann. Physik 11, 268 (1931). 
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This time we immediately write the wave equation [18.1] 
with the Coulomb potential [18.24], using the dimension- 
less coordinates «= EH/E, and 9o=r/a, given by [18.27]: 


Vu + (« a 4 u=0. [18.67] 


In the following, until the treatment of the continuous 
spectrum, we also use dimensionless parabolic coordinates; 
these follow from formulas [18.65] and [18.66] by making 
the substitutions r>r/a,, z—z/a). For the sake of sim- 
plicity, we also call these dimensionless quantities J, and A,. 


g. Separation of the wave equation [18.67] in parabolic co- 
ordinates 


In this case the square of the element of length is (dimen- 
sionless quantities!) 


ds? = do? + dz?+ o?dq@? 
=Ti 1. (A, dA, + A.dA,)? + 7 (dA — das)? + Aad? 


A, Ads, Ga 1 \ 
= 4 = Bious : 
m i, Ts “a a +7] dy [18.68] 
As it must be with saree coordinates, there are no 
mixed terms. It can be shown’ that the Laplacian operator 


in orthogonal, curvilinear coordinates 7,, v,, 2, is of the 
general form 


v1 1 90 ("se du ; 1 0 fese, du 
gee 04, \ €, Ob) eyenen oa, \ €, ee 


, 1 9d fee, du 
-f ae, 
C1 Cos 0%5,\ Cs dX 


a ; [18.69] 
with 


3 
ds? = > exdai., dV =e, ¢,¢,d2,dx,da;. [18.70] 
k=1 


7 See W. Pautt, Lectures in Physics: Electrodynamics (M.1.T. Press, Cam- 
bridge, Mass., 1972). 
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In our case this leads to 


| cep) (2 8u 
2 eo) je ee = = SS 
ih {sr (A a a @ _ 


Al 1\ 0? 
—|—+—]—_}. WS 
4 ( . } =| 


With this result we can write the wave equation [18.67] 
in parabolic coordinates as 


P) ou Lt 1 \ de 
Bean) on (asa) +3 a 
+. {et +As)+ fu =O,  fisi72) 


This equation can be separated with the trial solution 
UW = f,(A1)*fo(Ag) EF; [18.73] 


the result is 


a (aae)+(-F5 $4 Ve )ine [18.74] 


dA, 4 
df, of (754% i 
a (sgt) + ee sf\h=o, [18.75] 


where the parameter # is the arbitrary separation con- 
stant. These two equations have the form of Eq. [17.16] 
for the plane harmonic oscillator, which we now write as 


y+ |—F+ (m+ ™F2)— Fly =o, i=1,2 [18.76] 


(A/4 >n;+(m+1)/2, k>n,). We found 
f =constant x @”? e-*!? Din, (&) 

= constant X 2”? e-2/2 F(— n;, m +1, @) [18.77] 

to be the only useful solution of this equation (see [17.27])- 


The numbers n,, m,, and m are quantum numbers asso- 
ciated with the parabolic coordinates /,, A,, and @. 
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h. Discrete energy spectrum (e€<0) 


We can bring [18.76] into agreement with [18.74] and 
[18.75] by setting 


e=Apvp/—eE (/— é> 0) [18.78] 

and choosing f such that 
i _—_ a 
s+ A)=v=e(m, +42) 
1 1 , [18.79] 
a p)=V—e (n+ 7] (n; = 0, 1, 2,...) | 
From this, 

l= V/— e(n, + ,.+ m+1) [18.80] 


follows immediately. With 
N=n+n,+m+1 (c=, 2 nc) [18.81] 


this can also be written in the form 
3 SS [18.82] 


which corresponds to the expression [18.47] already found 
for the energy eigenvalues. In addition we obtain 


Ni — Ns, 


B= V/— &(n,— 2) = [18.83] 


The solution [18.73] can now be written in detail (a —>A,/n): 


A ae A, Ay hy mem A, 
on n2 mtn, n 


< Lean (*) etime 18,84] 


u = constant X exp 


Remark: The separation of the wave equation in para- 
bolic coordinates is also possible when there is an external 
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field F, which adds a term 
constant x (Fz) 


to the potential V in [18.1]. The only consequence of this 
is that the separated Eqs. [18.74] and [18.75] have addi- 
tional terms 


+constant x Aif,F and —constantx Af, F, 


respectively. If F is small, these equations can be solved 
approximately by means of the Born approximation (see 
Sec. 24). An interesting result is that the energy eigenvalues 
have the form 


1 
is + constant x (n, —n,.)nF 


in first approximation (first-order Stark effect).® 


1. Continuous spectrum (e> 0) 


We now return to the usual units; that is, we use aA, 
instead of A;. As with [18.52], we introduce the wave 
number & and, with the help of [18.53], 


\/ Sees = Fhe, [18.85] 


(the upper sign is for index 1), we write the condition [18.79] 
in the form 


5 (1+) =— thay (0, +5) 


2 2 
[18.86] 
JI : m+1 
Z(1—f)= thao, +™F2) 
Addition of the two equations yields 
B 
v, —V, = — [18.87] 


® See, for example, A. SOMMERFELD, Atombau und Spektrallinien, vol. 2. 
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where, in contrast to the n;, the v; are not integers. In 
this notation, the differential equation [18.72] is 


d 12] d (2 du vat | ) 
SA 3A, )| | ae oe ae 


2 fe Gene =a u=0, [18.88] 


and its solution is (e—-Fikd,) 


u = constant x exp ie ral (k2A,A,)"™?2 F(a, m+1,—ikd,) 


0 


xF (2 aE = m+, ik’, | etme [18.89] 


In this formula we have replaced the still undetermined 
v, by —a. 

Remark: All of our formulas and expressions are written 
for positive a), that is, for attractive forces. In the case 
of repulsive forces (which only makes sense for e>0), 
we must make the substitution a,>—a,. 


Chapter 6. Collision Processes 


We concern ourselves here with the scattering of particles 
by particles, without taking spin interactions into consid- 
eration. The solution of this problem follows from the for- 
mulas which we obtained in the last chapter for the case 
of a continuous energy spectrum (e> 0). 

The starting point for the present considerations is the 
solution [18.89]. The incident particle current will be rep- 
resented by a plane wave propagating in the positive z di- 
rection; of course, the wave will be symmetric about the 
2 axis, and we can therefore restrict ourselves to the case 
m=0. For the special case of «=0, we shall show that 
the solution [18.89], 


P (= i ikis| |i 


a) 
ka 


u = constant x exp = (A, — A,) 


can be written asymptotically as an incident plane wave 
plus the scattered wave. For the general case of «+0, 
the solution contains terms which represent a plane wave 
propagating in the negative z direction and an incoming 
spherical wave; this follows immediately from the asymp- 
totic formula [17.65] for the hypergeometric function F’. 
Therefore, this case is incompatible with the physical prob- 
lem as we have stated it. 

Remark: In the relativistic case it is not possible to sep- 
arate the wave equation in parabolic coordinates. In par- 
ticular, this simple decomposition of the asymptotic solu- 
tion is also not valid. 
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19. ASYMPTOTIC SOLUTION OF THE SCATTERING PROBLEM 


The asymptotic formula [17.65] for the hypergeometric 
function F(t/ka, 1, ikA,) is 


~~ ; —tlkag 1 
se £(1—i/ka) iain {a+ Fite| 
1 y ; —1+i/kay - 
+ Fea CREA) TaeBe [19si} 


in this expression we have used the first approximation 
given in [17.62]. Using [18.65] and [18.66], we again intro- 
duce polar coordinates r and #?: 


Ay = + + z 
2, = r(1— cos?) = 2r sin? 40 ‘ [19.2] 
Aim A, = 22 


Written out in detail, the solution [I] is now 


exp[— 2/2kay| ee - 
A — 9h 63 Ta =a i/kay) exp i kez = eq ee sin2 10)) 


ml i _ ,T'(1—ifkay) 
2kr(kay)* sin? $3 I(i/kay) 
exp [i(kr + kag" log (2kr sin? $8)) | 
- 2kr sin? 49 i 


In obtaining this expression we have used 


(— iy exp|- | 
[19.4] 


(kA,)~#*4e— exp|— —_log (2kr sin? 3%) 
kay 
In order to further simplify the expression, we use the 
formula 


ras t/kay) _ oe I'(1 — ifkay) 
iT(ilkay) kay (1+ t/kay) 


1 ; 
= ba, oP [— 2io(0, kay)], [19.5] 
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which follows straightforwardly from the functional equation 


a a i 
Fi1+ —)j=—T|— 
( +z] i (=) ie) 
and the definition of o(l, ka)) given in [18.61]. With the 
normalization 


exp[—2/2ka| _ 


T(i—ijffeas) ? ee 


the asymptotic solution [19.3] can now be written in the 
form 


= a fon 21 a 
ie exp | (ie 2 i gi6 g(2kr sin 10))) (1 Dheatr sin?) d 
__ exp fifkr + ka a sin? 40) — 20(0, kay) Be 


IG), 
2k?ayr sin? 43 (ie 
By assumption this formula is only valid for 
kA, = 2krsin?40 >1 ; [19.9] 


that is, it is valid in all of space except for a parabolic 
region defined by [19.9] (Fig. 19.1). 


Figure 19.1 


To lowest order, the solution [19.8] is of the form 


exp a 


u~exp[ikz] + f(#) [19.10] 


that is, w~incident plane wave + outgoing spherical wave. 
In the case of a Coulomb field, with which we have cal- 
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culated so far, the following substitutions must be made 
for the phases: 
i 


reputed sin? 47) 


ke—> kz —- 
[19.11] 


kr—> kr + —- log (2kr sin? 4) 
0 


This logarithmic phase correction does not occur if the 
potential falls off more rapidly, 


lim V(r):r =0. 


In an exercise (see Sec. 46) we shall note that another 
asymptotic representation of the wave function, one for a 
fixed value of J, is possible for such a potential: 


nim SOREN sin (ir —13 +61). [19.12] 


In formula [18.62] we have already derived the analogous 
asymptotic expression for the Coulomb potential, which 
follows from [19.12] with the substitution [19.11]. 


20. THE SCATTERING CROSS SECTION. THE RUTHERFORD SCAT- 
TERING FORMULA 


Using the solution given by [19.8] and [19.10], we can 
calculate the differential scattering cross section dQ. The 
definition of dQ is 


(ee of particles scattered eh 

dQ _ts__\ solid angle dQ per unit time 

"4,  (mumber of incident page Pe" : 
( unit area per unit time 


[20.1] 


If the solution of the scattering problem is of the form 
[19.10], then we have 


dQ =|f(9)|2dQ. [20.2] 


We prove this formula by calculating the ratio of the seat- 
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tered-particle current to the incident-particle current. 
From the general formula for the particle current den- 
sity [7.10], 


x h 
i= Omi (y* grady—y grady*) , 


we obtain 
hk 


— [20.3] 


ty = 
m 


for the (incident) plane wave; for the wave scattered into 
solid angle dQ, we immediately obtain 


hk 
= aren't dQ [20.4] 


(remember kr>>1 and grad~d/dr). Thus, formula [20.2] 
is proved. 
From [20.2] we obtain 
1 


dQ = 4k*a2 int} dQ . [20.5] 


for the solution [19.8]. With 


a mo  2mE ht 2K 
~—\ he COUR eam ~—«C ee’ 


this yields the famous Rutherford scattering formula, 


[20.6] 


Ze Ze! 


= ———___—_. = ——_——__ dQ. 20.7 
sad 16H? sin*4d 4m?vt sinthd) (age 


This expression, as well as formula [19.8], can also be 
derived classically; ! it is a result of the particular Cou- 
lomb potential. With other potentials it is generally not 
possible to obtain the same results classically and wave 
mechanically.? 


2W. Gorpon, Z.'Physik 48, 188 (1928). 
2 Th. SEXL, Z. Physik 67, 766 (1931). 
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21. SOLUTION OF THE FORCE-FREE WAVE EQUATION 


Now we want to try to solve the scattering problem di- 
rectly, without reference to the general asymptotic for- 
mula [17.65]. As a result we shall understand the problem 
somewhat better. 

Toward this end, we first solve the force-free wave 
equation 

V2u + k?u = 0 [21a 


in polar coordinates. With the usual trial solution [18.4], 
u= vr) ¥ (8,9) , 


we separate this differential equation into an angle-de- 
pendent part and a part depending only on r. The part 


depending on r is (see [18.6], with V =0) 
2 U(1+1) 
r2 


== n=0, [21.2] 


oar) + (k— 
or, written in terms of the dimensionless coordinate 9 =kr, 


utd) 
= v,) + (1 z Jo=o. [21.3] 


By induction (J>I+1) we obtain 


M+ =—@! = (o-'v:) = — x sad ate [21.4] 
Using this, we can write the solution of 5G in the form 
v,= 0'(—1)! feoke , [21.5] 

ede 


From [21.3] we immediately obtain two independent solu- 
tions for 0,_,: 


amen exp [ie] outgoing spherical 
io wave 
a 21.6 
pes exp[— io] incoming spherical 


10 wave 
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From these we obtain, for example, the solution 


: 
po = 5 (60+ 6) = 2 : [21.7] 


eS 


which is the only solution regular at g@=0. With [21.5] 
we can now write 


Q d \ig 
£3 — 91) oa) ge [21.8] 


w= $(éP + éP) 
d 1 
=o'(—1)! (5) Wot. [21.9] 
yi= (EP— EP) /26 


We now want to write this solution in a somewhat sim- 
pler form for the two special cases e>>1 and 9 <1. For 
o>1, all that we need differentiate in [21.8] is ¢*, because 
the other terms are of higher order in 1/g: 


,exPlie] a vive 
oe As 


__ exp {ile —U(a/2)]} 
a a . 


P= (-1)'(+4) 


With this we immediately obtain 


__ sin (e —Ia/2) 
—= 7 


Vr [21.10] 
On the other hand, for @<1, all that we need differen- 
tiate in [21.8] is 1/0. The result is 


o? 
Pe TX8X5X...X (2E4 1) 


aa xb x... x (21— 1) 
i a a | 


[21.11] 
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Connection with the cylinder functions 
(09) 
There is the following relation between the solutions ¢” 


and y, and the cylinder functions: 


@ a: tt (0), [21.12] 


y= 2 Tis3 (0) - [21.13] 


The functions H® +; are the Hankel functions of the first 
and second kinds, and J,,, is the Bessel function. It is 
well known that these functions have the following integral 


Figure 21.1 


representations (see Fig. 21.1): 


1 
ike) = On fexp [to cosw] exp in ( — | dw , 


Wo 


27 
Wr 
Wa 


a il F 7 
HY (0) = — |exp[ie cosw] exp lim ( ~3) dw. 


Substituting [21.12] or [21.13] into [21.3] does indeed lead 
to the differential equation for the cylinder functions: 


es (ee 
do? 0 de 0? iio . 121.0 
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22. EXPANSION OF A PLANE WAVE IN LEGENDRE POLYNOMIALS 


We want to expand the plane wave e’** in terms of Le- 
gendre polynomials (#=cos#): 


eke —explikr cost |=¢e= > fi(o)Pi(v). [22.1] 


The coefficients f,(@) are to be determined. In order to 
do this we need the integral representation 


+1 
ui 1 @ ‘1 tox 
S,= i (2) afe ez (1— x?)'da = pile) 4 [22.2] 
-1 
which we want to prove by deriving a recursion formula 
analogous to [21.4]. We have 
+1 
1 ee— ee ging 
=— ton = SO Ss = A 
So fe dx Dio F Po( 0) 


=i 
We can write S,, in terms of 8;: 


cian 


all o\ 1 iy 
Sa = G1)! (2) sfe (1— #7?) az 


ih 
oil 


a (5) 5 | Ge e+ (20) de 


2 


+1 


fot d Lj 
ee ||| ee 2(1 — 7?)'d 
(8) 3 |r" Ce 


-1 
d d l 
= — 9'— (9 8,) = —— 8+ &.-. 
. do (o-*S:) do Geer F 
Thus, since this is the same as the recursion formula [21.4], 


y, is indeed equal to the integral S,. 
We need the following well-known relations for the Le- 
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gendre polynomials (see [18.14], [18.16], [18.21]): 


—1)'/a\ 
Pca = P(e) a-e, 


+1 
freee) da = d.y 


ak 


S41’ P,1)=1. 


Using these relations, we immediately obtain 


ee +1 
Weiss 
file) —[eeP eae =o (ig)! fewer — x?)t dx 


Sil =1 


2 
2i+1 
from [22.1]; the last expression is the result of n-fold par- 
tial integration. By comparing with [22.2] we arrive at 


a 


t'2y,(0) = 
and with this we obtain the desired expansion: 
exp[io cosB] = > (21+ 1)i'p,(o) Pi(eos®) . [22.3] 
’ u 
23. SOLUTION OF THE SCHRODINGER EQUATION WITH AN 


ARBITRARY CENTRAL POTENTIAL 


In the wave equation 


2 
V2u + ku — = Virju=0, [23.1] 
we no longer use the Coulomb potential, which has special 
properties; instead, we simply require 
limrV(r) = 0, [23.2] 


that is, that the potential V(r) fall off faster than the 
Coulomb potential. With the usual separation, 


U = 0,(0) Yim(P, Y) » [23.3] 
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and with the abbreviations 


2m 


kr=0, =>V(r)= U(r), [23.4] 
the equation 
it de Wi+1) U(r) 
0 do? Ov) + {1— 0 Vv, ke —— 0 [23.5] 


follows from [23.1]. 


a. Method of partial waves 


We know two special cases of the solution of [23.5]: 
o<1: v,e)~o' if U(0) is finite [23.6] 


this can be verified, for example, by a power series trial 
solution (see [18.30])); 


Oe lie (0) ~—sin (0-13 +6, (i) [23.7 ] 
(see, for example, [19.12]). The latter, asymptotic repre- 
sentation was introduced in an exercise (see Sec. 46) in 
which 6, was defined to be zero for U=0. In the follow- 
ing, the normalization of v, in this formula will be chosen 
such that c=1 (this normalization is different from that 
in Section 46, where ¢ equals k\/2/z). Here we want to set 
up a different representation which is more practical for 
our scattering problem.? In analogy to formula [19.10], 
we take 


u = exp[io cos?) + f(a) ——— ae ao [23.8] 


This trial solution satisfies two conditions: that the solu- 
tion not contain an incoming spherical wave, and that it 
only consist of the incident plane wave at large distances 


2? This rather arbitrary normalization is meaningless here, because it only 
changes the coefficients in [23.9], 
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from the scattering center. In order to determine /(6), we 
now require that the asymptotic formula [23.8] agree with 
the general solution of the differential equation [23.1] at 
large distances, and we write this solution in the form 


u = > 4,0,(0) P(cos®) [23.9] 


(py does not appear because the problem is axially sym- 
metric). By expanding the plane wave in [23.8] according 
to [22.3] and, in [23.9], using the asymptotic representa- 
tion [23.7] for v,(@), normalized with c=1, we obtain 


> (20+ 1)é' sin (0 —Il- 3) Pa(cosd) + kf(P)ere 


= Da,sin (0-1 + 0449] P,(cos#) [23.10] 


by equating the two solutions. We have also used [21.10] 
here. Now, we write this equation in a different form: 


¥ (U4 18 = {ot ( i) — e843, P,(cos) + kf (9) e 
I 
= = Las tere eft (— 4)t— en e715} P (cos) . [23.11] 


Setting the coefficients of ¢'* and e equal to zero, we 
obtain 


a,= (21+1)itexp[-+ 76,], [23.12] 
arr a (21 + 1)(exp[2i6,] — 1) P,(cos#) 
al} 
=a (21+ 1) sind, exp[id,] P,(cos#) . [23ers] 


b. The differential and total scattering cross sections 
The differential cross section [20.2] can now be written as 


dQ = |f(8) |?d2 


i : : 
=a | > (21+ 1) sind, exp[t6,] P,(cos#) |?}dQ. [23.14] 


Sec. 23 | SOLUTION OF THE SCHRODINGER EQUATION 119 


We obtain the total cross section by integrating over solid 
angle (d2 =sind dd dq): 


Q = [dQ = 2x] |f(F) |?sin ddd. 


All of the cross terms vanish because of the orthogonality 
of the Legendre polynomials. Thus, using the normalization 


+1 
integral | Pi(x)dx =2/(21+1), we can immediately write 
t 
= 
= a 21 +-1) sin?d,(k 23.15 


c. Calculation of the phase shift 6,(k) 


We now want to calculate the phase shift 6,(k) between 
the unperturbed wave function, [21.10], and the wave 
function perturbed by a potential, [23.7]. We multiply the 
differential equation 


2 (1+1) 
2S ry) + (i 2 )m=o [23.16] 
by —r*v, and add to it Eq. [23.5] multiplied by 0’y,: 
d d d 
& Lorwer 5, (r00) = (ron) 5 (rah = Uir)rtpie 


d 
{trv 5 (re) — (rv,) a, iro} 


R 
—0 =| D¢ryrepinr [23.17 | 
r=R 
a 
Using the trigonometric identity 
sin a-cos(a-+6) — sin(a +6) cosa = — sind, 
we obtain the exact formula 


sind ,(k) = —k{ U(r)r?y,(kr)v(kr) dr, [23.18] 


0) 
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in which we unfortunately do not know the v,. In certain 
cases we can circumvent this difficulty by substituting yp, 
for v, as a first approximation: 


sind, (k) = — kf Un widr [23.19] 


24. THE BORN APPROXIMATION 


Again, we let the solution of Eq. [23.1], 


V%u-+k?u= U(r)-u, where U(r)= “ VG). [24.1] 
have the form 
u = exp[ikrcos?] + u,+..., [24.2] 


where the perturbation u,, which represents the scattered 
wave, is to be small compared with the incident plane wave. 
Substituting, we obtain 


V2u, + k2u, = U(r) exp[ikr cosd] [24.3] 


to lowest order; the term U(r)u, has been neglected. Of 
course, we could continue this method and write 


V2u.+ k®*u,= U(r)-u,, ete. 


However, it turns out that the Born approximation is 
practical only if the lowest approximation is sufficient. 
The solution of Eq. [24.3] is worked out in the theory 
of the Hertz dipole* (Fig. 24.1); it is 
1 { U(r’) expfikr’ cos 3’ 
(=) ae eee 


Ud at 
fe SY bere yay”. [24.4] 


It contains the radiation boundary condition which is also 


“See W. PauLi, Lectures in Physics: Electrodynamics (M.I.T. Press, Cam- 
bridge, Mass., 1972). 
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needed here. We do not need this solution in all general- 
ity; we restrict ourselves to two special cases. 


Figure 24.1 


P: field point; @: source point; O: origin 


We first investigate what happens when P—0: 


ee 1)] 


j= — 
U1 ( ) An yx! 


dv’. [24.5] 
If we integrate over solid angle (dV’=r'? sin# dr' dd? dg’) 
and omit the primes, we obtain 


u,(0) = — 7 | Oe) sinkr exp[ikr]dr. [24.6] 


This expression represents the scattered wave at the origin. 
By assumption, its amplitude is much less than 1 (we have 
normalized the incident wave to 1), which will certainly 
be true if 


Al U(r) |dr <1. [24.7] 


This is a sufficient condition for the applicability of per- 


turbation theory.® 
Now, let us investigate the behavior of the solution [24.4] 


5 Therefore, the Born approximation is to be used if the kinetic energy 
of the incident particle is large compared t» the interaction energy. It 
thus complements the method of partial waves, which is valid at low energies. 
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at large distances: 


r >>d=range of the potential. [24.8] 
With . 
; r’=r—r'cosé, [24.9] 
we obtain 


exp[ikr] 1 
r 47 


foe exp [ikr’ (cos #’— cos6)] dV’, 


or, by expanding the plane waves in Legendre polynomials 
according to [22.3], 


= exp [tkr] i \ lady! pt iN (Nit 
=— a rakes }r2dr Da (20+ 1) p,(kr’)(— 7) 
x (21 +1) pr(kr’) i P,(cos) P(cos9) dQ’. —‘ [24.10] 


Because of the orthogonality of the Legendre polynomials 
(also for the angles 3’, 6 of Fig. 24.1), only the term with 
V=l is nonzero: 


= expel U(r')r? dr’ > (2141)? 


Py, t 
: dQ’ 
x {y,(kr’)}? re P,(cos8’) P;(cos8) . [24.11] 
Using the relation 
dQ’ 
P(cos#) = (21+ 1)[P.(cos 9) P,(cos®) es [24.12] 


which holds for the angles of Fig. 24.1 we then have 


Tee sie > (21+ 1) P,(cos 9) 


U 
ro) 


x [ OU yet tery dr’. [24.13] 


u 
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This result matches the one which we already had, except 
that a factor ¢' has been neglected here, because 6, is 


very small (see [23.8], [23.13], and [23.19]). 


25. SCATTERING OF LOW-ENERGY PARTICLES 


To conclude our treatment of collision theory, we briefly 
consider a calculational method that uses the following two 
assumptions: (1) U(r)~0 for rza (a need not be pre- 
cisely defined); (2) ka <1, that is, the energy of the inci- 
dent particle is srnall. For r2 a, the term U/(r)v, in the 
differential equation 


1 ‘ rit (ie U(1 4-1) 


TG (2 


—\o= U(r), [25.1] 


can be neglected, and we have the two independent solu- 
tions of the force-free equation ([21.5]-[21.9]), 


qa) (2) 
ge aa [25.2] 
and 
Gee coe 
25.3 
a [25.3] 


The asymptotic forms of these solutions for kr 1 are 
y~ - sin (ur = 13) ; [25.4] 
= 008 (ir 13) [25.5] 


The general solution of [25.1] has the form 
= Aypkr) + Byilkr) . [25.6] 


Comparison with the asymptotic formula [23.7], 


ee ~ sin (ir—1%+6,), 25.7] 
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yields 
A=cosd;,  B=—sind,. [25.8] 


In contrast with what we did before (r>>a, kr>1), we 
now assume ka<kr<1. Because of k<1/r, [25.1] re- 
duces further to 


d? ot) 
~~ 


apa (re) (rv;)=0. [25.9] 


The differential equation has the solution 
Yr t r -—I-1 , 
eC, (=) + 5 (=) (€,/e, independent of k). [25.10] 


On the other hand, the general solution [25.6] with the 
(constant!) coefficients [25.8] is still valid; we just have 
to substitute the approximations [21.11], as a result of 
kr<1. Then we obtain 

(kr)? 
1X3 x. See) 
41x3x...x (21—1) sind,:(kr)-*, [25.11] 


0, = 080," p(kr) — sind, - (kr) = cos, 


Comparison of the two solutions yields 
1X8xX...x (2l—1)X1X3xX... x (21+ 1)(ka)-2*" tand, 


c 
=—= 4 (free parameter, independent of k), 


Cy 
or 
(ka)244 

a . 

and; = aan a2l 1) Ora)" [25.12] 
since 6, is small, we can write 

sind, ae (Kea)?! 
i {ix3x..x@l-py@l+1) t8l 


This relation shows that the amplitude of the Ith partial 
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wave in [23.13] is of the order of magnitude of (ka)?!. 
If this expression is substituted in the formula for the total 
cross section [23.15], 


sin?d 
U 
then we see, for example, that the isotropic part of the 
scattering (1=0, s wave) is independent of the energy of 
the incident particle at low energies. 


Chapter 7. Approximate Methods 
for Solving the Wave Equation 


26. EIGENVALUE PROBLEM OF A PARTICLE IN A UNIFORM FIELD 
a. The Airy function 


We imagine a particle in a uniform field: 


Exon = — F'-q. [26.1] 
The one-dimensional wave equation, 
d2u . 2m 
at ee + #E)u=0, [26.2] 
can, in terms of new variables 
pees 
2mF 
mt -// “ae 2 ao), = REG, [26.3] 
be written in the simpler form 
ale + 7u = 26.4 
da? — : [ . ] 


We solve this equation, which is familiar from the theory 
of diffraction, with the trial solution 


u =| e**f(t) dt. [26.5] 


0 


ut=[enrpat ; 


LU ={(5 e*) jdt -[5 (ertf) dt — for dt 


126 


Using 
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and [26.4], we obtain 


ae 
j 


as long as the path of integration C is chosen such that 


{* 
ef—f'=—0, Le f =constant x exp E [26.6] 


3 ? 


fae f)dt=0. [26.7] 


oC 


Now, the solution can be written, with the usual normali- 


zation, as 
y= eee lot Le [26.8] 
p= — > [exP 3 dt - : 
Cr 


Condition [26.7] is satisfied and convergence of the solution 
is guaranteed if, for |t}|+co, we require 


cos 39 < 0 GScle*), 
that is 
3 
22n + 5 <3P< ++ 22n. 


For example, the good intervals are 


ee 
Soy” NCE 6 I eG? GE 


The regions corresponding to these intervals are shaded 


C, 


Figure 26.1 
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in Fig. 26.1. It is evident that 


ie 


that is, there exist two independent particular integrals. 
The solutions wu, are called the Airy functions and they are 
denoted by A”. They are well known to us from the 
theory of diffraction.1 If # is real, then A” and A® are 
complex conjugates. Later we shall also need the average, 


A =}(A®+4+ A®), [26.9] 

Finally, let it be noted that the Airy functions are related 
to the cylinder functions: 

exp[tz/6] _ oo 

V/3 


By means of the saddle point method, we shall now 
derive asymptotic expressions from the integral [26.8]. 


A® —_ CO” ./e- HY (5 a) : Cos 


b. The saddle point method 


In order to evaluate the integral 


Jes [f(t)] dt, [26.10] 


we first expand f(t) about the ‘saddle point” ¢, which is 
defined by 


f'(t,) = 0: 
the expression that we use is 
(¢—%,)? 


f(t) = F(to) + 


gf’ (to) + - [26.11] 


In order that Re(f) fall off as rapidly as possible, we choose 


‘See W. Paul, Lectures in Physics: Optics and the Theory of Electrons 
(M.I.T. Press, Cambridge, Mass., 1972). 
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the path of integration to pass through ¢, (following the 
method of steepest descent of Cauchy and Riemann); then 
we have Im(f) constant. Let 


Fi (to) = |f"(to) | e** . [26.12] 
Because the imaginary part is constant, we then have 


(t — ty)? = — 0? e** 


t—t,= 9 exp i] 


dt = do expli7 >", 
and, as a result, 
f(t) = f(to) — 3 |f" (to) 0? - [26.13] 


Let us be bold and omit the higher terms of the expansion 
[26.11] when we evaluate the integral [26.10]. In this ap- 
proximation we then obtain 


| exp [f(t)]dté=exp [f(t.)]exp [ : . z | exp [— 3/f"(t) |e? ]de 
pte ee |/ 27 
=> exp [f(to) ] exp [ 9 eo) | ) 
and, finally, 


extn [f( t)] di = aia exp [f( (to) )] = "( is) [26.14] 


0 


with t, defined by f'(t)=0. According to [26.12], an ap- 
propriate sign still must be chosen for the square root. 
For an exact justification of the above calculations see, 
for example, Courant and Hilbert, Methods of Mathematical 


Physics. 
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c. Approximate solution for a particle in a homogeneous field 


First, we must determine {: 
eB ' 2 " 2 
ft=at+s, f(t) =e7+t=0, f"(to) = 2tr, t= — 2; 


tj=ti/2, for. 220, [26.15] 
ij=-+4/|%|, “for a@=10. [26.16] 


For the case «<0, because of 
i(t) =— $= +32, 


there are two terms in the solution [26.8]: 


AP = — Fe jv|-texp[+ 2|a|?] 


a ygleltexpl— lal, [26.17] 
AP = + Fole|texpl+ Fl2l4] 
inal 
Cet amen ati [26.18] 


This is because the path of integration goes over two saddle 
points, as shown in Fig. 26.2. Of course, the right-hand 
saddle is only half covered and, for that reason, its con- 


Path of integration: Im(/) = constant. 


x : saddle point 
Figure 26.2 
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tribution to the solutions is only half of what it would 
have been otherwise. 

In the case #>0, we have only one saddle point on each 
of the paths of integration C, and C, (see Fig. 26.3). The 


Figure 26.3 


solutions (with the famous exp[%(z/4)]) are: 
A® = : xtexp[i (2 2t — }x)], [26.19] 
Vi 
(2) 1 
A® — —x* exp[—i(zai—}x2)]). [26.20] 
4/7 
Therefore, the average values [26.9] of the two functions are 


2< 0: A=55 |v|-texp[— #|a|4], [26.21] 


VX 
aft 


at cos (¥ x? — 42). [26.22] 

It must be emphasized again that this solution represents 
an approximation. The Taylor expansion [26.11] is only 
appropriate in the neighborhood of the saddle, although 
the integration [26.8] runs from minus infinity to plus 
infinity. Of course, it is true that the major contribution 
to the integral comes from the region of the saddle. We 
further note that the first terms of AY and A® in [26.17] 
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and [26.18] are critical, because the error is of the order 
of magnitude of the expression itself. However, these terms 
cancel in [26.21]. 


27. THE WKB METHOD 


G. Wentzel and L. Brillouin have given a procedure with 
which it is possible to obtain the approximate solution of the 
wave equation when it is separable in the individual inde- 
pendent variables.2 H. A. Kramers and his students gave 
the mathematical foundation of the method and extended 
it. The formalism for obtaining the solution is analogous 
to the method in optics which leads from the general wave 
equation to the eikonal equation. 

We solve the one-dimensional and time-independent wave 
equation, 


2 
wt <5 [H—V(«)]Ju=0, where E—V(«)>0, [27.1] 


with the trial solution 


u = exp 


| [27.2] 


In this way, from the second-order differential equation of 


2 L. BRILLOUIN, Compt. Rend, 188, 24 (1926); J. de Physique 7, 353 (1926); 
G, WENTZEL, Z. Physik 38, 518 (1926); also see H. Jerrreys, Proc. London 
Math. Soc. (2) 28, 428 (1923), 

*>H. A. KrRaMERS, Z, Physik 39, 828 (1926). 

‘In optics (see Ref. 1 on p. 127) we solve the time-independent wave 
equation, 

dy 

= + n®key=0, 
with the trial solution 

y = exp[tkS]. 


We assume S to be a slowly varying quantity; namely, only the largest power 
of k) in the expression 


oi = - k2 (=) ik = 
da? “Ol aa * dat ‘4 
is used. From this follows the differential equation for the eikonal which 
is characteristic of ray optics, 
BS ote 
(=) en, 
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the first degree for uw, Eq. [27.1], we obtain one of first 
order and second degree for dS/dz = S’ (Riccati differential 
equation): 

S'* = 2m[E — V(x)] + th8’. [27.3] 


To solve this equation, we expand S in powers of h/i (this 
power series is of the asymptotic type): 


h 
S= 8428, +... [27.4] 
Substituting [27.4] into [27.3] then yields 
Si = V2m[E — V(a)], [27.5] 
which leads to the solution 
So =| Vint V(a)|da. [27.6] 
If we set 
(= So ) 
then, to zeroth order (k=0), we obtain the relation 
p(a) = V2m[E — V(x)], (27.7) 


which is known from classical mechanics. 
The trajectory of the particle is classical in the interval 
(%,, #,) in which 
E—V(#%)>0. 
It is also possible for there to be many such intervals; 
however, for the present we assume only one (Fig. 27.1). 


Ly oe z 
Figure 27.1 
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From [27.3], the next approximation is 


28:8) + So =0, 
S, = —4logS&, + constant . 


For the solution to the wave equation [27.1] we now obtain 


i= a pl foe ar}; [27.8] 


2 


E—V>0: oscillatory (classically attainable), 


E—V<0: damped (classically not attainable). 


The + is to remind us that we are concerned with two 
linearly independent solutions, corresponding to the two 
signs of the square root [27.7]. This solution is not usable 
in the vicinity of 


x) = V2m[E — V(x)]=0; [27.9] 


that is, it is not valid in the neighborhood of the classical 
turning point,’ where, indeed, #=V(#). 

In this region we can make a connection with the Airy 
functions that we have already seen. For this, the ex- 
pansion 


rt) Bi V(x) =(v — @,)F,+..., F,=—V'(a,)> 0, [27.10] 


must be substituted into Eq. [27.1]: 


=! — nee, [27.11] 


The solutions of this equation are Airy functions of the 
variable 
/2mF, 


E= (%—2)- he 


[27.12] 
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In order to find the correct combination of A™ and A®, 
we first consider the region x<,, in which u must be 
damped (classically not attainable). With p=—1|p|, 
[27.8] reads 


constant a) rp 
i X exp |— 7 \p|da|, COB}. [27.13] 


VIPI 
From [27.10] and [27.12] we obtain 


1 zy 5 ¥ Lay 
5 fiplae = |? [jo —a,|sde— 141". 


Then, according to [26.21], we have, for 7<42,, 


uw = constant x |¢|-* exp [— % |&|?] = constant x A(é), 
and the solution for «>a, is given by [26.22] as 


u = constant x é-* cos (3 é# — 42) 


constant aL : 
= ————  XCOs |= ae — 42). [27.14] 
VP a 


Boe 


Similar formulas are valid for the second turning point, 2. 
With 


PUo) _ V(v)= —(#—2@,)F.4+..., 


2m 
/2mF 
F=+V(q)>0, 9=/7@2a,—a), 
we obtain 
constant 1 f 
dite conte x exp|—; [Ip |az| = constant 
V|P| “ 


x|n|-texp[—3ln|#] =constant x A(7) [27.15] 
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in the region #>4,, and we then have 


u = constant x 47? cos (3? — £2) 


constant 1 : 
= ———=— X COs [= dr— 42] [27.16] 
/p fal ‘ 


x 


for «<@,. These expressions represent the wave function 
in the region from —oco to +oo, The wave function van- 
ishes exponentially at these limits. It changes its structure 
basically at the two classical turning points, 7, and 2. 
For a more detailed discussion of the behavior of these 
solutions, reference should be made to the original papers 
by Kramers and his co-workers. In order to obtain a 
unique description with the two pairs of solutions—[27.13], 
[27.14] and [27.15], [27.16]—we must also require 


as 
1 n 
; 5 | peorae ig fe [27.18] 


be 


in the common interval (v,, #). If we denote the famous 
phase integral by 


J = ved = 2(V2m(E — V(a)) de , 


Eo 


we no longer obtain the old Bohr quantum condition, 
J= n° 2ah ; [27.19] 
rather, we obtain the refined condition 
J = (n+ 4)-27h (t==0; 125.) [27.20] 


Of course, because of our derivation, even this expression 
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is an approximation, which becomes better as n increases. 
The error is of order 1/n; however, the constants are cor- 
rect. In some cases, as, for example, with the harmonic 
oscillator (see [15.20]), the quantum condition [27.20] is 
obeyed exactly for arbitrary n; this depends on the par- 
ticular form of the potential.’ 

Until now, we always assumed only one interval in which 
E—V(z)>0; this was an essential assumption for our 
derivation. If there exist more than one such interval, 
then new effects appear because of the fact that the wave 
function does not vanish in the classically unattainable 
intermediate regions (see [27.13] and [27.15]). Therefore, 
it is possible for the wave function to ‘“‘leak’’ from a re- 
gion (%,, 2) to another region (%,, %) in spite of their 
being separated by a classically insurmountable potential 


aos ae ae Dy x 


Figure 27.2 


barrier. This typically wave-mechanical phenomenon is 
called the tunnel effect; it plays a large role in numerous 
applications of quantum theory. We shall concern our- 
selves with the effect in an exercise (see Sec. 41). 


5 See, for example, E.C. Kemble, The Fundamental Principles of Quantum 
Mechanics, pp. 574 ff. 


Chapter 8. Matrices and Operators. 
Perturbation Theory 


28. GENERAL RELATIONSHIP BETWEEN MATRICES AND OPERATORS. 
TRANSFORMATION THEORY 


We recall the discussion of Section 16 and repeat the 
definition of a matrix element given there. Let 


U1(J) 5 --+y Un(Z), «+ 


be a complete, orthonormal set of functions. Then we can 
write 


Pun = Su(k|F ln), - [28.1] 


where J 
(k|F|n) = futteyag [28.2] 
represents the matrix element of the operator F with re- 


spect to the u,. 
According to [28.2], the Hermiticity condition [12.5], 


|eurruaa = |ui(Fu,) dq [28.3] 


means 
(k| F|n) = (n| F|k)* [28.4] 


for the matrices. Defining another matrix analogously, 
uy = Zus(miGk), —(k1G\m)=[uk(Gun)ag, (28.5 


138 
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we calculate the matrix of the operator GF: 


(GE )u = G(Fu) , [28.6] 


ie 2G 
= Zim D (m|@|k)(k|F\n) = =D Uin(n |GF |). [28.7] 


The last step in [28.7] was carried out in the usual sense 
of matrix multiplication, 


(m|GF\n) = > (m|@|k)(k|F\n), [28.8] 


k 


which shows that multiplication of operators corresponds 
to multiplication of the associated matrices. 

We now want to investigate how the matrices change if 
we go from the wu, to another complete, orthogonal set of 
functions v,. We can expand v, in terms of the u,: 


y= > u.(n|8|A) = Suz, [28.9] 


where 


(n|8|A)=|o.undg . [28.10] 


With each function f = a,u,+ a,u,+..., the operator S de- 
fined in [28.9] associates a function g=4a,v,+4,0,+... with 
the same expansion coefficients. We call S a transformation 
operator; it is fundamentally different from the Hermitian 
operators F considered so far. 


a. Unitary transformations 


As a result of the completeness relation [9.9], we have 


fotea=inigiarmisiz). (8.4 


n 


With this relation, the orthogonality and normalization 
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condition for the v, becomes 


> (n|8|A)* (| 8|B) = (A|8'S|B) = ous 
a , [28.12] 
or S'sS=1 


where S' denotes the Hermitian conjugate matrix obtained 


from 8S: 
(A | S*|n) = (n|8|A)*. [28.13] 


A rule for the Hermitian conjugate matrix is 
(FQH..)t= x, GE. [28.14] 


We can require, as a condition for the completeness of 
the v,, that each u, be expandable in terms of the v,: 


tig =D v4(A | S41) , (4) Stim =[usriaa [28.15] 


If the v, are to be complete, then the condition analogous 
to [28.12] must hold: 


(| SS*|m) = dam s SS*= 1; [28.16] 


this means that S"', the operator inverse to S, exists. 
For square matrices of finite rank this is always the case; 
that is, [28.16] always follows from [28.12] and vice versa. 
However, this is not the case with rectangular matrices 
and, therefore, it is also not the case with matrices of 
infinite rank, since then it no longer makes sense to speak 
of a square matrix. 

A matrix which satisfies [28.12] and [28.16] is called 
unitary, and the associated transformation is called a wni- 
tary transformation. A unitary transformation preserves 
orthogonality, normalization, and completeness. 

We now want to investigate how the matrix F defined 
in [28.2] transforms when we go from the u, to the ,. 
Using [28.9], and in analogy with 


(k| Fl) = | ut(Fu,) dq, [28.17] 
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we write 

(A| F"|B) = fetceenyad 

=| > DY (A|S*|n)ut(Pun)(m |S |B) dg 

= > > (A|S*| n)(n| F|m)(m|8|B)=(A|S*FS|B). [28.18] 


Thus, 
F'=S'FS (similarity transformation) [28.19] 


and, since S'= 8 for a unitary transformation, 

dey vi eh [28.20] 
It is easily shown that the Hermiticity of F is preserved 
by this transformation; that is, #’ is Hermitian if F is. 
b. Formulation of the eigenvalue problem 


If we write a matrix element of the Hamiltonian operator 
(n|H|m) = | utiun) ag [28.21] 


with respect to eigenfunctions of the wave equation, 


Hu= Hu, 
we obtain 


(n|H|m) = Bf wit dg = Eaeea; [28.22] 


that is, the Hamiltonian matrix defined in this system is 
diagonal and its diagonal elements are the energy eigen- 
values H#, of the wave equation. Thus, solving the wave 
equation corresponds to diagonalizing the Hamiltonian ma- 
trix (written with respect to some orthonormal set of func- 
tions); that is, we make a principal axis transformation 
with a unitary transformation matrix S which must satisfy 


(n| S*HS|m) = Endnm y [28.23] 
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according to [28.19]. Using SS'=1, we can also write 
(n| HS|m) = (n|SE|m) , [28.24] 


or 
Dd (n| A|k)(k|S|m) = (n|S|m)E,, . [28.25] 


In this fashion it is, in principle, possible not only to cal- 
culate the energy eigenvalues but also to calculate the 
associated eigenfunctions wu, of the wave equation; never- 
theless, the calculation can be carried out in only a few 
cases, for example, with the harmonic oscillator and in 
perturbation theory. 


c. Extension of the matrix method to continuous spectra 


The operations of matrix calculus can be generalized to 
the case in which continuous variables replace the dis- 
crete indices. If, for example, we consider the matrix 
product F-G, 


> (n| FIRM k|E|m) , [28.26] 


and assume,k to be a continuous variable, then we need 
only replace the sum by an integral: 


>for, vony ky) O7K « [28.271 


In place of the Kronecker 6, we have the Dirac function: 
Onn! > O-M(k)O'(k — k') [23.28] 


In this way we can, for example, rewrite conditions [28.12] 
and [28.16]: 


StS=1, Y(A|S*|n)(n|S|B) =07(A)6(A—B); [28.29] 


n 


SSt=1, Jo |S|A)o(A)dA(AS*|m) = dam. [28.30] 
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Since the “density function” @ can also have a discrete 
Spectrum, it is even possible to have a juxtaposition of 
discrete and continuous spectra. 


29. GENERAL FORMALISM OF PERTURBATION THEORY IN THE MA- 
TRIX REPRESENTATION 


We recall the example in Section 16 of a linear oscillator 
with a perturbing potential, and we make a more general 
assertion: We can carry out a perturbation calculation if 
the Hamiltonian matrix has the form 


H= H+ V, [29.1] 


where H™ is diagonal, and where V is not diagonal but is 
small compared with H by the criterion 


[(m|Vin)|<|(m|H|im)—(n|H|n)| (m#n). [29.2] 
Using [28.25] we can write 
E9(n|S|m) +S (m| V | k)(k|S|m) = (n|S|m)By, [29.3] 


where the H® are the eigenvalues of H®. This relation is 
still exact. Now, however, because of [29.2], we make an 
expansion: ? 


E, = EO + EO + EO 4...., [29.4] 
S= 8O+ SO 4 SO sane So — vie [29.5] 


1 Because of [29.2], we can first set (n{F|k) in [29.3] approximately equal 
to zero for k#n: 


{EQ + (n|V |n) — Em)(n |S |m) ~ 0. 
We can number the eigenfunctions of H such that 


En = EO + (n|V |n) 
and 
(n |S |m) = 0 (m #7). 


Thus, because of SSt= S™S =1, we have 
(n |S |n)~ 1. 


This justifies the expansions [29.4] and [29.5]. 
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(We have avoided introducing a parameter «. Instead of 
writing 
S=14 8+ e824 28+..., 


we stipulate that the quantities Z® and S®° in the expan- 
sions be, respectively, one order of magnitude smaller than 
Eo» and Ss) 


a. First approximation 


Substituted into [29.3], this yields 


ED bam + Ex (n|S®|m) + (n| V|m) 
= EO Sam + EO 6am + (n| 8° |m)Eo, 


Ey (n| S| m) + (n| V|m) 
= ED Sum + (n|S|m)E, [29.6] 


to first order. In solving this equation, we want to distin- 
guish between two cases: 
1. nem: 

(Ey — E,,)(n| 8 |m) = — (n| V|m). [29.7] 


In case H® is degenerate, this equation can only be sat- 
isfied by 


4 


(n|V|m)=0 for EYP =ER, nxm; [29.8] 


thus, if there is degeneracy, this condition must be satis- 
fied in order that it be possible to carry out the pertur- 
bation calculation. The approximation will also be bad if 
the eigenvalues are very close together; V must satisfy 
condition [29.2] for this reason. From [29.7] we obtain 


(n|V|m 
(n|S@|m) = — oe a [29.9] 
2.m=n 
E® = (n|V|n). [29.10] 


Thus, in first order, the energy eigenvalues are displaced 
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by amounts equal to the diagonal matrix elements of the 
perturbation. 
Since the transformation matrix S must be unitary, 


S'S = 1, 
we have 
sm = St 0, [29.11] 


in first order. This means that (n|S|n) is purely imaginary 
but otherwise completely arbitrary, corresponding to the 
fact that we can always make a phase transformation. 


b. Second approximation 
In second order, [29.3], [29.4], and [29.5] yield 
Eb m+ B2(n|8®|m) + E(n|$®|m) + (n| V|m) 
+ > (n| V]k)(k| SO | m) = BOS pm + EL + EES nm 
ae i |S] m) Ee + (n| S| m) EB, + (n|SO|m)ED ; 
with [29.6], this becomes 
(BS — BS)(n| S|m) + ¥ (m| Vk) (Fe| S|) 
= Oamly + (n|S®|m) EY. [29.12] 


Again, we distinguish between two cases: 
1. nA~m: 
n| V\n) — (m| V|\m)} (n| Vm 
ea? 
(| Vk) (| V|m) 


(n| $|m) =" 


i 2 (Ey — EE? = B®) [29.13] 
ken oS m 
kAm 
Pre na=—=—™: 
(mi VIR ELY |») | (m| Vk) 2 
E® = > So a a y Eo E®© Eo [29.14] 


k 
kfn ia 


(always negative for the lowest eigenvalue). 
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c. Degeneracy of H® 


In case H® is degenerate, we have already seen that V 
must satisfy the condition [29.8] before we can perform a 
perturbation calculation at all. In order to achieve this 
with g-fold degeneracy (EO =H =...=E® =E,), we must 
pick out the corresponding g-dimensional subspace and solve 
equation [28.25] exactly in this subspace, 


$ (n|H|k)(k|S|A) = 2B, (n|S|A), [29.15] 


in which case [29.8] is satisfied. Equation [29.15] represents 
a set of g systems of equations (labeled with the index <A), 
each system containing g linear, homogeneous equations 
(labeled with the index n) in the g unknowns (n|S|A). 
The matrix of the coefficients of each of these systems of 
equations is 


H—£,:1 
Ey+ Vii— Ey Vis° cane i cp Vi, 
Vou E,+ Ve,— Ey 
= “4 a ; » [29.16] 
UG pie ee oe a E,t+ Vi =e) 


and, as is well known, there are nontrivial solutions of the 
system only under the condition 


det || H—H,-1|=0. [29.17] 


This is an algebraic equation of the gth degree for E,, 
whose solution determines the energy eigenvalues. Once 
this equation has been solved, there is no further obstacle 
to solving [29.15]. 

In conclusion, we remark that this eigenvalue problem 
is nothing more than the famous principal axis transfor- 
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mation of a gXg matrix H, the Hermiticity of which re- 
quires that all of the eigenvalues be real. Equation [29.17] 
is known as the secular equation because it first arose in 
connection with the calculation of secular perturbations of 


planet trajectories. 


30. TIME-DEPENDENT PERTURBATION 


If we have a time-dependent perturbation, then of course 
we are only interested in a solution of the time-dependent 
wave equation 

-, OW 
ih-* = Hy + Vithy. [30.1] 


We expand the solution of this equation in the usual way 
(see, for example, [7.7]), in this case with respect to the 
unperturbed eigenfunctions uw, (as proposed by Dirac): 


y= > a,(t)u, exp |- ; Bo : [30.2] 


where 
Ayu, => EU, . [30.3] 


Substituting [30.2] into [30.1], and remembering [30.3], 
we obtain 


ih > = Um @XP |- ~ Be J 


= ; Ba [30.4] 


= D4, (t) Vm exp 


We multiply this equation by u, and integrate over dq: 


7 Ady aa 
ih ia exp - i E, ] 


= Dd a,,(t)(n|V|m) exp - ad] [30.5] 
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or 

da, a (0) (0) 

ih = = > a,,(t)(n | V|m) exp + 7(En SE 


We seek a solution for given initial values, 
a,(0) =a”. 
We also expand this solution: 
a,(t) = ay + an (t) + an (t) +... 
Then the initial conditions are 
a)(0) = a?(0) =... = 0. 


For the sake of brevity we introduce 


(n|Q|m) = (n|V|m) exp 


;(H?— B2)4), 
and, by integrating [30.6] over ¢, we obtain 


t 
a(t) = — => a] (n|Q|m)at, 
h yg 
0 


a2) =— FS m1 QinaP eat 


=~ FZ ak S/n Olniayar fe] @imye'yar’ 
m t 


| Chap. 8 


[30.6] 


[30.7] 


[30.8] 


[30.9] 


[30.10] 


[30.11] 


These formulas are valid for time-dependent V,? as well 
as for time-independent V.* However, if V is constant in 


* For example, they can be used to caleulate the induced emission and 
absorption of light, in which case the perturbing operator F is given by the 


incident radiation field. 


’ Since we only integrate from 0, we can also assume that a constant per- 
turbation V is “‘turned on” at time {=0, for example, by the transition of 


a@ system to an excited state. 
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time, we can evaluate the integrals: 


exp [(i/h)(E,” —E))t] —1 


an) = — Fala] Vm) re ao , [30.12] 
a®(t) = Ya® S (n|V [0 (| V |m) 
"exp l(é/a)( 22 — Hy4]—1 
(B® — H)\(H— BS) 
_ exp G/AH2— BPI] gy ay 


(HO — BYE — Be) 


The approximation is not necessarily bad if H° =E® (res- 
onance denominator) because, as a result of having inte- 
grated from 0 to t, the numerator vanishes then also. 
Thus, the term in [30.12] containing the resonance denom- 
inator is 


a 
—am(n| V |) , 


and the approximation is usable for times such that 
l(n|Vim)tl<h. 


On the other hand, a solution valid for a longer time can 
only be obtained by treating the degenerate subspace 
exactly. 

We now treat the important special case in which one 
of two energy values, H®, lies in the continuous spectrum. 
This is the case, for example, if an atom makes a transition 
from an excited state to the ground state and, instead of 
emitting a y-quantum in the process, ejects an electron 
from an outer shell;4 the free electron has a continuous 
spectrum (or, in the case of distant walls, very densely 
distributed energy eigenvalues—a quasi continuum). The 


4 This radiationless transition ig known as the Auger effect if the excitation 
of the atom is a result of removing an inner electron. Details can be found, 
for example, in E. H. S. BurHop, The Auger Effect and Other Radiationless 
Transitions, 
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initial values that we use in [30.12] are 
a—=—O0fornA~Am, a, =1 for n=m; [30.14] 


that is, we consider the ejection of an electron which is 
known to come from state m. We use k as the variable 
in the continuous spectrum, and we obtain 


exp {(i/h)[ HO (k) — Bm lt} — 1 


aW(k, t) = — (k|V|m) E(k) — B® 


[30.15] 


We now calculate the transition probability for a transi- 
tion from the state m to an arbitrary state k within the 
time interval from 0 to ¢t, where state & lies in the interval 


EO —4AR< B(k)< HP + AE. [30.16] 
The result is 
W(k, t) =| aD(k, t) |2dk 


4 sin? (t/2h)[ E(k) — BO} 
[HO(k) — BY} 


~ [Chi V|my|e- fa [30.17] 


Let P(Z) be the number of states on the energy shell 


” 


dk = P(E)dE. [30.18] 
With the abbreviation 
2 - [B(k) — B®}, [30.19] 
we then obtain 
+00 
W(k, t) = P(B)| (| V]m)|2 4 x a dx. [30.20] 


We justify the limits of integration —co to +00 by the 
re quirement 


AE-t 
= [30.21] 
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The integral equals z; therefore, 


Qa 
W(k, t) = P(E)|(k|V|m) . [30.22] 
This formula is of especial importance, because it gives 
the transition probability between two completely arbitrary 
states. Because of the importance of the formula, Fermi 


called it the “Golden Rule.” 


Chapter 9. Angular Momentum and Spin 


31. GENERAL COMMUTATION RELATIONS 


The eigenvalue problem for angular momentum and the 
corresponding transformations have an important place in 
wave mechanics. We have been partially concerned with 
them in an exercise treating the eigenvalue problem of the 
symmetrical top (see Sec. 47). We restrict ourselves here 
to the most essential points. Consult the literature for a 
detailed treatment of these problems (for example, P. A. 
M. Dirac, Quantum Mechanics). 

For a particle with angular momentum 


1 (measured 
P=(P,, P., P;)== 
Pry Eer £0) 7) 52 enna 
il 1 P) P) 31.1 
P= ; (LPs an Pz) =: (2 de, == ihn =| ee : [ ] 


(cyclical permutations of the indices) 
we have the commutation relations 
[P., PoJ= P,P,— P,P, =iP,,... [31.2] 
(cyclical permutations of the indices). 


These commutation relations also hold for an arbitrary 
number of particles, where 


P) ra) 
==> [af =a, 
BFE (Pap Sapp 


> (aps? — ws’py), .... [31.3] 


152 
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The angular momentum operator is introduced purely 
formally into wave mechanics for describing the transfor- 
mation properties of the wave function under rotation of 
the coordinate system. On the basis of the kinematical 
properties of the rotation group, the following more gen- 
eral commutation relations can be given which are inde- 
pendent of the particular form [31.3] of P: 


lee; C} =0 (k = 1, 2, 3), [31.4] 


[P,, A,] = 0, eee 


[Py Al = —[P,, AltA, ... }, [31.5] 


where C is a scalar function of p® and x”, and A, is the 
kth component of a vector A( Bp” x”), As examples, we can 
take C=H (Hamiltonian), C=|P|?, A=x”. The eigen- 
values and matrix elements = the general operator P can 
be derived from these relations by purely algebraic cal- 
culation. 

Of course, we can also proceed analytically by expressing 
the differential operators [31.1] in polar coordinates (see 
the exercise on this topic, Sec. 45). We find that the 
eigenvalue equation 

Pea A Y [31.6] 


is identical to Eq. [18.9] for the spherical functions Y,(0, 9). 
Thus, P? has the eigenvalues 
A=IU(1+1) (Werle ce Oe [31.7] 


If we further choose the 2/+1 linearly independent Y, so 
that they are simultaneously eigenfunctions of P,, which 
is possible because, according to [31.4], P? commutes with 
the P., 


o 
ee Ye, 


1 See, for example, M. Born and P. JornpAn, Elementare Quantenmechanik ; 
P. A. M. Dirac, Quantum Mechanics. 
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then we arrive at the spherical harmonics Y, ,, (see [18.18]), 
with eigenvalues 


pm; m=—l,—1+1,...,+1. [31.8] 


32. MATRIX ELEMENTS OF THE ANGULAR MOMENTUM 


Here, we give a summary of the matrix elements of the 
angular momentum operator. As we have already empha- 
sized, these can be derived purely algebraically from the 
commutation relations [31.5]. 

We now denote the quantum number / by 7. By so doing, 
we mean to leave open the possibility that 7 assumes half- 
integral values, in contrast to the particular form [31.1], 
which, according to [31.6] and [31.7], allows only integral 
values of 1. 

Here, we write the matrix elements for fixed 7 in the 
form 

(j, m'| |j, m), 
where m labels the matrix elements subject to the con- 
dition —j<m<-+j. The matrix then has the following 
appearance: ? 


m=j m=j—1 oe. M=—j 
m=j—11(3,3—-11 13,97) 3, 7-1) 14,9-0) 
m=j—2 Oi 2 7) 


m=—j |(j,—J| 19,9) 
* This notation means that, for example, there are the folowing matrices 
for j=: 
1 
0 


’ = 


0 
Pit r= | 
0 


0 0 
’ P,—iPy= 
sa) 


Corresponding to these, there are the matrices 


1 
|e il 
21 0 


1 


ie 
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In the representation in which P? and P, are diagonal 
(as is well known, commuting operators can be simulta- 
neously diagonalized), we find 


(3, m'| P?|j, m) =JJ+1)bun's [32.1] 
(J; m' | P; lj, m) =a 0 mm ’ : [32.2] 
and, for the non-Hermitian matrices 


Pitt, , 
we further find 


Q, m+1/P,+ iP, (7, m) = j, m| BP, — iP, |j, m+ 1) 
=V(j—m)(j+1+m). [32.3] 


All other matrix elements vanish. 

For every vector A—in particular, for the coordinate 
matrices—we can derive the following general expressions 
from the commutation relations [31.5]: 


(j+1, m+1|A,+ 7A, |j, m) 

=FG+1 |Al)V GE m+ 2G+ m+) 
(j+1, m|A;|j, m) 

= (j+1/Alf)V G+ m+1)\(j—m+ 1) 

(j, m+1|A,+7A, 


j, m) aaa 
= (j/AlLNV Gj Fm)\jtm+1) 

(j—1, m41|A, + 1A, |), m) 
=+(j—1|A|j)V (GF m)j-m—1) 
(j—1,m|A,|j, m) = (7 -—1| A] VG + mj — m) 


[32.4] 


For all other pairs of values j, m in the initial and final 
states the matrix elements vanish. The expressions (j’|A|j’) 
are numbers which are independent of m. 
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33. SPIN 


In order to explain experimental observations—first in 
connection with anomalous Zeeman splitting—a spin was 
ascribed to the free electron [A-3].° This fact can be included 
in the theory by means of the general representation [32.1], 
[32.2], and [32.3] of the angular momentum. This repre- 
sentation indeed does not rest on the special definition 
[31.3], which cannot be applied to spin; instead, it is based 
solely on the relations [31.5]. On this basis, it is possible 
to introduce spin into the theory in a purely formal way.* 
In the relativistic treatment of the electron (Dirac equa- 
tion), it turns out that the spin formalism discussed in 
this section is already contained in the equations and be- 
comes evident in the limiting case of small velocities.® 

We now want to consider the generalized description of 
a particle with spin. By the spin s of a particle we shall 
always mean an angular momentum whose magnitude (as 
opposed to its components) is always a fixed number. We 
denote the spin operators by s,. Then, for example, the 
commutation relations analogous to [31.5] hold: 


, $182 = 8,8) at 183, eee [33.1] 


(cyclical permutations of the indices) . 


Since we consider |s|? as a fixed number, which, according 
to [32.1], must be of the form 


|s|*= si+ 83+ 83 =s(s+1), 


with s half-integral or integral, we can introduce one of 
the components s,, for example, s,, into the wave function 
as a new independent variable: y=y(q,s;,¢). However, 


°*G. EK. UHLENBECK and S, GoupsMiT, Naturwiss. 13, 953 (1925); Nature 
117, 264 (1926). 

“W, PAULI, Z. Physik 48, 601 (1927). 

*P. A. M, Dirac, Proc, Roy. Soc. (London) A 117, 610 (1928); A 118, 351 
(1928). 
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since s, can only assume the values —s, ..., +s (see [32.2]), 
we can also write 


wd, §,;t)= = C183) pu; t) 3 <w<s), [83.2] 
7 


where, for example, the C,,(s,) are defined by 


1 for s5,=4p 
C (s => ; 33.3 
we | 0 otherwise 
and satisfy the orthogonality relation 
i fora — 
¥ 0% 83)* Cy(83) = ee [33.4] 
8,=—8 0 for be se re 


The effect of the operators s, on a wave function is most 
easily seen in the matrix notation: 


82° Yulds 1) = > pul’ | lH) » [33.5] 
B 
with the matrix elements given by (see [32.3]) 


(u+1/8,+ i8,|u) = V(s Fu)(st+1+p) 


[33.6] 
(“| 83 |“) = 


Thus, for example, we have 


$3°Yu= Yu's 
81° Yu= 3 (81+ 182) ppt 3 (81 — 182) Pp 
= tyaV (s+ e)(8 + 1-2) 
+ dyn V(s—w)(s+14+y), 


where the first term vanishes for ~——s and the second 


term vanishes for w=s. 
The most important and most fundamental application 
of this formalism is to the spin of the electron. Since 


s=%}, which means pw=—}34,+%4, [83.7] 
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everything is particularly simple: 


|s|?= si+ 83+ 83 = 8(s+1)=}, 


0 1 _ fo 0 4 0 
a+ in|) ol a—in=|, ’ Soe|| «yas 


a; 
2 
With a new operator defined by 


o=2s 5 ' [33.8] 
we obtain ° 


[33.9] 


These operators obey the commutation relations 


[O15 C,] = 10, — 0,0, — Zia, «.. [saat] 


(cyclical permutations of the indices) . 
In addition, they satisfy the relations 


5 GOP — C10 — Vu: 
(cyclical permutations of the indices) 
or 
0,02-+ 6,0, = [0,, 02], =0,... [33.12] 
(cyclical permutations of the indices) . 


Because of the latter relations, the o, are said to anti- 
commute. 

Now, in addition to the position coordinates and the time, 
the wave function of an electron will also contain the dis- 
crete spin variable s,; the spin variable represents the 
additional degree of freedom of a spinning electron. In 
correspondence with the two eigenvalues s,=+4 and —3, 


*The matrices o,, o,, and o; are called the Pauli spin matrices [note by 
the students]. 
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we divide wy into two terms, y, and y»,, as in [33.2], and 


write 


(9; 833 t) = (33.13) 


Wi(q; t) | 
w2(q; t) 


The normalization condition is 


[vvear=fiv.tav+fiyteay = 10 


where |y,|/?dV and |y,|?dV can be thought of as the prob- 
abilities that the electron has its spin respectively parallel 
or antiparallel to the positive z direction. 

Obviously, this formalism can be immediately applied to 
the usual calculational methods of wave mechanics (trans- 
formation theory, perturbation theory, etc.). In general, in 
addition to p and q, the Hamiltonian function will also 
contain the s,; for example, if there is a magnetic field H, 
the Hamiltonian will contain an additional term equal to 
a constant times 


8, H, se 8, H,+ 8H, = 


The generalization of this formalism to NW particles with 
spin is immediate. The wave function will now also contain 
the spin variable s;” of each of the W particles (a =1,..., V), 
where s has one of the values —s®,..., +s. For an 
electron we have simply s\?=+4. A total angular mo- 
mentum, 


J= S f(x x p®) =f hao}, 
a=1 — 


can be defined which obeys the usual commutation rela- 
tions [31.5]. 


34. SPINORS AND SPACE ROTATIONS 


Here we want to investigate how the spin formalism that 
has been presented behaves with respect to rotations in 
space. 
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For this purpose we consider a 22 unitary matrix S, 
for which we specifically require 
det S=1, 

(From unitarity, SS'= 1, it only follows that |det S|?=1; 


that is, there is a trivial phase factor undetermined in 
det S=e, where «a is real.) 


Say Sis 
R S ) det S= Sas Sas — Sas Sa 1 . [34.1] 
Bay Ses ate 
From the condition 
Sa— st, with S- | ee [34.2] 
— — 1 a , 
= aiid ae — 8 Si : 
we then obtain 
Se oa So. = — aa | Si. |?+[S,2[*=1. 
With these results we can write 
= Sy Si. 34 3 
=~|- sh stl) 343] 


We shall show that the transformation represented by the 
matrix S is a rotation. This representation of a rotation 
in space was already known before quantum mechanics.’ 
The matrix S can be written, for example, as a function 


Figure 34.1 


7 See F. KLEIN and A. SOMMERFELD, Uber die Theorie des Kreisels. 
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of the three Euler angles (Fig. 34.1): 
Si.= cos Sexp 5 (y +2] 
Sig = 7 8in cope Re — »| am 
2 2 ‘ 
We now introduce a two-component mathematical form, 


eee é,||, [34.5] 


called a spinor by Ehrenfest, which is transformed by 8 
as follows: 


&=E-S8 &=§ Sit ae wie eeu [34.6] 


= & = €,81.+ E, Soo = &,81.+ oe Si 


We see that our two-component wave function [33.13] 


Yi 
Wo 


y= [34.7] 


transforms contragrediently to the & with respect to space 
rotations: ® 


yi= 8 -y, yp S291 — Spy.= Shy, — SrePe2 | . [34.8] 


1 —- 
Wa Sa Sor Y1 ae SuYe= Shy ap Si Ye 


We note that |é7, || transforms like || &, —é,||, and that 


that is, the formulas remain correct if we make the fol- 
lowing substitutions: 


= 
in transforms like } 


We 


E 
—Y 


bi é, vir . [34.9] 


&>—& Ys) 


8 The contragredient nature of » with respect to & is characterized by the 
invariance of the form (&¥). 
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At this point we make the important observation that 
every unitary transformation, including our S, leaves the 
form 

W = (€6*) = 4,814 &.83 [34.10] 
or 
0 = (pty) = pip t yay [34.11] 


invariant. This can be demonstrated immediately by sub- 
stituting &’=&S and &'*=— 8'é*. Conversely, it can be said 
that this requirement, which is basic to the physical mean- 
ing of the theory, together with the requirements of uni- 
modularity (det S=1) and linearity, determines the form 
of our particular transformation matrix S. 

Now, we want to show that the transformation S actually 
represents a rotation. In connection with o=(a, 02, 63) 
(see [33.10]), we introduce two vectors x and d (d is called 
the spin density): 


w= (fc6*), d = (y*ay). [34.12] 
This notation means, for example, 


=> De.caphs- 


a=1,2 B=1,2 
Using 
: 0 <i 0 0 
a tin= 2) ol)? i= 2) ol? 
we obtain 


0, + it, = 2,85 d,+ id, = 2yi yp. 
0-40, — Bere d, — id, = 2y3y, }. [34.13] 
L3 = bf — &,69 = YrY:—Y22 


We now ask how x and d transform when we transform & 
and y with 8. To answer this, we calculate the expression 


+43 + 2 = (2,+ i@,)(@, — tm.) + a3 . 
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Using [34.13] and [34.10], we immediately obtain 

x + a3 + 23 = N2 = constant. [34.14] 
Similarly, we find 

d; + a; + d; = o?= constant. [34,15] 


Thus, since S leaves the forms |x|? and |d|? invariant, we 
have shown that it induces a rotation when acting on & 
and y. In both cases we have exactly the same rotation. 
As a result of the linear transformation property of such 
a rotation, we can also write 


w= > H%An, ad; = >a, An. [34.16] 
k k 


Using [34.6], [34.8], and [34.12], we immediately obtain 
SoS => Y eA: . [34.17] 
k 


By virtue of [34.17], an S transformation, which, accord- 
ing to [34.4], is characterized by three real parameters, is 
associated with a rotation A, 


SA. [34.18] 
It is apparent from the form of [34.17] that the inverse 
association is double-valued: 

AS, 

A->—S8. 
If, for example, é, and & both change sign (that is, S= —1), 
then x and d are unchanged: 


ee — y=—y, d=d. 
If 
S,>A, Su B, 
then, in the sense of matrix multiplication, 


8,8, AB. 
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We show this with the help of [34.17]: 
8, OG; a > Or. Ans» 
k 
Su 9; Si = = 2 By ’ 


S(Sucite s,s = ze 6.87 Bis = > o( AB), - 
U 


In addition, we recall the rule 
(S Sey ‘aa ihe st ) (S; 8,,)' 8, Sy =, 


The content of this section is very closely connected 
with group theory, which we have not been able to intro- 
duce here. There are detailed textbooks on the relation- 
ship of wave mechanics to group theory ;° for a deeper study 
of the subject, reference should be made to these books. 

Remark: By dropping the restriction that S be unitary 
(but keeping the requirement det S=1), the above consid- 
erations can be extended to the Lorentz group; ! we ar- 
rive at relativistic spin theory in this way. 

4 

®* B. L. VAN DER WAERDEN, Die Gruppentheoretische Methode in der Quan- 
tenmechanik; H. WEYL, Group T'heory and Quantum Mechanics; E. P. WIGNER, 
Group Theory and its Application to the Quantum Mechanics of Atomic Spectra; 
W. PAULI, Continuous Groups in Quantum Mechanics, CERN, 1956. 


1° This is the group of Lorentz transformations, which are well known to 
be characterized by invariance of the form 


eae pie ess 7 


Chapter 10. Identical Particles with Spin 


35. SYMMETRY CLASSES 


We first consider two identical particles labeled by su- 
perscripts (1) and (2). Because the particles are identical, 
the associated Hamiltonian operator must be symmetric 
(that is, invariant under interchange of the particles). 
For example, for two electrons in the Coulomb field of a 
nucleus and in an additional magnetic field, we have 


7 p™ p”” Zer Ze 


2m) 2m rh pe 
e? eh 
a m + uy H(6 + 0); i ee P35. 1 | 
Let 
yr= p(x, 8); x, 83°) [35.2] 


be a solution of the wave equation for this Hamiltonian. 
The additional solution 


yn = p(x, 3. x2, sP) = Pry [35.3] 


follows immediately from the symmetry of the Hamiltonian 
operator. Here P,, is the operator that interchanges the 
two particles (the exchange operator). The linear combina- 


tions 
PoP | (symmetric solution) , [35.4] 
Pi2Ys= Ys 
‘siete. a ; (antisymmetric solution) [35.5] 
PizYa= SS 
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have the important property that the matrix elements 
between symmetric and antisymmetric solutions are always 
exactly zero: 


(a|H|s) -{f,2 P: a pe Ea aadar . [35.6] 


If we interchange the particles in this equation, then the 
left-hand side does not change (symmetry of H!), whereas 
the right-hand side changes sign; thus, we must have 


(a|H\s)=0. [35.7] 


Therefore, a symmetric solution can never develop into an 
antisymmetric one, and vice versa, which is to say that 
there are two separate classes of particles between which 
no transitions are possible for any interaction: ? 


symmetric : : 
t 1 
paitien bosons, integral spin 
[35.8] 
antisymmetric . j : 
half- 
pastels fermions, half-integral spin 


If there are more than two particles (N>2), then, as 
can be shown with the group-theoretical method of rep- 
resentations, there is a class of solutions symmetric in all 
of the particles and a class antisymmetric in all of the 
particles: 

Py,(x x”, CUE ee sy) = Ws (x™, Bae Sua Peal age?) , [35.9] 
PyA(x™, ae. a8 aa”, 85) =vegpaoc Wiss”; cage ie), [Sori 


where P is an arbitrary permutation, and 


ép=+1 for P even, 


ép=—1 for P odd (for example, for one interchange of 
two particles). 


‘The names fermion and boson refer to the fact that particles with half- 
integral spin obey Fermi statistics, and particles with integral spin obey Bose 
statistics, See W. PauLi, Phys. Rev. 58, 716 (1940). 
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There are yet other symmetry classes, which are distin- 
guished from those already mentioned in that they cannot 
be extended from WV to N+1. Of course, if an additional 
particle collides with those already present (for example, 
in an atom), then there are transitions possible between 
these classes. Therefore, there can be present either only 
a mixture of all classes, or only the symmetrical and/or 
antisymmetrical classes. Experiment shows that nature 
oddly makes no use of the first possibility: until now only 
bosons and fermions have been found. 


36. THE EXCLUSION PRINCIPLE 


We now want to investigate the properties of the bosons 
and fermions somewhat more precisely. For this, we con- 
sider particles which are isolated to first approximation; 
that is, we neglect the forces of interaction between the 
particles. 

Let us begin again with two uncoupled particles in states 
m, and n, and write their associated eigenfunctions as 

Un (x, oe) and CPR eae of) ; 
We know that the wave function for uncoupled particles 
is the product of the individual wave functions: 
a 3) “tin, (x, s?) 
sm Un (es a) stile as 83”)}: Cans) [36.1] 


U, = {Un (x 


(@ 9M) 


2 
Ue = {n(x ; sf (2) -) 


) Un lx", 8s 
= Un,(¥, 88°) *Un,(H, 83)}*Cnyn,y [36.2] 
V2 tor 7, 


a [36.3] 
1fjf/2 for n,~n,. 


NyNs 


If the two functions are identical, we obtain from this 


Up = Un (H, 85°) Un (¥s 83) 5 [36.4] 


U,=O. [36.5] 
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For fermions, the state for which both particles are in the 
same state does not exist! That is exactly the prediction 
of the exclusion principle, which was postulated for elec- 
trons even before the invention of wave mechanics.2 We 
now know that it also holds for all other fermions (for ex- 
ample, protons and neutrons). 

From a superficial consideration of the exclusion prin- 
ciple, it might be thought that a sort of action-at-a-distance 
is being postulated, as a result of which even two widely 
separated particles are aware of one another (‘‘sign a con- 
tract”). However, this is not so, because the exclusion 
principle is only valid as long as the wave packets of the 
two particles overlap (Fig. 36.1). If the wave packets do 


Separated wave packets. Overlapping wave packets. 

Particles can be followed Particles cannot be 

individually. followed individually. 

Exclusion principle inoperative. Exclusion principle operative. 
Figure 36.1 


not overlap, then, everywhere in space, we have 


83 )=0, [36.6] 


Un (x, ie) Un, (x 


or, for nonstationary states, 
US, 8s tee x te ne) — 0 [36.7] 


In that case, since according to [36.6] we can have u, =U, 
only if uw, =0, the above symmetrization is of no conse- 
quence. 


*W. Pauul, Z. Physik 31, 765 (1925). 
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The generalization to many particles must simply be 
calculated combinatorially. When we symmetrize we must 
sum over all N! permutations: 


Us = oD Pita”, s$”)} (t,j7=1, 2,...,N), [36.8] 


ti O° > ep Piau, (2, a2)} 
E 


a) 
Un, (% 1 ; gi Un(x”, o.) ‘Ste Un, (x™, aS”) 
(x, sf) 50 ¢ ead : 
=c: : . . [36.9] 
qd) (d 
Un (% , 83 ’) En aoe Un, (x, ag”) 


If any two of the functions uw, are equal, this determinant 
is identically zero. Again, this is exactly the content of 
the exclusion principle for an arbitrary number of fermions 
whose wave functions overlap. 


37. THE HELIUM ATOM 


We obtain the helium spectrum to good approximation 
with the assumption that the Hamiltonian operator is the 
sum of two terms, one symmetric in the space coordinates 
and one symmetric in the spin coordinates. This case, in 
which the spin-orbit interaction is small compared to the 
Coulomb interaction, is also known as Russel-Saunders 
coupling. 

Let us first consider the spins of the two electrons alone. 
For the spin eigenfunctions we take (see [33.3}) 


d for 3,=+ 3 
C+(83) 0 for s,=—+4 37.1] 
C = 0 for s,=+43 _ 
=) ealetafor «= — i 


From this we can form only one antisymmetric combination: 


OCR = (C+(s)+ C_(8) — C89) C4 (8). [87.2] 
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Using [33.5] and [33.8] it is easy to verify that 
(on + oF )- C's a2 )=0. ages) 


Thus, we obtain only one eigenvalue (that is, zero); the 
antisymmetric state is a singlet. 

On the other hand, we have three possibilities for forming 
a symmetric combination: 


Oe, &) = oe -O4(99) 


Co(ss', 8) = => a1 C4(89”) C-(88") 


[37.4] 
+ O_(s?)-C+(s3")] 
Chile ee) = OGD): 0a) 
The relations 
bio + o) Cf =1-CE, m,=—1 
bio + o) 08 =0, m,=0 [37.5] 


4(o3’+ 03) C1,=(—1)-C11,  m=—1 


hold for these combinations; that is, the symmetric state 
is a triplet. We have characterized the three triplet states 
by a quantum number m, for the s; component of the 
total spin. We also have 


{F(o + o)}? 08, = 1{3 + 6-0} 02,=205,. [37.6] 


Now, we must correctly combine the symmetry classes of 
the spin eigenfunctions with the symmetry classes of the 
space functions 


us ue) -0(x) + u(x) -0(x)} ; [37a 
ut=— Se (cit) . aa) ae ule) -0(x™)} 3 
/? a i 7.8] 


The total wave function for the electrons (fermions!) must 
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be antisymmetric (with respect to simultaneous inter- 
change of space and spin). There are two possibilities for 
doing this: 


W(x. 9. ax”, 8 Bios ny x) « Ca(g) 9 
( 523) 53 3 ) as We (83 » §3 | [37.9] 


ree 


VetR iss sia ee, nagar”, a) Cyagay wee) 
triplet ;- [37.10] 
=41,0,—1. 


Therefore, the states of the helium atom fall into two 
classes : 


singlet: para helium, 
triplet: ortho helium. 


In the ground state, 
“=v, 80 that uw°=0, or V°=0, [37.11] 


which means that there is only a singlet term. 

In the approximation considered here, in which the Ha- 
miltonian operator is symmetric in the space coordinates 
by themselves, it follows from [35.7] that ortho and para 
states do not combine. This is a result typical of wave 
mechanics that was not understandable earlier. 

In order to determine the energy eigenvalues, we con- 
sider the Coulomb interaction between the electrons, 


: 2 
V(x, x) = V(x, x?) = = ; [37.12] 


12 


as a perturbation (an exact solution has not yet been 
found). With the Coulomb integral 


Ae = fiw) [2 | o(x) Ve, x2 _) d3a™) d3yl2) [373] 
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and the exchange integral * 


dh a Ve) (x?) O(a) oe) . Vis; xi?) dea™ dia , 
[37.14] 


we find the following shifts in the energy eigenvalues 
using [29.10]: 


AEF ying = Jo + di ; [37.15] 
hire J, =e [37.16] 


The difference 2J, between the triplet and singlet term is 
large compared with the spin-orbit coupling (splitting of 
the triplet state); it is of the order of magnitude of the 
electrostatic energy. 

Better approximations can already be found in the ear- 
liest articles about the helium spectrum.‘ 


> The frequency of electron exchange is connected with the exchange in- 
tegral as follows: 

If we consider the state 

u(x, x, § = 0) = a (ute, 2) + we), 2) , 
2 
at time f we have 
u(x), x?) i= aby fos (x@, x) exp - a E*-t a u(x), x{)) exp |— io E®-t ' 
/2 | h h 

With 


ee 075 
we obtain 


1 2 i a 
u(x), (2) 4) — — {uta x')) ex a74 ux), x@)b exp |—— #2. 
vi p ree + Ue, ue p h t\» 


from which follows 


h 1 
: (=, og =| = —— |u? — u? |= Jute), 2 t = 0) I. 


2d Jt 


The electrons have exchanged their positions in time ¢ = 2h/2J,; the corre- 
sponding angular frequency is 


“W. HEISENBERG, Z. Physik 38, 411 (1926) and 39, 499 (1927). 
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38. COLLISION OF TWO IDENTICAL PARTICLES: MOTT’S THEORY ® 


We treat the collision of two identical fermions (charge e, 
spin 3). Let x be the relative separation of the two par- 
ticles: 

i ee) [38.1] 


Earlier we calculated the wave function for seattering by 
the Coulomb potential (see, for example, [19.8] and the 
remark on p. 106). Here, we can therefore write 


u(x) = P+ S-f(6), [38.2] 
P= exp([tkz + ty logk(r—2z)], [38.3] 
= “exp [tkr — ty log kr], [38.4] 
1 e? 
la ka, =r, [38.5] 
(0) = : 


~ Imv? sin? $0 
-exp[— iy log (1— cos 8) — 2i0(0, —1/y)] [38.6] 


(m=reduced mass). Because the particles are identical, 
we cannot distinguish between 6 and x—6 or between x 


Figure 38.1 


and —x. Classically we would add the intensities at 6 and 
2—6; however, wave mechanically we must write, for un- 


5N. F. Mort, Proc. Roy. Soc, (London) A 126, 259 (1930). 
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polarized particles, 


W(x) = £{3|u(x) — u(— x) |? +1 |u(x) + u(—-)|#}. [38.7] 
triplet : singlet 

The weighting factor of 3 arises because of summation 

over the unobserved spin orientations, m,= —1, 0, +1, 

in the final state. Equation [38.7] contains an interference 
term which is typical of wave mechanics: 


W(x) = | u(x) [2+ |u(— x) |? — d{u(x)u*(— x) + w*(x)u(—x)}. 


The particles with antiparallel spin orientations are still 
distinguishable after the collision [A-4]. 
From the singlet and triplet scattering amplitudes 


[#(8) & f(a — 8) |? = [7(8) |? + [f(x — 6) |? 
+ (f(8)f*(a — 8) + f*(8)f(~—6)), [38.8] 


interference term 


we obtain 
dQ = 2{3 | f(8) — fla — 8) |? + [f(0) +(x — 6) |?}dQ 
ef l, " ul a 
~ 4m2v4 |sint36 © cost30 sin? $6 cos? 40 [38.9] 
’ x Cos (> log os — A dQ 


for the differential cross section. The interference term in 
this formula is characteristic of wave mechanics. The inter- 
ference maxima come closer together as the velocity is 
decreased. The special case 6 =2/2 is peculiar, in that all 
of the four terms in formula [38.8] are then equal. 

For two identical bosons with spin zero (for example, 
« particles), [38.7] and [38.9] are to be replaced by 

W(x) = | w(x) + u(— x) |? 

and 


dQ = [f(6) + f(a — 0) |?dQ. 


According to [38.8], there is an interference term here also. 
The interference term was experimentally demonstrated 
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with, for example, electrons, protons, and « particles; thus, 
the fermion or boson character of these particles was 
verified. 


39. THE STATISTICS OF NUCLEAR SPINS 


With a simple and elementary method we shall now show 
that the nuclear spin [A-3] has an important influence on the 
statistics of rotational states (for example, in a diatomic gas). 

We consider a molecule X, that is composed of two 
identical atoms X. The eigenfunctions of such a rotator 
(‘dumbbell’) are known to be spherical harmonics (see 
Sec. 31), and the eigenvalues are 


h2 


gery 


i+ 1), 
where A is the moment of inertia of the molecule. For 
even 1, the spherical harmonics are even: 

Y(8, p) = (—1)' Y(8, o'); Y=a-F, punt. 
For zero nuclear spin we have Bose statistics, and only 
these symmetrical states occur. However, if the nuclear 


spin J is different from zero, the situation is different. 
A nuclear spin J has the 2J-+1 orientation possibilities m,, 


—-I<m,<4+1. 


Therefore, the molecule has (27-+1)? possible nuclear spin 
states m‘?, m?, which fall into three classes: 


OP - 62, 21 +1 states (1-41)(21 +1) 
Cx) Cnt Cal’ On mem; I(2I+1) states in all 


C8 O8~ OR Omt mim; I(21-+1) states . 


‘The separations of these rotational states (1! = 0,1, 2,...) are very small 
in comparison with the excitation energy of a molecular vibration. There- 
fore, to each vibrational level belongs an entire band of rotational levels, and 
this is the reason for the typical band spectra of molecules. 
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As before, we must correctly combine these spin functions 
with the space functions: 


a. I integral, bosons: (symmetric in spin + space) 


number of states with even! I+1 
~ number of states with odd l vi 


b. I half-integral, fermions: (antisymmetric in spin + space) 


number of states with even T if 
~ number of states with odd 1 J+1- 


Here Q is exactly the intensity ratio of adjacent lines 
(L=0, 1, 2,...) in a band and, therefore, it can be deduced 
from molecular spectra. From a knowledge of Q, the value 
of J and the symmetry character of the nucleus can be 
determined. 

Historically, this played a role in determining whether 
the atomic nucleus is composed of protons and electrons 
or of protons and neutrons.’ 


7™See W, Pau, “‘Zur alteren und neueren Geschichte des Neutrinos” in 
Aufsdtze und Vortrdge tiber Physik und Erkenntnistheorie (Vieweg, Braunsch- 
weig, 1961). 


Exercises 


40. FUNDAMENTAL SOLUTION FOR INTERVAL 


1. By the following two methods, determine the funda- 
mental solution for a particle that can move freely between 
two parallel walls which are separated by a distance L: 


a. By appropriate superposition of eigensolutions 


Yn(X, t) = Un(@) exp 


= 5B, 


taking the completeness relation 
2 ale )ux(x') = 6(a@ — a’) 


into account. 


b. Using the fundamental solution for a completely free 
particle and using the method of images to satisfy the 
boundary conditions. 

Show also that the expressions obtained with the two 
methods are identical. For this, use a property of the 
? function defined by 4 


8,(2| 1 aS exp [2niz] exp[intn?] 


n=—@ 


= | oe S cos 2nz-exp[inrn®], 


n=1 


1See E. T. WHITTAKER and G. N. WATSON, 4 Course of Modern Analysis 
(Cambridge University Press, New York, 1962), pp. 462 ff. 
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the needed property being 
ae ee 9 
0,(2|7) = (— tt) oe |= : ( 


14 


2. Using the fundamental solution K(a, x’, t) found in the 
previous part, determine the motion of a wave packet 


according to 
L 


p(x, t) =| da' f(x’) K(a, x’, t) , [40.1] 


assuming that the wave packet has the form f(«#) for t—0. 
Choose for f(x) the series obtained by reflection of the 
Gaussian distribution 


eee |) iii 
f(a — 4%) = (aoV 220) t-exp tl re + i (x va)mr , 


Show that [40.1] can also be represented by 


+00 
wa, t) = > {ys(2nL + (w — 2), t) 
—yo(2nL —(a—a,),t)}, [40.2] 
where y,(#—2, t) is the solution for a completely free 
particle which belongs to the initial distribution f. Using 


[40.2], discuss the interferences in the probability distribu- 
tion, P(x,t)=ypy*, by separating P into a classical term, 


i Dd {lyr l2+ | yal, 

and an interference term, 

P.= Dd lyayn + payaut voye + payin, 
where wy. is defined by 

y= Wo(2nL + (w — ap), t). 

Show that the remainder is negligible for times that are 
not too large (o,+ h*t?/4m’o3 <1’). In addition, represent 
P, with %, functions and show, by means of the property 


of the #, function given in the first part, that P. goes to 
the uniform distribution 1/Z for t->co [A-5]. 
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41. BOUND STATES AND TUNNEL EFFECT 


1. From the condition that the wave function, along with 
its first derivative, be continuous at the discontinuities of 
the potential, determine the energy eigenvalues (bound 
states): 


a. for a one-dimensional potential well (Fig. 41.1a); 

b. for a three-dimensional, spherically symmetric well 
(Fig. 41.1b), under the assumption that the wave function 
too is spherically symmetric (S state). 


Discuss the number of eigenvalues as a function of the 
dimensions of the well. 


Figure 41.1 


2. With the same condition as in the first part, calculate 
the transmission coefficient |w,/u,|? for a rectangular bar- 
rier as a function of the energy H# (>0). This is an ex- 
ample of the tunnel effect (Fig. 41.2). 


Figure 41.2 
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42. KRONIG-PENNEY POTENTIAL 


Derive the equation which determines the energy eigen- 
values E for a periodic rectangular potential (the Kronig- 
Penney potential). Use the condition that the wave func- 
tion, along with its first derivative, is continuous at the 


V(z) 
Ve 


0 b a z 
Figure 42.1 


discontinuities of the potential, as well as the following 

general property of the eigenfunctions of a periodic potential: 
p(x) =exp[tkx]u,(x) ; u,(x) periodic . 

Discuss the behavior of H(k) in the limiting case b->0, 

ab(mV,/h?) = P = constant. 

43. SPHERICAL HARMONICS 


The spherical function Y,(#, y), which yields a harmonic 
polynomial of degree | when it is multiplied by 7", satisfies 
the differential equation 


1 0) - OL, il sone: . a 
sin? (sino + aap Se +U0+1)¥,=0. 


It can be expanded in terms of the 21+1 linearly inde- 
pendent functions 


Vim, p) = P(x) em, 


where «=cos# and —l<m<-+H; the Y,,, satisfy the dif- 
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ferential equation 


(1 — x?)y"— 2ay’+ {u+ 1)— : hy =O. 


1— 2? 
There is only one solution for a given value of m? that 
is finite at e—+1 and «=—1. 

1. With the aid of the substitution 


y = (1—2?)m2y, 
and starting with the cases m=+1, show that the dif- 


ferential equation for P” is solved by 


(—1)! amt 
2H dame 


P(x) = (1— w?)m/2 (1— a)? [43.1] 
(the numerical factor here is the conventional one). Also, 
by applying this formula for both positive and negative m, 
prove the relation 
m —m (I + m)! 
1 = Oye ar ; where C= (—1)™ (I—m)! C 
The value of C,,, can be verified, for example, by com- 
paring with the coefficient of the highest power of x in [43.1]. 
2. Starting with the identity 


(1— #?)41 4 2(A+1){(24+1)(1— a2) 
—24(.1—a?)47}=0, [438.2] 


d2 
dav? 


where A is an arbitrary real number, derive the recursion 
formula 


pe Spe (211i at PP [43.3] 
by differentiating /+-m times. Using 


qQitm+1 5 
er ee ae) 
Qitm-1 


qit™ 1 ane 1 ) ——— (1— 2)! 
X 1 e are | a glia) aroma ’ 
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verify the additional identity 
a P+ (L+ m)V1— 02pm — Pm, =0, 
which can be brought into the form 


(20+ LaPP= (l—m +1) Pht (0+ m) Pes [43.4] 


by applying [43.3]. 
3. Calculate the normalization integral 


| {PP a)}2da = NP. 


For this purpose, derive a recursion formula from [43.4] 
which connects N7?,, with Ny. For the case m=I, there 
is an additional recursion formula that follows from [43.2]. 
The result is 

(l-+m)! +m)! 2 


(my! W4E1- 


[43.5] 


4. Using the normalized spherical harmonics, 
ih 


Fal, @ = ae —- —— P7(cos 8) eim?, 
Nj 


derive the matrix elements 


a 


(UU, m+1|e,+%e,|l,m) and (l’, mles|l, m) 
of the unit vector e whose components are 
é, + te, = sinde'? and e,=— cost. 


They are different from zero only when l’=1-+1 or l’/=I—1. 


44. FUNDAMENTAL SOLUTION FOR HARMONIC OSCILLATOR 


Determine the fundamental solution (where t= a) 


K(a, a’, t) = 2 )h,(v) exp[—i(n+4)t] [44.1] 
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of the harmonic oscillator by evaluating the sum over n 
in closed form. For this purpose it is convenient to use 
the representation 


+0 


al 
exp[— #2] = oa |emP [— v? + 2iva]dv, [44.2] 


-eo 


from which 


H, (2) = (= 1)" exp [a*] (s3) exp (ae? | 


= (—1)" exp[a?] oa [rier exp[— v?+ 2ive]dy [44.3] 


—@ 


follows for the Hermite polynomials; thus, 


1 et) it 
h, = ———==|—— Ke am = 
- Viale. 2 ex [+3 yz 
x feioy exp[— v?+ 2iva]dv. [44.4] 


mo (ee) 


Substitution of this representation into [44.1] yields a 
double integral for K(a#, x,t), in the integrand of which 
the summation over » can be carried out. The double 
integral can then be easily evaluated by means of [44.2]. 
The result is 


1 . (vw? + av'2) cost — 2a’ 
KH ap t) = ——————. exp Car 5 
V 2x1 sint 2 sint 


The completeness relation 


tea a’, 0) = zie )h, (a) = 6(4@ — 2’) 


also follows directly from the integral obtained for K. 
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45. ANGULAR MOMENTUM 


1. Show that the components of angular momentum can 
be represented in polar coordinates by 


ae a .cos? 9 
1 3 sin 9 dg)’ 
cos? 9 
— iP, = e-i# 
i ae (- er sin 3 2). 


In addition, show that application of these operators to 
Y, .(O, 9) = Pre? yields 


(an + iP) L palts ~) = Y aU, ¢) ) 
(P, — tP2) Yigal, p) = — (0+ m)(l — ne 1) ¥ ienl@, @) « 
The equation (see Sec. 43) 


(i+ m)! 


Eg 1) ea 


Pp" [45.1] 


can be used in proving the second formula. 
2. The matrix elements 


(l,m+1|P,+ iP, 


) and (l,m—1|P,—iP,\l, m) 


are to be calculated by using the value of the normaliza- 
tion integral | {P™(x)}2da (Sec. 43). Verify that P,—<P, 
is Hermitian conjugate to P,+iP,. 

3. By means of the substitution t+—(1—~2*)/t, verify 
[45.1] with the help of the complex integral representation 


ee l= (i+m)! {1 — (x 
Pe) S01 =a)" sat ot ae vai 


Here, it is essential that not only [+m be integral but 
that l be integral also. [Regularity of the integrand for 
t =1—a@ and t=—(1+2).] 
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46. PARTIAL WAVES 


1. If the potential V,(r) in the wave equation 
a 2m 
Pr (r0,) + k?v, — rE Vag te— 0; 


falls off faster than 1/r for large 7, then the asymptotic 
behavior of the wave function v, for large r is given by 


»,(r) ~~ sin (ir — 1 ; a 4.) [46.1] 


in the force-free case, (2m/h?)V,(r) =1(l+1)/r?, we have 
6,=0. Show that normalization according to 


posvwr dr = 6(k’— k) [46.2] 


o 
is equivalent to the requirement 
C=/2)x. [46.3] 
For this purpose use [46.2] in the form 


R  k+(Ak/2) 
jim Jrar)oovat'— at [46.4] 
0 k—(Ak/2) 


and use the relation 


R 


1 d da 
fore rar ae fre) (ros) — (ron) £ (ron) , [46.5] 


r=R 


which follows from the continuity equation. The limit 
R—->oco must be taken after the integration over k’. 

2. Show that this result is also valid in the case of the 
Coulomb potential if a correction of the order of logr is 
added to the phase in [46.1]. 
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3. For the wave function for the Coulomb potential in 
parabolic coordinates, caleulate the diagonal element of the 
coordinate <= 4$(A4,—A,), Eq. [18.66] (expressed in terms 
of the unit of length a). For doing this, the generating 
function for the Laguerre polynomials LZ” , ean be used, 


as can the normalization integral that was already eal- 
culated in Section 17 (see Eq. [17.42)]). 


47. THE SYMMETRICAL TOP 


Discuss the eigenvalue problem for the symmetrieal tep 
which has moments of inertia 4=B and C. The Hamil- 
tonian function of classical mechanies, expressed in terms 
of the Euler angles 3, y, y and their eanonieally conjugate 
momenta Py, Pz, Py, 18 


Po , (Pe—cosdp,)? py 
hoa? eee 
9A Asin*® | 2¢ 


This leads to the wave equation 


; even eee a if by an\* 
same ne e+ 2 — cos o = ) y 


2A sin? d VO) | 2A sin?d lg Na 
_i.,. 6, 
oe te we 


A trial solution 
y = 0(8) exp[i(my — m’g)] 


leads to the differential equation 


. Galea. OG 
sin? a (sino z 


— (m’? + m?+ 2 cos fmm’ — Asin? 3)6 == 0, ee 


where 
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Substituting #=(1—cos#)/2, 1—x= (1 +cos9)/2, we ob- 
tain 


d 
a(1—z) aa — x) “| 


1 
st i Ai — | (ie m')? + (A+ mm')e\9 0 [47.3] 

from [47.1]. With 
( — glmtm')/2 (1 — ap )om—m' 12 ee : [47.4] 


the equation for f is the hypergeometric differential equation 


a(1—x)f'+ {y— («+ B+ 1)a}f'—apBf=0, [47.5] 


with 
a+B=2m+1 
aB=—A+m(m+1) }. [47.6] 
y=m+m'+1 


If y is neither zero nor a negative integer, then an integral 
of [47.5] which is regular at 2=0 is given by the series 


a(a+1) BB +1), 
1-20 y(y +1) 


a 
F(a, B, y; amit oPa 4 aE nas 
Here, we allow m and m’ to have integral and half-integral 
values, of both signs, but such that m-+m’ is integral. 
Show that a solution valid for integral values of y of 
both signs is given by the integral 


F(a, B, y, %) 
__ exp[iaa] (1—«) ee Cee ae 
—— rome (2-1 (1 — t)? (1—tw)-Fdt. 
oO 


The path C encircles the points 0 and 1 in the positive 
sense, but the point 1/x is outside of C. 
We have 


F=Iy)F(, B,y, a) for y=1, 2,.... 
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The eigenvalues of 4 come from the requirement that F 
be a polynomial and [47.4] be finite, which means that a is 
a negative integer or zero. (Interchanging a and # does 
not yield a new solution.) Show that this is the case if 
and only if 


A= j(j +1) [47.7] 
(j is positive by definition), and 
j+m and j+m’ are integral and nonnegative. [47.8] 
We then have 
ea=—jtm, B=jtmt+1l, y=m4+m4+1. [47.9] 


Show, by means of the two substitutions t=1/s and 
t= (1/x)((~—s)/(1—s)) into the integral representation, that 
the function 6,,,,, in [47.4] with f, =F. satisfies 


m 


ymin pm") Fm)! 


(ee — a . . Conia’ fle 
nerd Gam yiG— my me? Taal 
and derive 
m____—*(j—m)! m+m' de —m!" m! 
Pra pw (Ge) Or aye) 
1 ; : 
a 


dintat : ; 
x (<:) {air (1 — x) }, [47.11] 
Moreover, show that 
(1— a) =P Py — B, y— a, y, @) = F(a, B, y, @) 
implies 
Oe =O . [47.12] 


Finally, calculate the normalization integral 


1 1 
W=[(0n.n(0)}*da =f (1— 2) (Pde. [47.13] 
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In this expression, the (Ff)? can, for example, be replaced 
by the product of the two expressions for F given in [47.11]; 
then a partial integration can be performed, and the Euler 
integral of the first kind, 


iT 


Ip) IQ) 
(a eS qld ee 0; 0 47.14 
Je ( #) £ Pip a) ; p> q> [ ] 


0 


can be used. 
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Appendix. Comments by the Editor 


[A-1] (p. 12). The problem raised in this paragraph con- 
stitutes one of the most fascinating aspects of the history 
of quantum theory. It is the question of whether the prob- 
abilistic nature of quantum theory is compatible with the 
requirement that it give a complete description of nature. 
This problem was heatedly discussed at the Fifth Solvay 
congress in 1927 from which the “Copenhagen interpreta- 
tion’? emerged victorious. While the latter is a credo in 
this compatibility, the opposite view defended by de Bro- 
glie, Einstein, Schrédinger, and later by Bohm is a con- 
viction in the deterministic foundation of a complete theory. 

The ‘intuitive pictures” to which Pauli refers in this 
paragraph make allusion to the “guiding wave” of de Bro- 
glie, which is a special ease of the “hidden variables” of 
the “deterministic scheme.’ The problem is most clearly 
exposed in Pauli’s article ‘‘Remarques sur le probleme des 
paramétres cachés dans la mécanique quantique et sur la 
théorie de l’onde pilote’’ in the volume dedicated to de 
Broglie’s sixtieth birthday, Louis de Broglie physicien et 
penseur (Albin Michel, Paris, 1953). He writes (p. 35) 
“J’eus le plaisir de discuter cette théorie au congrés Sol- 
vay 1927 avec L. de Broglie. Peu aprés, de Broglie aban- 
donna en faveur de l’interprétation complémentaire de la 
mécanique quantique de Heisenberg et Bohr; les raisons 
en sont exposées en détail dans son ‘Introduction 4 la mé- 
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canique ondulatoire’ (1929).” Later, on p. 37, Pauli pro- 
ceeds to discuss from a general point of view “‘la tentative 
de compléter la mécanique quantique de facon 4 en faire 
un schéma déterministe 4 aide de paramétres cachés; la 
théorie de l’onde pilote n’en est qu’un example spécial.”’ 
And he concludes (p. 42) ‘‘Ce sont ces raisons physiques, 
qui n’ont rien 4 faire avec des préjugés philosophiques 
quant 4 V’interprétation et & la justification de théories 
physiques en général, qui me font juger que l’interpré- 
tation de la mécanique quantique basée sur Vidée de la 
complémentarité est la seule admissible. Loin de consi- 
dérer comme définitif état actuel de la mécanique quan- 
tique dans le domaine relativiste, je crois pourtant que le 
développement de cette théorie ne fera que nous éloigner 
davantage de la possibilité d’une interprétation causale et 
déterministe.” 

As to Einstein’s belief mentioned in this paragraph it is 
most clearly expressed in his article ‘‘EKinleitende Bemer- 
kungen ueber Grundbegriffe’’ in the same volume. He 
writes (p. 6) “Es gibt so etwas wie den ‘realen Zustand 
eines physikalischen Systems, was unabhangig von jeder 
Beobachtung oder Messung objectiv existiert und mit den 
Ausdrucksmitteln der Physik im Prinzip beschrieben wer- 
den kann.”’ And p. 8, “Man fihlt sich daher gezwungen, 
die Beschreibung eines Systems durch eine Y-Funktion als 
eine unvollstandige Beschreibung des realen Zustandes zu 
betrachten.”’ 

For a more recent discussion of the problem of hidden 
variables see J. S. Bell, Rev. Mod. Phys. 88, 447 (1966). 
For a very recent experimental test, see S. J. Freedman 
and J. F. Clauser, Phys. Rev. Letters 28, 938 (1972). 


[A-2] (p. 28). The d(x) and all its derivatives are tempered 
distributions. 


A distribution is defined as a linear continuous func- 
tional F[p] of indefinitely differentiable functions g(#) with 
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bounded support. This means 
Flaigi+ A2Q2] = AF lgi)+ AF [2] ’ 
lim Flp,]= Fle), 


for any sequence 9,, 9, ... such that lim y;=q, and the 
nonzero values of g(x) are all contained in a finite do- 
main of &. 

A tempered distribution is defined with functions y(z) 
all derivatives yp” of which tend asymptotically strongly 
to zero, 


Jim, |e p(x) = 0 


for any nonnegative J and m. 

A locally integrable function f(x) which increases suffi- 
ciently slowly at infinity (ie., for which there exist posi- 
tive numbers A and « such that |f(a”)|<Aa* when |«|—> 00) 
defines a tempered distribution 


fry] = [f0) y(o) ae 


Two locally integrable functions f, and f, define the same 
distribution if /,—f,— 0 except on a set of measure zero. 
Therefore all representations of the Dirac function 6(x) of 
which examples are given on p. 29 define the same dis- 
tribution 


é[y] = / (x) p(”) da = »(0) . 


If a locally integrable function is differentiable, the deriv- 
ative of its distribution is defined by 


Fly] =|revyte da = — [fey da = —f[yp’]. 


In quantum mechanics the scalar product of an eigen- 
function f of an observable with continuous spectrum with 
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a general wave function y(2), | f(x) p(«) daw, defines a tem- 
pered distribution f[y] because any wave function y(x) is 
square integrable. 

For a concise exposition of the properties of distributions 
see, e.g., A. Messiah, Quantum Mechanics (North Holland, 
Amsterdam, 1961), vol. 1, appendix A. 


[A-3] (pp. 156, 175). Pauli’s role in the history of the 
spin both of the electron and of nuclei is most significant. 
As early as 1924 he proposed the hypothesis of a ‘‘resulting 
angular momentum” of nuclei. Indeed, in Naturwiss. 12, 
741 (1924) he discusses the “satellites of some spectral 
lines” in terms of what today is called the hyperfine inter- 
action; he writes ‘‘Diese Energieunterschiede werden dabei 
als vom zusammengesetzten Bau des Kernes herrthrend 
aufgefasst und es wird angenommen, dass der Kern im 
allgemeinen ein nicht verschwindendes resultierendes Im- 
pulsmoment besitzt.” 

It is remarkable that at the same time Pauli was not 
ready to accept the idea of an angular momentum of the 
electron, although in Z. Physik 31, 373 (1925) (submitted: 
2 December 1924) he drew attention to the ‘‘strange two- 
valuedness” of the quantum states of the electron. This 
conclusion was based both on the “doublet structure of 
the alkali spectra” (spin-orbit splitting) and on the ‘‘vio- 
lation of Larmor’s theorem” (anomalous Zeeman effect). 
He writes in the quoted paper ‘‘Die Dublettstruktur der 
Alkalispektren, sowie die Durchbrechung des Larmor- 
theorems kommt gemass diesem Standpunkt durch eine 
eigentiimliche, klassisch nicht beschreibbare Art von Zwei- 
deutigkeit der quantentheoretischen Eigenschaften des 
Leuchtelektrons zustande.’’ Pauli explained the reason for 
his hesitation to accept the idea of the electron spin in 
his 1946 Nobel Prize Lecture (reprinted in Collected Scien- 
tific Papers by Wolfgang Pauli, Interscience, New York, 
1964, vol. 2, p. 1080) as follows: “‘Although at first I strongly 
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doubted the correctness of this idea because of its classical 
mechanical character, I was finally converted to it by 
Thomas’ calculation [L. H. Thomas, Nature 117, 514 (1926) 
and Phil. Mag. 3, 1 (1927). Compare also J. Frenkel, Z. 
Phys. 87, 243 (1926)] on the magnitude of doublet split- 
CMS cick: 

It was indeed because of the discrepancy by this Thomas 
factor of two that Pauli did not believe in Kronig’s cal- 
culation of the spin-orbit splitting (unpublished, begin- 
ning 1925). More details on this subtle episode of the his- 
tory of science are found in the articles by R. Kronig and 
by B. L. van der Waerden in Theoretical Physics in the 
Twentieth Century. A Memorial Volume to Wolfgang Pauli 
(Interscience, New York, 1960). 

With regard to the history of the nuclear spin an amusing 
account can be found in the article by S. A. Goudsmit in 
Physics Today 14, No. 6, p. 18 (1961); he writes “For a 
number of years, whenever I met Pauli, he would remark 
cryptically that he ‘could afford not to be quoted.’ It 
was only in the late thirties that I found out to what he 
referred.” 


[A-4] (p. 174). This is true for spin-independent forces 
only. In this case the scattering does not give rise to spin 
flip so that identical particles with antiparallel spins remain 
distinguishable in the collision. 


[A-5] (p. 178). This problem is treated in great detail 
by M. Born in Dan. Mat. Fys. Medd. 30, No. 2 (1955). 
This paper was dedicated to Niels Bohr on the occasion 
of his seventieth birthday. 
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discrete energy spectrum of, 
95-97, 104-105 
energy eigenvalue of, 96, 104, 105 
principal quantum number of, 
96, 97 
radial quantum number of, 95 
separation of, in parabolic 
coordinates of, 102-106 
separation of, in spherical 
coordinates, 88-89 
Hyperfine interaction, 196 
Hypergeometric function, 
confluent, 81-85, 94-96, 107 
general, 82 
integral representation of, 82 


Incident plane wave, 107, 111, 120 
Initial value problem, 33-37, 177- 
178, 182-183 
Integral representation, 
cf confluent hypergeometric 
function, 83, 85, 987 
of cylinder functions, 114 
Euler, of I function, 81 
of Laguerre polynomials, 77 
of spherical harmonies, 184 
Interaction, 
Coulomb, 54, 169, 171 
between particles, 51 
spin-orbit, 169 
Interaction potential, 60 
Interference effects (phenomena), 
12, 14 
Interference of probabilities, 20, 
178 
Interference term, 174-175, 178 


Jeffreys, H., 132 
Jordan, P., 69, 153 


Kemble, E. C., 137 
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Klein, F., 160 

Kramers, H. A., 132 

Kronig, R., 197 
Kronig-Penney potential, 180 


Laguerre polynomials, 76, 84, 96, 
186 
Laplacian, 
in orthogonal curvilinear 
coordinates, 102 
in parabolic coordinates, 103 
in spherical coordinates, 88 
Legendre polynomials, 91, 92, 116, 
122 
Light, coherence property of, 18-19 
Light quanta, 1, 15 
Linear harmonic oscillator, 55 
Linear operator, 41 
Linear superposition of 
(stationary) states, 24, 27, 34 
Lorentz group, 164 


Magnetie field, 38, 42, 46, 47, 74, 
97, 159, 165 

Matrix, 

Hermitian, 63 

unitary, 140, 160 

Matrix element, 63, 138 

of angular momentum, 154 
of harmonic oscillator, 64 
Matrix mechanics 

Matrix method, 

for continuous spectra, 142 
of harmonic oscillator, 69-72 
Matrix representation, 68, 69 
Measurement, 


momentum, 16-17 

position, 14-15 

Messiah, A., 196 

Method of images, 177 
Microscope, 14—15 

Mixture, 21-22 

Moment of inertia, 175, 186 
Momentum conservation, 14 
Momentum measurement, 16—17 


Mott, N. F., 173 
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Neutron, 168, 176 
Nonrelativistic wave equation, 4, 
37 
Normalization, 5 
in the continuum, 27-31 
of Hermite polynomials, 58 
for plane harmonic oscillator, 79 
of spherical harmonics, 91, 182 
for spinning electron, 159 
for symmetrical top, 188 
Nuclear spin, 175-176, 196-197 
Nucleus (atomic), 93, 176 


Observable quantities, 43 
Operator, 
differential, 3, 153 
exchange, 165 
Hamiltonian, 39, 42 
Hermitian, 40, 41, 42-43 
linear, 41 
transformation, 139 
Orthogonality relation, 25, 26, 28, 
58 
for spin, 157 
Ortho helium, 171 
Orthonormal functions, 31 
Ortho states, 171 


Parabolic coordinates, 101, 186 
Para helium, 171 
Para states, 171 
Parseval, formula of, 5 
Partial waves, method of, 117-118, 
185-186 
Particle, 
in box, 23-27 
in magnetic field, 38, 39 
Particle current density, 111 
Pauli, W., 156, 164, 166, 168, 176, 
193-194, 196-197 
Pauli spin matrices, 158 
Periodic potential, 180 
Perturbation theory, 67, 68, 121, 
142 
in matrix representation, 143-147 
time-dependent, 147-151 
Phase correction, logarithmic, 110 
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Phase integral, 136 
Phase shift, 119-120 
Phase velocity, 2 
Polar coordinates, 

plane, 73 

spherical, 88 
Position measurement, 14-15 
Potential, 42, 54, 67 

central, 88, 90, 116—120 
Coulomb, 92, 97, 100, 102, 109- 
110, 11], 173 

periodic, 180 

range of, 122 

vector, 38 

Potential barrier, 137 

Potential well, 179 
Probabilistic nature of quantum 

theory, 193 

Probability, reduction of, 19, 21 
Probability current density, 25, 48 
Probability density, 11, 19-20, 24, 

52 

Probability distribution, 43, 178 
Proton, 168, 175, 176 

Pure case, 21, 22 


Quantum number, 
angular momentum, 96 
principal, 96 
radial, 95 


Radial differential equation, 92, 117 
Range of potential, 122 
Recoil momentum, 18 
Recursion formula, 
for Hermite polynomials, 60 
for Laguerre polynomials, 82 
for spherical harmonics, 181 
Reduced mass, 93, 173 
Relativistic particle mechanics, 1 
Relativistic spin theory, 164 
Relativistic wave equation 
(sealar), 3 
Representation, matrix, 68, 69 
Rotation group, 153 
Rotation in space, 160-164 
Rotator (“dumbbell”), 175 
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Russel-Saunders coupling, 169 
Rutherford scattering formula, 111 


Saddle point method, 128 
Scattered wave, 107, 111, 120, 121 
Scattering center, 118 
Scattering cross section, 110-111, 
118-119 
Scattering formula, Rutherford, 111 
Scattering of low-energy particles, 
123-125 
Schrodinger E., 39, 55, 79, 193 
Schrodinger equation, 39, 55, 68, 
88. See also Wave equation 
Secular equation, 147 
Selection rule for harmonic 
oscillator, 70 
Separable equation, 72, 73 
Sexl, Th., 111 
Similarity transformation, 141 
Singlet, 170, 171 
Singlet scattering amplitudes, 174 
Sommerfeld, A., 101, 105, 160 
Spectral analysis, 8 
Spherical harmonics, 89-91, 154, 
175, 180-182 
Spherical wave, 
incoming, i112 
outgoing, 108, 112 
Spin, 
history of, 197 
of electron, 156, 157 
Spin density, 162 
Spin eigenfunctions, 169 
Spin matrices, 158 
Spin operators, 156 
Spinor, 161 
Spin-orbit interaction, 169, 172 
Spin variable, 158 
Standing waves, 23, 55 
Stark effect (first order), 105 
State, 11, 21 
antisymmetric, 170 
bound, 179 
ground, 171] 
ortho, 171 
para, 17] 
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singlet, 170 

stationary, 23, 24 

symmetric, 170 

triplet, 170, 172 
Stationary state, 23, 24 
Statistics, 

Bose, 166 

classical and quantum, 19-22 
Fermi, 166 

of nuclear spins, 175-176 

of rotational states, 175 
Superposition principle, 11, 40 
Symmetrical top, 186-189 
Symmetric particles (bosons), 166 
Symmetric state, 170 
Symmetry classes, 165-167, 170 
Symmetry of Hamiltonian 

operator, 165 


0 function, 177 
Time-dependent forces, 46, 148 
Time-dependent perturbation 

theory, 147-151 
Total cross section, 119, 125 
Transformation, 

principal axis, 146 

similarity, 14] 

unitary, 140 
Transformation operator, 139 
Transformation theory, 138-143 
Transition probability, 150-151 
Translational key between classical 

and wave mechanics, 3, 38 

Transmission coefficient, 179 
Triplet, 170, 171 
Triplet scattering amplitudes, 174 
Tunnel effect, 137, 179 


Uhlenbeck, G. E. and Goudsmit, S.,156 


Uncertainty, 
of momentum, 15 
in position, 17 
Uncertainty relation, 8, 14, 19, 20 
Unitary matrix, 140, 160 
Unitary transformation, 140, 141 


Valence electron, 97 
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van der Waerden, B. L., 164, 197 Zeeman splitting, anomalous, 156 
Variance of Gaussian distribution, Zero-point energy, 57, 72 
9, 11 


Vector, in Hilbert space, 67 
Vector potential, 38 
Velocity, 

group, 2 

particle, 2 

phase, 2 
Virial theorem, 49 


Wave, 
incident plane, 107, 111, 120 
scattered, 107, 111, 120, 121 
standing, 23 

Wave equation, 
for central potential, 88, 116 
force-free, 112-114 
for harmonic oscillator in plane, 

v2 

for linear harmonic oscillator, 56 
nonrelativistic, 4, 37 
relativistic scalar, 3 
for symmetrical top, 186 
time-dependent, 52, 147 
time-independent, 39, 55, 132 
with uniform field, 126 

Wave function, 3, 34 
antisymmetric, 165, 171 
generalized, 51 
symmetric, 165 
two-component, 161 

Wave kinematics, 14 

Wave nature of atom, 19 

Wave normal, | 

Wave number, 10 

Wave packet, 3, 5, 11, 17, 178 
with Gaussian distribution, 8 
normalized, 5 
overlapping, 168 

Wave train, 16, 17 

Wave vector, 1 

Wentzel, G., 132 

Whittaker, E. T. and Watson, G. N., 
86, 177 

Wiener, E. P., 164 

WKB method, 132-137 
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“Wolfgang Pauli i was a great theoretical 
physicist who made many major con- 
‘tributions to the formulation of the. 
key concepts of present-day theoreti- 
= cal’ physics. It is therefore of consider- _ 
able interest to see his approach to the . 
presentation of these concepts through 
the various lecture courses he gave to — 
- students at the Swiss Federal Institute 
of Technology in Zurich. The books 
_ reviewed here are excellent transla- 
_ tions by S. Margulies and H. R. Lewis 
_ of lecture notes taken by different stu- 
dents and collaborators of Pauli. The 
notes have been carefully edited by. 

G P. Enz who provides a pistes and 
an appendix to each volume. .... A 
foreword to this series of back: by — 


Victor Weisskopf sets out the rationale 


for this publication of notes of lec- 
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never out of date,’ his style ‘is com- 
mensurate to the greatness of its sub- 
ject in its clarity and impact’ and his 
lectures ‘show how physical ideas can 


-be presented clearly and in good 


mathematical form, without being 

hidden in formalistic expertise.’ There 

can be no argument with these views. 
. There is no doubt that the books 


could provide excellent ‘additional 


reading’ and would be particularly 
valuable for the well-motivated and 
enthusiastic student. Equally they 
could provide a good basis for optional 
courses in theoretical physics either 


for physics- or mathematicsebased 


undergraduates. Together they are a 
remarkable testimony to the elegance 


of Pauli’s approach to physics,.his 


powers of exposition and his pene- 
trating insight. | should certainly wish 
to associate myself with the con- 
cluding sentence of Weisskopf’s fore- 
word—‘May these volumes serve as an 
example of how the concepts of 


_ theoretical physics were conceived and 


taught by one of the great men who 


created them.’ ''—Nature : 
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